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In recent papers, we and colleagues have introduced a way to visualize the full vacuum Riemann 
curvature tensor using frame- drag vortex lines and their vorticities, and tidal tendex lines and their 
tendicities. We have also introduced the concepts of horizon vortexes and tendexes and 3-D vortexes 
and tendexes (regions on or outside the horizon where vorticities or tendicities are large). In this 
paper, using these concepts, we discover a number of previously unknown features of quasinormal 
modes of Schwarzschild and Kerr black holes. These modes can be classified by a radial quantum 
number n, spheroidal harmonic orders (Z,ra), and parity, which can be electric [(—l) 1 ] or magnetic 
[( — Among our discoveries are these: (i) There is a near duality between modes of the 
same (n, Z,m): a duality in which the tendex and vortex structures of electric-parity modes are 
interchanged with the vortex and tendex structures (respectively) of magnetic-parity modes, (ii) 
This near duality is perfect for the modes' complex eigenfrequencies (which are well known to be 
identical) and perfect on the horizon; it is slightly broken in the equatorial plane of a non-spinning 
hole, and the breaking becomes greater out of the equatorial plane, and greater as the hole is 
spun up; but even out of the plane for fast-spinning holes, the duality is surprisingly good, (iii) 
Electric-parity modes can be regarded as generated by 3-D tendexes that stick radially out of the 
horizon. As these "longitudinal," near-zone tendexes rotate or oscillate, they generate longitudinal- 
transverse near-zone vortexes and tendexes, and outgoing and ingoing gravitational waves. The 
ingoing waves act back on the longitudinal tendexes, driving them to slide off the horizon, which 
results in decay of the mode's strength, (iv) By duality, magnetic-parity modes are driven in this 
same manner by longitudinal, near-zone vortexes that stick out of the horizon, (v) When visualized, 
the 3-D vortexes and tendexes of a (Z,m) = (2,2) mode, and also a (2,1) mode, spiral outward 
and backward like water from a whirling sprinkler, becoming outgoing gravitational waves. By 
contrast, a (2,2) mode superposed on a (2, —2) mode, has oscillating horizon vortexes or tendexes 
that eject 3-dimensional vortexes and tendexes, which propagate outward becoming gravitational 
waves; and so does a (2,0) mode, (vi) For magnetic-parity modes of a Schwarzschild black hole, 
the perturbative frame-drag field, and hence also the perturbative vortexes and vortex lines, are 
strictly gauge invariant (unaffected by infinitesimal magnetic-parity changes of time slicing and 
spatial coordinates), (vii) We have computed the vortex and tendex structures of electric-parity 
modes of Schwarzschild in two very different gauges and find essentially no discernible differences in 
their pictorial visualizations, (viii) We have compared the vortex lines, from a numerical-relativity 
simulation of a black hole binary in its final ringdown stage, with the vortex lines of a (2,2) electric- 
parity mode of a Kerr black hole with the same spin (a/M — 0.945) and find remarkably good 
agreement. 



PACS numbers: 04.25.dg, 04.25.Nx, 04.30.-w 

I. MOTIVATIONS, FOUNDATIONS AND 
OVERVIEW 

A. Motivations 

This is the third in a series of papers that introduce a 
new set of tools for visualizing the Weyl curvature tensor 
(which, in vacuum, is the same as the Riemann tensor), 
and that develop, explore, and exploit these tools. 

We gave a brief overview of these new tools and their 
applications in an initial Physical Review Letter [Tj . Our 
principal motivation for these tools was described in that 
Letter, and in greater detail in Sec. I of our first long, 



pedagogical paper [2! (Paper I). In brief: We are moti- 
vated by the quest to understand the nonlinear dynamics 
of curved spacetime (what John Wheeler has called ge- 
ometrodynamics) . 

The most promising venue, today, for probing ge- 
ometrodynamics is numerical simulations of the collisions 
and mergers of binary black holes [3] . Our new tools pro- 
vide powerful ways to visualize the results of those simu- 
lations. As a byproduct, our visualizations may motivate 
new ways to compute the gravitational waveforms emit- 
ted in black-hole mergers — waveforms that are needed as 
templates in LIGO's searches for and interpretation of 
those waves. 

We will apply our tools to black-hole binaries in Paper 
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IV of this series. But first, in Papers I— III, we are apply- 
ing our tools to analytically understood spacetimes, with 
two goals: (i) to gain intuition into the relationships be- 
tween our tools' visual pictures of the vacuum Riemann 
tensor and the analytics, and (ii) to gain substantial new 
insights into phenomena that were long thought to be 
well understood. Specifically, in Paper I [2], after intro- 
ducing our tools, we applied them to weak-gravity situ- 
ations ("linearized theory"); in Paper II [4], we applied 
them to stationary (Schwarzschild and Kerr) black holes; 
and here in Paper III we will apply them to weak pertur- 
bations (quasinormal modes) of stationary black holes. 



B. Our new tools, in brief 

In this section, we briefly summarize our new tools. 
For details, see Sees. II, III, and IV of Paper I [2 , and 
Sees. II and III of Paper II 0]. 

When spacetime is foliated by a family of spacelike 
hypersurfaces (surfaces on which some time function t is 
constant), the electromagnetic field tensor F^ v splits up 
into an electric field E% = and a magnetic field = 
\tij%Fjy which are 3- vector fields living in the spacelike 
hypersurfaces. Here the indices are components in proper 
reference frames (orthonormal tetrads) of observers who 
move orthogonally to the hypersurfaces, and e^t is the 
Levi-Civita tensor in those hypersurfaces. 

Similarly, the Weyl (and vacuum Riemann) tensor 
C^Xp splits up into: (i) a tidal field £y = C^q, which 
produces the tidal gravitational accelerations that ap- 
pear, e.g., in the equation of geodesic deviation, A a-- = 

—Sj^Ax k [Eq. (3.3) of Paper I]; and (ii) a frame-drag 
field By = ^ipqCjQpqi which produces differential frame- 
dragging (differential precession of gyroscopes), AQj — 

By k Ax % [Eq. (3.11) of Paper I]. 

We visualize the tidal field £ by the integral curves 
of its three eigenvector fields, which we call tendex lines, 
and also by the eigenvalue of each tendex line, which we 
call the tendicity of the line and we depict using colors. 
Similarly, we visualize the frame-drag field B by frame- 
drag vortex lines (integral curves of its three eigenvec- 
tor fields) and their vorticities (eigenvalues, color coded). 
See Figs. [2] and [3] below for examples. Tendex and vor- 
tex lines are analogs of electric and magnetic field lines. 
Whereas through each point in space there pass just one 
electric and one magnetic field line, through each point 
pass three orthogonal tendex lines and three orthogonal 
vortex lines, which identify the three principal axes of £ 
and B. 

A person whose body is oriented along a tendex line 
gets stretched or squeezed with a relative head-to-foot 
gravitational acceleration that is equal to the person's 
height times the line's tendicity (depicted blue [dark gray] 
in our figures for squeezing [positive tendicity] and red 
[light gray] for stretching [negative tendicity]). Similarly, 
if the person's body is oriented along a vortex line, a gyro- 



scope at her feet precesses around her body axis, relative 
to inertial frames at her head, with an angular velocity 
equal to her height times the line's vorticity (depicted 
blue [dark gray] for clockwise precession [positive vor- 
ticity] and red [light gray] for counterclockwise [negative 
vorticity] ) . 

We color code the horizon of a black hole by the 
normal- normal component of the tidal field, £nn, to 
which we give the name horizon tendicity, and also by the 
normal- normal component of the frame-drag field, Sat at, 
the horizon vorticity; see, e.g., Fig.|9]below. These quan- 
tities are boost-invariant along the normal direction N to 
the horizon in the foliation's hypersurfaces. 

A person hanging radially above the horizon or falling 
into it experiences head-to foot squeezing (relative accel- 
eration) equal to the horizon tendicity times the person's 
height, and a differential head-to-foot precession of gy- 
roscopes around the person's body axis with an angular 
velocity equal to the horizon vorticity times the person's 
height. 

For any black hole, static or dynamic, the horizon ten- 
dicity £nn and vorticity Bnn are related to the horizon's 
Newman- Penrose Weyl scalar ^2, and its scalar intrinsic 
curvature 1Z and scalar extrinsic curvature X by 

£ NN + iB NN = 2^ 2 = ~(n + iX) + 2(w-\(T) ; (1.1) 

[5], and Sec. Ill of [4]. Here p, a, /i, A are spin coef- 
ficients related to the expansion and shear of the null 
vectors I and n used in the Newman-Penrose formalism 
[with (I + n)/y/2 = u the normal to the foliation's hy- 
persurfaces, {I — n)/y/2 = N the normal to the horizon 
in the foliation's hypersurfaces, and e*2 = (m + m*)/y/2 
and £3 = (m — m*)/(iy/2) tangent to the instantaneous 
horizon in the foliation's hypersurfaces]. For station- 
ary black holes, p and a vanish, and £nn — ~\R> and 
Bnn = 

For perturbations of Schwarzschild black holes, it is 
possible to adjust the slicing at first order in the per- 
turbation, and adjust the associated null tetrad, so as 
to make the spin coefficient terms in Eq. ( |1.1| ) vanish 
at first order in the perturbation; whence £nn — 
and Bnn — — \X> For perturbations of the Kerr space- 
time, however, this is not possible. See App. [E] for de- 
tails. Following a calculation by Hartle [6 , we show in 
this appendix that for Kerr one can achieve 7Z + iX = 
—4(^2 + A^cr^ 1 )) on the horizon, accurate through first 
order. Here, and throughout this paper, the superscripts 
M (or subscripts (^) indicate orders in the perturbation. 

For the dynamical black holes described in pQ and 
for the weakly perturbed ho les i n this paper, we found 
that the spin terms in Eq. (1.1) are numerically small 



compared to the other terms, so £nn — and 
Bnn — — \X- I n addition, in a recent study of the ten- 
dexes and vortexes of approximate black hole initial data, 
Dennison and Baumgarte [7 found that these spin terms 
vanish to a high order in the small velocities of their 
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black holes, giving further evidence that these terms are 
typically negligible. 

Because X is the 2-dimensional curl of a 2-dimensional 
vector (the Hajicek field) [8 , its integral over the black 
hole's 2D horizon vanishes; and by virtue of the Gauss- 
Bonnet theorem, the horizon integral of 7Z is equal to 8tt. 
Correspondingly, for fully dynamical black holes as well 
as weakly perturbed black holes, the horizon integrals of 
£nn and Bnn have the approximate values pQ 



[ B NN ~0, [ 



NN 



-47T 



(1.2) 



C. Overview of this paper's results 

1. Slicing, coordinates and gauges 

Throughout this paper, we use slices of constant Kerr- 
Schild time t (which penetrate smoothly through the 
horizon) to decompose the Weyl tensor into its tidal and 
frame-drag fields; and we express our quasinormal per- 
turbations, on the slices of co nst ant t, in Kerr-Schild spa- 
tial coordinates (Sees. |II A| and |IIB| and also Paper II 
[4]). In the zero- spin (Schwarzschild) limit, the Kerr- 
Schild slices become slices of constant ingoing Eddington- 
Finkelstein time t and the spatial coordinates become 
those of Schwarzschild. Our choice of Kerr-Schild is dic- 
tated by these coordinates' resemblance to the coordi- 
nates that are typically used in numerical-relativity sim- 
ulations of binary black holes, at late times, when the 
merged hole is settling down into its final Kerr-black-hole 
state; see, e.g., Fig. [l5| below. 

For a perturbed black hole, the slices and coordinates 
get modified at perturbative order in ways that depend 
on the gauge used to describe the perturbations (i.e., the 
slicing and spatial coordinates at perturbative order) ; see 
Sec. ECl 

For spinning black holes, we perform all our computa- 
tions in ingoing radiation gauge (Sec. II A| and App. C|. 
For non-spinning (Schwarzschild) black holes, we explore 
gauge dependence by working with two gauges that ap- 
pear to be quite different: ingoing radiation gauge (App. 
[C| , and Regge- Wheeler gauge (App.|A|. Remarkably, for 
each mode we have explored, the field-line visualizations 
that we have carried out in these two gauges look nearly 
the same to the human eye; visually we see little gauge 
dependence. We discuss this and the differences in the 



gauges, in considerable detail, in Sec. |II C| and App. |D) 

For a Schwarzschild black hole, we have explored 
somewhat generally the influence of perturbative slicing 
changes and perturbative coordinate changes on the tidal 
and frame-drag fields, and on their tend ex a nd vortex 
lines, and tendicities and vorticities (Sec. II C). We find 
that the tendicities and vorticities are less affected by 
perturbative slicing changes, than the shapes of the ten- 
dex and vortex lines. We also find that while coordinate 
changes affect the shapes of the tendex and vortex lines, 



the tendicity and vorticity along a line is unchanged, and 
that in the wave zone a perturbative change in coordi- 
nates affects the tendicity and vorticity at a higher order 
than the effect of gravitational radiation. 

For this reason, in this paper we pay considerable at- 
tention to vorticity and tendicity contours, as well as to 
the shapes of vortex and tendex lines. 



2. Classification of quasinormal modes 

As is well known, the quasinormal- mode, complex 
eigenfrequencies of Schwarzschild and Kerr black holes 
can be characterized by three integers: a poloidal quan- 
tum number / = 2,3, an azimuthal quantum number 
m = —I, —I + 1, ... + i, and a radial quantum number 
n. For each {n, /,m} and its eigenfrequency uo n i mi there 
are actually two different quasinormal modes (a two-fold 
degeneracy). Of course, any linear combination of these 
two modes is also a mode. We focus on those linear com- 
binations of modes that have definite parity (App. |C|). 

We define a tensor field to have positive parity if it 
is unchanged under reflections through the origin, and 
negative parity if it changes sign. A quasinormal mode of 
order (n,/,m) is said to have electric parity [or magnetic 
parity) if the parity of its metric perturbation is (— I) 1 [or 
(— The parity of the tidal-field perturbation is the 
same as that of the metric perturbation, but that of the 
frame-drag field is opposite. In much of the literature the 
phrase "even parity" is used in place of "electric parity" , 
and "odd parity" in place of "magnetic parity" ; we avoid 
those phrases because of possible confusion with positive 
parity and negative parity. 

In this paper, we focus primarily on the most slowly 
damped (n = 0) quadrupolar (/ = 2) modes, for vari- 
ous azimuthal quantum numbers m and for electric- and 
magnetic-parity. Since we discuss exclusively the n = 
modes, we will suppress the n index and abbreviate mode 
numbers as (/, m). 



3. The duality of magnetic-parity and electric-parity modes 

In vacuum, the exact Bianchi identities for the Rie- 
mann tensor become, under a slicing-induced split of 
spacetime into space plus time, a set of Maxwell-like 
equations for the exact tidal field and frame-drag field 
[Eqs. (2.15) of Paper I [2 in a local Lorentz frame; Eqs. 
(2.13) and (2.4) of Paper I in general]. These Maxwell- 
like equations exhibit an exact duality: If one takes any 
solution to them and transforms £ — > B, B — >> — £, they 
continue to be satisfied (Sec. II B 1 of Paper I [2]). 

This duality, however, is broken by the spacetime ge- 
ometry of a stationary black hole. A Schwarzschild black 
hole has a monopolar tidal field S and vanishing frame- 
drag field B; and a Kerr black hole has a monopolar 
component to its tidal field (as defined by a spherical- 
harmonic analysis at large radii or at the horizon), but 
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only dipolar and higher-order components to its frame- 
drag field. 

When a Schwarzschild or Kerr black hole is perturbed, 
there is a near duality between its electric-parity mode 
and its magnetic-parity mode of the same (Z,m); but 
the duality is not exact. The unperturbed hole's dual- 
ity breaking induces (surprisingly weak) duality-breaking 
imprints in the quasinormal modes. We explore this du- 
ali ty breaking in con siderab le detail in this p aper (Sees. 
Il| |mAl |IIIB3[ and |III C 2[ and Apps. [C] and [E2| . 

If one tries to see the duality between electric-parity 
and magnetic-parity modes, visually, in pictures of the 
perturbed hole's tendex and vortex lines, the duality is 
hidden by the dominant background tidal field and (for 
a spinning hole) the background frame-drag field. To 
see the duality clearly, we must draw pictures of tendex 
and vortex lines for the perturbative parts SS and SB 
of the tidal and frame-drag fields, with the unperturbed 
fields subtracted off. We draw many such pictures in this 
paper. 

We have made extensive comparisons of the least 
damped (n = 0) electric-parity and magnetic-parity 
modes with (I = 2, m = 2). These two (2,2) modes 
(for any chosen black-hole mass M and spin parameter 
a) have identically the same complex eigenfrequency, i.e., 
they are degenerate (as has long been known and as we 
discussed above). This frequency degeneracy is an un- 
broken duality. 

Pictures of the perturbative vortex and tendex lines 
and their color-coded vorticities and tendicities show a 
strong but not perfect duality: For a non-spinning hole, 
the perturbative vortex lines and their vorticities for the 
magnetic-parity mode (e.g., Fig.[2| look almost the same 
as the perturbative tendex lines and their tendicities for 
the electric-parity mode (Fig. [l2|); and similarly for the 
other pair of lines and eigenvalues. As the hole's spin is 
increased, the duality becomes weaker (the corresponding 
field lines and eigenvalues begin to differ noticeably) ; but 
even for very high spins, the duality is strikingly strong; 
see bottom row of Fig. [12] below. The duality remains 
perfect on the horizon in ingoin g radia tion gauge for any 
spin, no matter how fast (Sec. Ill A and App. py), and 



there is a sense in which it also remains perfect on the 
horizon of Schwar zschi ld in Regge- Wheeler gauge (last 
paragraph of App. A 5). 



Digression: Electromagnetic perturbations of a 
Schwarzschild black hole 



As a prelude to discussing the physical character of 
the gravitational modes of a black hole, we shall discuss 
electromagnetic (EM) modes, i.e., quasinormal modes of 
the EM field around a black hole. The properties of 
EM modes that we shall describe can be derived from 
Maxwell's equations in the Schwarzschild and Kerr space- 
times, but we shall not give the derivations. 

Because the unperturbed hole has no EM field and the 




FIG. 1: (color online), (a) Some magnetic field lines in the 
equatorial plane for the (1,1) quasinormal mode of the elec- 
tromagnetic field around a Schwarzschild black hole, with 
Eddington-Finklestein slicing. The horizon is color coded 
by the sign of the normal component of the magnetic field. 
The configuration rotates counterclockwise in time, (b) Some 
magnetic field lines for this same quasinormal mode, in 3 di- 
mensions. 



vacuum Maxwell equations exhibit a perfect duality (they 
are unchanged when E — » B and B — >> —E), the EM 
modes exhibit perfect duality. For any magnetic-parity 
EM mode, the magnetic field pierces the horizon, so its 
normal component Bn is nonzero, while En vanishes. By 
duality, an electric-parity EM mode must have En ^ 
and Bn = 0. For a magnetic-parity mode, the near- 
zone magnetic fields that stick out of the horizon can 
be thought of as the source of the mode's emitted EM 
waves. We make this claim more precise by focusing on 
the fundamental (n = 0), magnetic-parity, I = 1, m = 1 
mode: 

Figure [2] shows magnetic field lines for this (1,1) mode, 
on the left (a) in the hole's equatorial plane, and on the 
right (b) in 3 dimensions with the equatorial plane hor- 
izontal. On the left, we see a bundle of magnetic field 
lines that thread through the horizon and rotate coun- 
terclockwise. As they rotate, the field lines spiral out- 
ward and backward, like water streams from a whirling 
sprinkler, becoming the magnetic-field component of an 
outgoing electromagnetic wave. The electric field lines 
for this mode (not shown) are closed circles that repre- 
sent the electric part of electromagnetic waves traveling 
outward at radii r ^> 2M and inward at radii r ~ 2M. 
This mode's waves, we claim, are generated by the near- 
zone, rotating magnetic field lines that thread the hole 
(Fig. An analogy will make this clear. 

Consider a rotating (angular velocity <r), perfectly 
conducting sphere in which is anchored a magnetic 
field with the same dipolar normal component Bn oc 
9?[F n (<9, cj))e- iat ] as the horizon's B N for the (1,1) quasi- 
normal mode (the red [light gray] and blue [dark gray] 
coloring on the horizon in Fig. [T]). At some initial mo- 
ment of time, lay down outside the conducting sphere, a 
magnetic-field configuration that (i) has this B n at the 
sphere, (ii) satisfies the constraint equation V • B = 0, 
(hi) resembles the field of Fig. [I] in the near zone, i.e., at 
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r < A = c/cr and at larger radii has some arbitrary form 
that is unimportant; and (iv) (for simplicity) specify a 
vanishing initial electric field. Evolve these initial fields 
forward in time using the dynamical Maxwell equations. 
It should be obvious that the near-zone, rotating mag- 
netic field will not change much. However, as it rotates, 
via Maxwell's dynamical equations it will generate an 
electric field, and those two fields, interacting, will give 
rise to the outgoing electromagnetic waves of a / = 1, 
m = 1 magnetic dipole. Clearly, the ultimate source of 
the waves is the rotating, near-zone magnetic field that 
is anchored in the sphere. (Alternatively, one can regard 
the ultimate source as the electric currents in the sphere, 
that maintain the near-zone magnetic field.) 

Now return to the magnetized black hole of Fig. [I] 
and pose a similar evolutionary scenario: At some initial 
moment of time, lay down a magnetic-field configuration 
that (i) has the same normal component at the horizon 
as the (1,1) mode, (ii) satisfies the constraint equation 
V • B = 0, and (iii) resembles the field of Fig. [I] in the 
near zone. In this case, the field is not firmly anchored 
in the central body (the black hole), so we must also 
specify its time derivative to make sure it is rotating at 
the same rate as the (1,1) quasinormal mode. This means 
(by a dynamical Maxwell equation) that we will also be 
giving a nonvanishing electric field that resembles, in the 
near zone, that of the (1,1) mode and in particular does 
not thread the horizon. Now evolve this configuration 
forward in time. It will settle down, rather quickly, into 
the (1,1) mode, with outgoing waves in the wave zone, 
and ingoing waves at the horizon. This is because the 
(1,1) mode is the most slowly damped quasinormal mode 
that has significant overlap with the initial data. 

As for the electrically conducting, magnetized sphere, 
so also here, the emitted waves are produced by the 
rotation of the near-zone magnetic field. But here, by 
contrast with there, the emitted waves act back on the 
near-zone magnetic field, causing the field lines to gradu- 
ally slide off the horizon, resulting in a decay of the field 
strength at a rate given by the imaginary part of the 
mode's complex frequency. 

This back-action can be understood in greater depth 
by splitting the electric and magnetic fields, near the 
horizon, into their longitudinal (radial) and transverse 
pieces. The longitudinal magnetic field is Bjy and it ex- 
tends radially outward for a short distance; the tangential 
magnetic field is a 2-vector B T parallel to the horizon; 
and similarly for the electric field, which has En = 
and so is purely transverse. The tangential fields actu- 
ally only look like ingoing waves to observers who, like the 
horizon, move outward at (almost) the speed of light: the 
observers of a Schwarzschild time slicing. As one learns 
in the Membrane Paradigm for black holes (Sees. III.B.4 
and III.C.2 of [9]), such observers can map all the physics 
of the event horizon onto a stretched horiz on — a spac elike 
2-surface of constant lapse function a = y/l — 2M/r <C 1 
very close in spacetime to the event horizon. On the 
stretched horizon, these observers see E T = N x B T 



(ingoing- wave condition), and the tangential magnetic 
field acts back on the longitudinal field via 



dB N 
dt 



V • (aB T ) = . 



(1.3) 



Here ^ V • (aB T ) is the 2-dimensional divergence in the 
stretched horizon, and the lapse function in this equation 
compensates for the fact that the Schwarzschild observers 
see a tangential field that diverges as 1/a near the hori- 
zon, due to their approach to the speed of light. 

Equation (1.3) is a conservation law for magnetic field 



lines on the stretched horizon. The density (number per 
unit area) of field lines crossing the stretched horizon is 
5 at, up to a multiplicative constant; the flux of field lines 
(number moving through unit lenth of some line in the 
stretched horizon per unit time) is B T , up to the same 
multiplicative constant; and Eq. ( |1.3| ) says that the time 
derivative of the density plus the divergence of the flux 
vanishes: the standard form for a conservation law. 

The longitudinal magnetic field a/- (#,</>) is laid down 
as an initial condition (satisfying the magnetic constraint 
condition). As it rotates, it generates the ingoing-wave 
near-horizon transverse fields embodied in E T and B T 
(and also the outgoing electromagnetic wave s far from 



the hole); and the divergence of aB T , via Eq. (1.3), then 
acts back on the longitudinal field that produced it, push- 
ing the field lines away from the centers of the blue ( dark 
gray) and red (light gray) spots on the stretched horizon 
toward the white ring. Upon reaching the white ring, each 
field line in the red region attaches onto a field line from 
the blue region and slips out of the horizon. Presumably, 
the field line then travels outward away from the black 
hole and soon becomes part of the outgoing gravitational 
waves. The gradual loss of field lines in this way is re- 
sponsible for the mode's exponential decay. 



5. The physical character of magnetic-parity and 
electric-parity modes 

For a Schwarzschild black hole, the physical character 
of the gravitational modes is very similar to that of the 
electromagnetic modes: 

Just as a magnetic-parity EM mode has nonzero Bn 
and vanishing E^r, so similarly: for a Schwarzschild 
black hole, the magnetic-parity modes of any (/,ra) 
have nonzero (solely perturbative) horizon vorticity 
SBnn = &nn, ctnd vanishing perturbative horizon ten- 
dicity SEnn — 0; and correspondingly, from the horizon 
there emerge nearly normal vortex lines that are fully per- 
turbative and no nearly normal, perturbative tendex lines. 

Just as in the EM case, the near-zone magnetic fields 
that emerge from the horizon are the source of the emit- 
ted electromagntic waves, so also in the gravitational 
case, for a magnetic-parity mode the emerging, near- 
zone, vortex lines and their vorticities can be thought 
of as the source of the emitted magnetic-parity gravita- 
tional waves (see the next subsection). In this sense, 
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magnetic-parity modes can be thought of as fundamen- 
tally frame-drag in their physical origin. Figure [2] below 
depicts a (2,2) example. We will discuss this example in 

Sec. ens 

For a Schwarzschild black hole, the electric- parity 

modes of any (Z,ra) have nonzero perturbative horizon 
tendicity 5£nn ^ 0, and vanishing horizon vorticity 
SBnn = = 0; and correspondingly, from the hori- 

zon there emerge nearly normal perturbative tendex lines 
and no nearly normal vortex lines. The emerging, near- 
zone, perturbative tendex lines can be thought of as the 
source of the mode's emitted electric-parity gravitational 
waves. In this sense, electric-parity modes can be thought 
of as fundamentally tidal in their physical origin. 

There is a close analogy, here, to the tidal and frame- 
drag fields of dynamical multipoles in linearized theory 
(Paper I [2 j): Electric-parity (mass) multipoles have a 
tidal field that rises more rapidly, as one approaches the 
origin, than the frame-drag field, so these electric-parity 
multipoles are fundamentally tidal in physical origin. By 
contrast, for magnetic-parity (current) multipoles it is 
the frame-drag field that grows most rapidly as one ap- 
proaches the origin, so they are fundamentally frame- 
drag in physical origin. 

When a black hole is spun up, the horizon vortici- 
ties of its electric-parity modes become nonzero, and the 
horizon tendicities of its magnetic-parity modes acquire 
nonzero perturbations. However, these spin-induced ef- 
fects leave the modes still predominantly tidal near the 
horizon for electric-parity modes, and predominantly 
frame-drag near the horizon for magnetic-parity modes 

(Sec, [m]). 



6. The (2, 2) magnetic-parity mode of a Schwarzschild hole 

In this and the next several subsections, we summa- 
rize much of what we have learned about specific n = 0, 
I = 2 modes (the least-damped quadrupolar modes), for 
various m. We shall focus primarily on magnetic-parity 
modes, since at the level of this discussion the proper- 
ties of electric-parity modes are the same, after a duality 
transformation 5S — » 5B, SB — >• —SS. 

It is a remarkable fact that, for a magnetic-parity mode 
of a Schwarzschild black hole, all gauges share the same 
slicing, and the mode's frame-drag field is unaffected by 
perturbative changes of spatial coordinates; therefore, 
the frame-drag field is fully gauge invariant. See Sec. 



II C| This means that Figs. [2}j8] are fully gauge invariant. 

We begin with the (2, 2) magnetic-parity mode of a 
Schwarzschild black hole. Figure [2] depicts the negative- 
vorticity vortex lines (red) and contours of their vorticity 
(white and purple [dark gray]), in the hole's equatorial 
plane. Orthogonal to the red (solid) vortex lines (but not 
shown) are positive-vorticity, vortex lines that also lie in 
the equatorial plane. Vortex lines of the third family 
pass orthogonally through the equatorial plane. The en- 
tire configuration rotates counterclockwise, as indicated 




FIG. 2: (color online). Some vortex lines (solid, red lines) 
and contours of vorticity (shaded regions) in the equatorial 
plane for the (2, 2) magnetic-parity quasinormal mode of a 
non-rotating, Schwarzschild black hole, with complex eigen- 
frequency u = (0.37367 - 0.08896i)/M where M is the hole's 
mass. The horizon (central circle) is color coded by the hori- 
zon vorticity Bnn as seen by someone looking down on the 
black hole; this vorticity is entirely perturbative. The thick, 
solid red curves are one set of vortex lines in the equatorial 
plane — the set with negative vorticity. These lines include 
some that emerge from the horizon in the negative-vorticity 
(red) regions, and some that never reach the horizon. The 
other, positive-vorticity, equatorial vortex lines are orthog- 
onal to the ones shown, and are identical to those shown 
but rotated through 90 degrees around the hole so some of 
them emerge from the horizon in the positive-vorticity (blue) 
regions. The contours represent the vorticity of the red 
(negative-vorticity) vortex lines, with largest magnitude of 
vorticity white and smallest purple (dark gray); the contours 
mark where the vorticity has fallen to 50%, 25%, 10%, and 5% 
of the maximum value attained at the center of the horizon 
vortex. The two dotted circles are drawn at Schwarzschild 
radii r = A and r = nX = A/2. They mark the approximate 
outer edge of the near zone and the approximate inner edge 
of the wave zone. The arrow marks the direction of rotation 
of the perturbation. 



by the thick dashed arrow. The dotted lines, at radii 
r = A and r = ttX (where A is the emitted waves' re- 
duced wavelength), mark the approximate outer edge of 
the near zone, and the approximate inner edge of the 
wave zone. 

Just as the near-zone electromagnetic (1,1) pertur- 
bations are dominated by radial field lines that thread 
the black hole and have a dipolar distribution of field 
strength, so here the near-zone gravitational perturba- 
tions are dominated by (i) the radial vortex lines that 
thread the hole and have a quadrupolar distribution of 
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their horizon vorticity Bnn = SBnn, and also by (ii) 
a transverse, isotropic frame-drag field B^ = B^ = 

— \$Bnn that is tied to Bnn in such a way as to guar- 
antee that this dominant part of B is traceless. 

This full structure, the normal-normal field and its 
accompanying isotropic transverse field, makes up the 
longitudinal, nonradiative frame-drag field B L near the 
horizon. (As we shall discuss below, this longitudinal 
structure is responsible for generating the mode's grav- 
itational waves, and all of the rest of its fields.) Some- 
what smaller are (i) the longitudinal-transverse compo- 
nents of B (B~q and S^), which together make up the 
longitudinal-transverse part of the frame-drag field, a 2- 
vector B LT parallel to the horizon, and give the horizon- 
piercing vortex lines small non-normal components; and 
(ii) transverse-traceless components B§§ = ~B^ which 
make up the transverse-traceless part of the frame-drag 
field, a 2-tensor B TT parallel to the horizon, and are 
ingoing gravitational waves as seen by Schwarzschild ob- 
servers. (This decomposition into L, LT, and TT parts 
is useful only near the horizon and in the wave zone, 
where there are preferred longitudinal directions associ- 
ated with wave propagation.) 

As the near-zone, longitudinal frame-drag field B L ro- 
tates, it generates a near-zone longitudinal-transverse 
(LT) perturbative frame-drag field B LT via B's prop- 
agation equation (the wave equation for the Riemann 
tensor), and it generates a LT tidal field SS via 
the Maxwell-like Bianchi identity which says, in a local 
Lorentz frame (for simplicity), dS/dt = (V x B) s , where 
the superscript S means "symmetrize" [Eq. (2.15) of Pa- 
per I]. These three fields, B L , B LT , and SS LT together 
maintain each other during the rotation via this Maxwell- 
like Bianchi identity and its (local-Lorentz-frame) dual 
dB/dt = — (V x £) s . They also generate the transverse- 



traceless parts of both fields, B , and 5£ , which be- 
come the outgoing gravitational waves in the wave zone 
and ingoing gravitational waves at the horizon. 

In the equatorial plane, this out going- wave genera- 
tion process, described in terms of vortex and tendex 
structures, is quite pretty, and is analogous to the (1,1) 
magnetic-field mode of Fig. [I] and Sec. I C 4 



As 



moves outward into the induction zone and then the wave 
zone, the equatorial vortex lines bend backward into out- 
going spirals (Fig. |2| and gradually acquire accompany- 
ing tendex lines. The result, locally, in the wave zone, is 
the standard pattern of transverse, orthogonal red and 
blue vortex lines; and (turned by 45 degrees to them) 
transverse, orthogonal red and blue tendex lines, that 
together represent plane gravitational waves (Fig. 7 of 
Paper I). 

It is instructive to focus attention on regions of space 
with large magnitude of vorticity. We call these regions 
vortexes. Figure [3] shows that the equatorial frame-drag 
field consists of four outspiraling vortexes, two red ([light 
gray] counterclockwise) and two blue ([dark gray] clock- 
wise). 




-0.5 



FIG. 3: (color online). The vortex lines (solid black for clock- 
wise; dashed black for counterclockwise) and color-coded vor- 
ticities in the equatorial plane for the same magnetic-parity 
(2, 2) mode as in Fig.[2] This figure differs from Fig. [2] in ways 
designed to give information about the emitted gravitational 
waves: (i) It extends rather far out into the wave zone, (ii) 
It shows the angular structure of the vorticity for the domi- 
nant vortex lines in each region of the equatorial plane. More 
specifically: the color at each point represents the vorticity of 
the equatorial vortex line there which has the largest magni- 
tude of vorticity, with radial variations of vorticity normalized 
away (so the linear color code on the left indicates vorticity 
relative to the maximum at any given radius). The regions 
of large positive vorticity (blue [dark gray]) are clockwise vor- 
texes; those of large negative vorticity (red [light gray]) are 
counterclockwise vortexes. 



The solid black lines in the figure are clockwise vor- 
tex lines. In the clockwise vortexes of the wave zone, 
they have the large magnitude of vorticity that is de- 
picted as blue (dark gray), and they are nearly trans- 
verse to the radial wave-propagation direction; so they 
represent crests of outgoing waves. In the counterclock- 
wise vortexes (red [light gray] regions), these clockwise 
vortex lines have very small magnitude of vorticity and 
are traveling roughly radially, leaping through a red vor- 
tex (a wave trough) from one blue vortex (wave crest) to 
the next. These clockwise vortex lines accumulate at the 
outer edges of the clockwise (blue) vortexes. 

The dashed black lines are counterclockwise vortex 
lines, which are related to the red (light gray), counter- 
clockwise vortexes in the same way as the solid clockwise 
vortex lines are related to the blue (dark gray), clockwise 
vortexes. 

Outside the equatorial plane, this mode also represents 
outgoing gravitational waves, once one gets into the wave 
zone. We depict the strengths of those waves in Fig. [4] 
The blue (dark gray) regions are locations where one vor- 
tex line has vorticity at least 85% of the maximum at that 
radius; in this sense, they are clockwise vortexes. In the 
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Magnetic 



nying driving equation, Eqs. (E9) of Appendix [E 




D^ 1} + (S* + 2tt - 2a)tf i = 



(1.4a) 

(D - 2e)*i = (<S* + 7T - 4a)tf • (1.4b) 




25% 



FIG. 4: (color online). Some 3-dimensional clockwise vortex 
lines (shown black) and regions of large vorticity (vortexes, 
shown blue [dark gray] and red [light gray]) and small vorticity 
(shown off-white), for the same magnetic-parity mode of a 
Schwarzschild black hole as in Figs. [2] and [3] More specifically: 
the inner sphere is the horizon, color coded by its vorticity. 
The blue region is a clockwise vortex in which one vortex 
line has vorticity at least 85% of the maximum value at that 
radius, and similarly for the counterclockwise red region. The 
four off-white regions are locations where no vortex line has 
magnitude of vorticity in excess of 25% of the maximum at 
that radius. 



near zone, two (blue) clockwise vortexes emerge radially 
from the horizon parallel to the plane of the picture, and 
two (red) counterclockwise vortexes emerge radially to- 
ward and away from us. These are 3-dimensional versions 
of the four vortexes emerging from the horizon in the 
equatorial plane of Fig. [3] In the wave zone, the "85%" 
vortexes are concentrated in the polar regions, because 
this mode emits its gravitational waves predominantly 
along the poles. The waves are somewhat weaker in the 
equatorial plane, so although there are spiraling vortexes 
in and near that plane (Fig. [3]), they do not show up at 
the 85% level of Fig. [§ The off-white, spiral -arm struc- 
tures in the equatorial plane represent the four regions 
where the wave strength is passing through a minimum. 

Turn attention from the wave zone to the horizon. 
There the ingoing waves, embodied in B TT and 5S TT 
(which were generated in the near and transition zones 
by rotation of B L ), act back on B L , causing its vortex 
lines to gradually slide off the horizon and thereby pro- 
ducing the mode's exponential decay. 

Just as this process in the electromagnetic ca se is as- 
sociated with the differential conservation law (1.3) for 
magnetic field lines threading the horizon, dB^ / dt 



Here the notation is that of Newman and Penrose: D 
is a time derivative on the horizon, is the mode's 
d£ L + iB L (equivalently S£nn -\-H3nn m disguise), with 
SS L and SSnn vanishing for our mode; \£i is the LT field 



LT 



SB LT (as measured by Schwarzschild observers) 



in disguise; \I/o is the ingoing-wave 5£ TT + iSB TT (as 
measured by Schwarzschild observers) in disguise; 5* is 
a divergence in disguise; and e, ir and a are NP spin 
coefficients. Equation (|1.4b ) says that the ingoing waves 
embodied in 5£ TT +i& dri ve the evolution of the quan- 
tity ^i, and Equation (1.4a) is an approximate differen- 
tial conservation law in which this plays the role of the 
flux of longitudinal vortex lines (number crossing a unit 
length per unit time) and ^2 (i-e., Bnn) is the density 
of longitudinal vortex lines. This differential conserva- 
tion law says that the time derivative of the vortex-line 
density plus the divergence of the vortex-line flux is equal 
to some spin-coefficient terms that, we believe, are gener- 
ally small. (By integrating this approximate conservation 
law over the horizon %, we see that f y Bn ncIA must be 
nearly conserved, in accord with Eq. ( |1.2| ) above, which 
tells us that the horizon integral is nearly zero. In both 
cases, the integral conservation law ( 1.2 ) and the differen- 



( 2 ) V • (aB T ), so also here it is associated with an anal- 
ogous (approximate) conservation law and an accompa- 



tial conservation law (1.4a), it is numerically small spin 
coefficients that slightly spoil the conservation. In Eq. 
(4.3), for a magnetic-parity mode of Schwarzschild and 
Eddington-Finkelstein slicing, we make this conservation 
law completely concrete and find that in this case it is 
precise; there are no small spin coefficients to spoil it. 

Returning to the evolution of the (2,2) magnetic-parity 
mode: The ingoing waves, via Eqs. ( |1.4| ), push the longi- 
tudinal vortex lines away from the centers of the horizon 
vortexes toward their edges (toward the white horizon 
regions in Figs. [2] and [3] At the edges, clockwise vor- 
tex lines from the blue (dark gray) horizon vortex and 
counterclockwise from the red (light gray) horizon vor- 
tex meet and annihilate each other, leading to decay of 
the longitudinal part of the field and thence the entire 
mode. 

We expect to explore this evolutionary process in 
greater detail and with greater precision in future work. 

Turn, next, to spinning black holes. In this case, the 
(2, 2) magnetic-parity mode has qualitatively the same 
character as for a non-spinning black hole. The princi- 
pal change is due to the spin raising the mode's eigen- 
frequency, and the near zone thereby essentially disap- 
pearing, so the perturbed vortex lines that emerge from 
the horizon have a significant back-spiral-induced tilt to 
them already at the horizon. See Fig. [l2| below. 



9 




FIG. 5: (color online). Three-dimensional vortexes for the 
magnetic-parity, (2,1) mode of a Schwarzschild black hole. 
The colored surfaces enclose the region where, for each radius, 
the vorticity of at least one vortex line exceeds 90% of the 
maximum for that radius. In the blue and red (dark and 
light gray) regions, the clockwise and counterclockwise vortex 
lines, respectively, have the larger vorticity. 



7. The (2,1) magnetic-parity mode of a Schwarzschild hole 



tidal field.) 

Close scrutiny of the near-horizon region of Fig. [5] re- 
veals a surprising feature: Within the 90% vortexes (col- 
ored surfaces), the sign of the largest vorticity switches as 
one moves from the near zone into the transition zone — 
which occurs not very far from the horizon; see the in- 
ner dashed circle in Fig. [5] above). This appears to be 
due to the following: The near-zone vortexes are dom- 
inated by the longitudinal part of the frame-drag field 
5B L , which generates all the other fields including SB LT 
via its rotation [see discussion of the (2,2) mode above]. 
The longitudinal-transverse field SB LT is strong through- 
out the near zone and comes to dominate over SB L as 
one moves into the transition zone. Its largest vortic- 
ity has opposite sign from that of 5B L , causing the flip 
of the dominant vorticity and thence the color switch as 
one moves into the transition zone. (Note that a similar 
switch in the sign of the strongest vorticity occurs for the 
magnetic-parity (2,2) mode vortexes illustrated in Fig. 4j 
although there the transition occurs farther out, at the 
edge of the wave zone.) 

In Sees. | V A| and | V B] we explore in considerable detail 
this magnetic-parity (2,1) mode and also its near dual, 
the electric-parity (2,1) mode, focusing especially on the 
shapes of their vortexes. 



For the (2,1) magnetic-parity mode of a Schwarzschild 
black hole, there are two horizon vortexes in the hole's 
northern hemisphere (one counterlockwise, the other 
clockwise), and two in the southern hemisphere. From 
these emerge the longitudinal part of the frame-drag field, 
in the form of four 3D vortexes (Fig. [5|. 

These four vortexes actually form two spiral arms, each 
of which contains vortex lines of both signs (clockwise 
and counterclockwise). The surface of each arm is color 
coded by the sign of the vorticity that is largest in mag- 
nitude in that region of the arm. This dominant vor- 
ticity flips sign when one passes through the equatorial 
plane — from positive (i.e., blue [dark gray]; clockwise) 
on one side of the equator to negative (i.e., red [light 
gray]; clockwise) on the other side. The reason for this 
switch is that for m = 1 the e irn ^ angular dependence 
means reflection antisymmetry through the polar axis, 
which combined with the positive parity of the I = 2 
frame-drag field implies reflection antisymmetry through 
the equatorial plane. The (2,2) mode of the previous sec- 
tion, by contrast, was reflection symmetric through both 
the polar axis and the equatorial plane. 

By contrast with the (2,2) mode, whose region of 
largest vorticity switched from equatorial in the near zone 
to polar in the wave zone (Fig. [3]), for this (2,1) mode, the 
region of largest vorticity remains equatorial in the wave 
zone. In other words, this mode's gravitational waves are 
stronger in near-equator directions than in near-polar di- 
rections. (Recall that in the wave zone, the vortexes are 
accompanied by tendexes with tendicities equal in mag- 
nitude to the vorticities at each event, so we can discuss 
the gravitational-wave strengths without examining the 



8. The (2,0) magnetic-parity mode of a Schwarzschild hole 

The (2,0) magnetic-parity mode has very different dy- 
namical behavior from that of the (2,1) and (2,2) modes. 
Because of its axisymmetry, this mode cannot be gen- 
erated by longitudinal, near-zone vortexes that rotate 
around the polar axis, and its waves cannot consist of 
outspiraling, intertwined vortex and tendex lines. 

Instead, this mode is generated by longitudinal, near- 
zone vortexes that oscillate, and its waves are made up 
of intertwined vortex lines and tendex lines that wrap 
around deformed tori. These gravitational-wave tori re- 
semble smoke rings and travel outward at the speed of 
light. More specifically: 

Because of axisymmetry, the (2,0) magnetic-parity 
mode has one family of vortex lines that are azimuthal 
circles of constant r and 0, and two families that lie in 
surfaces of constant <p. Figure [6] is a plot in one of 
these surfaces. (The plot for any other <ft will be iden- 
tical to this, by axisymmetry.) This plot shows the vortex 
lines that lie in <S^, and by color coding at each point, 
the vorticity of the strongest of those lines. 

Notice that, at this phase of oscillation, there are clock- 
wise (solid) vortex lines sticking nearly radially out of the 
horizon's polar regions and counterclockwise (dashed) 
vortex lines sticking nearly radially out of the horizon's 
equatorial region. A half cycle later the poles will be 
red (light gray) and equator blue (dark gray). These 
near-zone vortex lines are predominantly the longitudinal 
part of the frame-drag field JB L , which we can regard as 
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FIG. 6: (color online). Vortex lines and vorticities for 
magnetic-parity (2,0) mode of Schwarzschild in a surface S</> of 
constant <j>. The line and coloring conventions are the same 
as in Fig. [3] (solid lines for clockwise, dashed for counter- 
clockwise; color shows vorticity of the vortex line with largest 
magnitude of vorticity, with radial variations removed and in- 
tensity of color as in the key on right edge of Fig. |3j). The 
central circle is the horizon, color coded by the horizon vor- 
ticity. 



working hand-in-hand with the near-zone, longitudinal- 
transverse tidal field 5£ LT to generate the other fields. 

As we shall see in Sec. V C| (and in more convincing 
detail for a different oscillatory mode in Sec. IV C), the 
dynamics of the oscillations are these: Near-zone energy 1 



1 We use the term energy in a generalized and descriptive sense 
here and elsewhere in this paper. We note, however, that with 
a suitable (nonunique) definition of local energy, we can make 
these notions more precise. For example, the totally symmet- 
ric, traceless Bel-Robinson tensor serves as one possible basis for 
this. In vacuum it is T^ vpa = 1/2(C ^ av pC ' P a J + C* au pC; a J) 
with * denoting the Hodge dual, and it is completely symmet- 
ric and obeys the differential conservation law V p T^ upc r = 0. 
Given a unit timelike slicing vector u we conveniently have 
W(u) = T^paU^u^uPu" 7 = l/2(E ij E i i + B^B^) > as a 
positive-definite superenergy built from the squares of the tidal 
and frame-drag fields in a given slice (see the reprint of Bel's ex- 
cellent paper |10| for motivation and definition, e.g., Penrose and 
Rindler 5^ for the spinor representation of the Bel- Robinson ten- 
sor, and e.g., [J for its relation to notions of quasilocal energy). 
As another example, magnetic-parity modes of Schwarzschild are 
describable by the Regge- Wheeler function Q(r*,t) which satis- 
fies t he Sturm-Liouville equation Q, r *r* — Q,tt — V(r*)Q = [Eq. 
( |A5| ) but with the e~ tujt time dependence absorbed into Q]. The 
integral conservation law associated with this Sturm-Liouville 

equation is d/dt J* (Q 2 ^ + Q 2 t + VQ 2 ) dr* = 2Q, r *Q,t |*. The 
quantity inside the integral can be regarded as an energy density, 



oscillates back and forth between the near-zone 5B L , and 
the near-zone SB LT and 5£ LT . As SB L decays, its vor- 
tex lines slide off the hole and (we presume) form closed 
loops, lying in S^, which encircle outgoing deformed tori 
of perturbed tendex lines that become the transverse- 
traceless gravitational waves. Only part of the energy in 
SB L goes into the outgoing waves. Some goes into the TT 
ingoing waves, and the rest (a substantial fraction of the 
total energy) goes into SB LT and 5£ LT , which then use 
it to regenerate 5B L , with its horizon-penetrating vortex 
lines switched in sign (color), leading to the next half 
cycle of oscillation. 

The vortex lines that encircle the gravitational-wave 
tori are clearly visible in Fig. [6] Each solid (clockwise) 
line is tangential (it points nearly in the direction) when 
it is near the crest (the maximum-vorticity surface) of a 
blue (dark gray), lens-shaped gravitational- wave vortex. 
As it nears the north or south pole, it swings radially 
outward becoming very weak (low vorticity) and travels 
across the red trough of the wave, until it nears the next 
blue crest. There it swings into the transverse, direction 
and travels toward the other pole, near which it swings 
back through the red trough and joins onto itself in the 
original blue crest. 

Each dashed (counterclockwise) closed vortex line be- 
haves in this same manner, but with its transverse por- 
tions lying near red (light gray) troughs (surfaces of most 
negative vorticity). Near the red troughs, there are blue 
azimuthal vortex lines (not shown) that encircle the hole 
in the <fi direction, and near the blue crests, there are red 
azimuthal lines. 

Figure [7] sheds further light on these gravitational-wave 
tori. It shows in three dimensions some of the pertur- 
bative tendex lines for the (2,0) magnetic-parity mode 
that we are discussing. (For this mode, two families of 
perturbative tendex lines, one red [counterclockwise] and 
the other blue [clockwise] , have nonzero tendicity and the 
third family has vanishing tendicity.) As is required by 
the structure of a gravitational wave (transverse tendex 
lines rotated by 45 degrees relative to transverse vortex 
lines), these perturbative tendex lines wind around tori 
with pitch angles of 45 degrees; one family winds clock- 
wise and the other counterclockwise, and at each point 
the two lines have the same magnitude of vorticity. 

A close examination of Fig. [7] reveals that the tori 
around which the perturbative tendex lines wrap are half 
as thick as the tori around which the vortex lines wrap. 
Each tendex-line torus in Fig. [7] is centered on a single 



and the quantity on the right hand side an energy flux. For the 
(2,ra) magnetic-parity mode, Eqs. ( |A12a| and ( |A12d[ ) express Q 
in terms of the time derivative of the longitudinal part of 53 with 
its angular dep enden ce Y 2 m re moved: Q = (r 3 /12) 35 Bf.f./dt. 



Others of Eqs. \A12\ and \A26\ relate Q, r * to the LT parts of 
53 and 5S. This could be the foundation for a second way to 
make more precise the notion of energy fed back and forth be- 
tween the various parts of 53 and 5S. 
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FIG. 7: (color online). Positive-tendicity (blue) and negative- 
tendicity (red) perturbative tendex lines of a (2, 0) magnetic- 
parity perturbation of a Schwarzschild black hole. These lines 
spiral around deformed tori of progressively larger diameter. 
The viewpoint is looking down onto the equatorial plane from 
the positive symmetry axis. Upper right inset: The negative 
tendex line spiraling around the outermost torus, viewed in 
cross section from the equatorial plane. Lower right inset: 
The negative tendex line spiraling around the small torus 
third from the center, viewed in cross section from the equa- 
torial plane. Also shown, in black in the main drawing and 
the large inset, are two of this mode's vortex lines, one from 
Fig. [6] wrapping around the outermost torus in a plane; the 
other an azimuthal circle wrapping around that torus in the 
direction. This figure was actually drawn depicting vortex 
lines of the electric-parity mode discussed in Sec. | VP] but by 
duality (which is excellent in the wave zone), it also repre- 
sents the tendexes of the magnetic-parity mode discussed in 
this section. 



node of the gravitational- wave field; the thick red torus 
in the upper right panel reaches roughly from one crest 
of the wave to an adjacent trough. By contrast, each 
vortex-line torus (Fig. [6] and black poloidal curves in Fig. 
[7| reach from crest to crest or trough to trough and thus 
encompass two gravitational- wave nodes. 

Each node in the wave zone has a family of nested 
tendex-line tori centered on it. The four tendex-line tori 
shown in Fig. [7] are taken from four successive families, 
centered on four successive nodes. The second thin torus 
is from near the center of one nested family; it tightly 
hugs a node and therefore has near vanishing tendicity. 
The two thick tori are from the outer reaches of their 
nested families. 



For further details of the (2,0) modes, see Sees. VC 
and [VD] below. 



9. The superposed (2, 2) and (2,-2) magnetic-parity mode 
of a Schwarzschild hole 

As we have seen, the magnetic-parity, (2, 2) mode of 
a Schwarzschild black hole represents vortexes that ro- 
tate counterclockwise around the hole, spiraling outward 
and backward (Figs. [5J [3] and [4] above). If we change 
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FIG. 8: (color online). Equatorial vortex structure of the 
superposed (2,2) and (2,-2), magnetic-parity, fundamental 
modes of a Schwarzschild black hole. The colors encode the 
vorticity of the dashed vortex lines. The vorticity of the solid 
lines is not shown, but can be inferred from the fact that 
under a 90° rotation, the dashed lines map into solid and the 
solid into dashed. 



the sign of the azimuthal quantum number to m = —2, 
the vortexes rotate in the opposite direction, and spi- 
ral in the opposite direction. If we superpose these two 
modes (which, for Schwarzschild, have the same eigenfre- 
quency), then, naturally, we get a non-rotating, oscilla- 
tory mode — whose dynamics are similar to those of the 
(2,0) mode of the last subsection. See Sec. 



IV 



for details. 

Figure [8] is a snapshot of the two families of vortex lines 
that lie in this mode's equatorial plane. The plane is col- 
ored by the vorticity of the dashed vortex lines; they are 
predominantly counterclockwise (red), though in some 
regions they are clockwise (blue). 

The red (light gray) regions form interleaved rings 
around the black hole, that expand outward at the speed 
of light, along with their dashed vortex lines. These rings 
are not tori in three dimensions because [by contrast with 
the (2,0) mode] the frame drag field grows stronger as 
one moves up to the polar regions, rather than weaken- 
ing. As the mode oscillates, the longitudinal near-zone 
frame-drag field B L , which drives the mode, generates 
new interleaved rings, one after another and sends them 
outward. 

During the oscillations, there are phases at which the 
longitudinal field B L threading the hole goes to zero, and 
so the hole has vanishing horizon vorticity. The near-zone 
oscillation energy, at these phases, is locked up in the 
near-zone, longitudinal-transverse fields B LT and SS LT , 
which, via the Maxwell-like Bianchi identities (and the 
propagation equation that they imply), then feed en- 
ergy into the longitudinal near-zone frame-drag field B L , 
thereby generating new horizon-threading vortex lines, 



12 



which will give rise to the next ejected interleaved ring. 

We explore these dynamics in greater detail in 
Sec. ITVCl 



D. This paper's organization 

The remainder of this paper is organized as follows: 
In Sec. |TI| we introduce the time slicing and coordi- 
nates used throughout this paper for the background 
Schwarzschild and Kerr spacetimes, we introduce the 
two gauges that we use for Schwarzschild perturbations 
(Regge-Wheeler-Zerilli and ingoing radiation gauges) and 
the one gauge (ingoing radiation) we use for Kerr, we 
discuss how our various results are affected by changes 
of gauge, and we discuss how we perform our computa- 
tions. In Sees. |III[ [TV] and [V| we present full details of 
our results for the fundamental (most slowly damped) 
quadrupolar mod es of Schwarzschild and Kerr: (2,2) 
modes in Sec. |III| superposed (2,2) and (2,-2) modes in 
Sec.[lVl and both (2,1) and (2,0) modes in Sec.[V| In Sec. 



|III D we compare vortex lines computed in a numerical- 
relativity simulation of a binary black hole at a late time, 
when the merged hole is ringing down, with the vortex 
lines from this paper for the relevant quasinormal mode; 
we obtain good agreement. In Sec. |VI[ we make a few 
concluding remarks. And in six appendices, we present 
mathematical details that underlie a number of this pa- 
per's computations and results. 



II. SLICINGS, GAUGES AND 
COMPUTATIONAL METHODS 

When calculating the tidal and frame-drag fields of 
perturbed black- hole spacetimes, we must choose a slic- 
ing and also spatial coordinates on each slice, for both 
the background spacetime and at first order in the per- 
turbations ("perturbative order"). The perturbative- 
order choices of slicing and spatial coordinates are to- 
gether called the chosen gauge. We will always use the 
same choice of background slicing and coordinates in this 
study, but we will use different choices for our gauge. 

This section describes the choices we make, how they 
influence the vortex and tendex lines and their vorticities 
and tendicities (which together we call the "vortex and 
tendex structures"), and a few details of how, having 
made our choices, we compute the perturbative frame- 
drag and tidal fields and the vortex and tendex struc- 
tures. Most of the mathematical details are left to later 
sections and especially appendices. 

In Sec. |II A[ we describe our choices of slicing and spa- 
tial coordinates. In Sec. |IIB[ we sketch how we calculate 
the perturbative frame-drag and tidal fields and visual- 

IICl we 



ize their vortex and tendex structures. In Sec. 
explore how those structures change under changes of 
gauge, i.e., changes of the perturbative slicing and per- 
turbative spatial coordinates. 



A. Slicing, spatial coordinates, and gauge 

Throughout this paper, for the background (unper- 
turbed) Kerr spacetime, we use slices of constant Kerr- 
Schild (KS) time t, which is related to the more familiar 
Boyer-Lindquist time by 



t = t + r* — r , where 

dr 



(2.1) 



(Eq. (6.2) of Paper II [4]). Here t and r are the 
Boyer-Lindquist time and radial coordinates, a is the 
black hole's spin parameter (angular momentum per unit 
mass), and A = r 2 — 2Mr + a 2 , with M the black- hole 
mass. Our slices of constant t penetrate the horizon 
smoothly, by contrast with slices of constant t, which are 
singular at the horizon. In the Schwarzschild limit a — >• 0, 
t and r become Schwarzschild' s time and radial coordi- 
nates, and t becomes ingoing Eddington-Finkelstein time, 
t = £ + 2Mln|r/2M- 1|. 

On a constant-t slice in the background Kerr space- 
time, we use Cartesian-like KS (Kerr-Schild) spatial co- 
ordinates, when visualizing vortex and tendex structures; 
but in many of our intermediary computations, we use 
Boyer-Lindquist spatial coordinates {r, 0, <f\ (which be- 
come Schwarzschild as a — >• 0). The two sets of coordi- 
nates are related by 



x + iy = (r + ia)e 1 ^ 



suit 



r cost 



[Eq. (6.7) of Paper II]. Here 



s: 



-dr 



(2.2) 



(2.3) 



[Eq. (6.5) of Paper II] is an angular coordinate that, 
unlike (/>, is well behaved at the horizon. In the 
Schwarzschild limit, the KS {x, y, z} coordinates become 
the quasi-Cartesian {x, y, z} associated with Eddington- 
Finkelstein (EF) spherical coordinates {r, #,0}. 

Our figures (e.g., [2jj8] above) are drawn as though 
the KS y, z} were Cartesian coordinates in flat 
spacetime — i.e., in the Schwarzschild limit, as though the 
EF {r, 0, (/)} were spherical polar coordinates in flat space- 
time. 

We denote by gffi the background metric in KS space- 
time coordinates [Eq. (6.8) of Paper II] (or EF spacetime 
coordinates in the Schwarzschild limit). When the black 
hole is perturbed, the metric acquires a perturbation h^ v 
whose actual form depends on one's choice of gauge — i.e., 
one's choice of slicing and spatial coordinates at pertur- 
bative order. 

For Schwarzschild black holes, we use two different 
gauges, as a way to assess the gauge dependence of our 
results: (i) Regge-Wheeler-Zerilli (RWZ) gauge, in which 
is a function of two scalars (Q for magnetic parity and 
Z for electric parity) that obey separable wave equations 
in the Schwarzschild spacetime and that have spin- weight 
zero (see App.[A|for a review of this formalism), and (ii) 
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ingoing radiation (IR) gauge, in which h^ v is computed 
from the Weyl scalar \I>o (or ^4) that obeys the separa- 
ble Bardeen-Press equation. The method used to com- 
pute the metric perturbation from \I/o is often called the 
Chrzanowski-Cohen-Kegeles (CCK) procedure of metric 
reconstruction (see App. |C| ). 

In App. |DJ we exhibit explicitly the relationship be- 
tween the RWZ and IR gauges, for electric- and magnetic- 
parity perturbations. The magnetic-parity perturbations 
have different perturbative spatial coordinates, but the 
same slicing. (In fact, all gauges related by a magnetic- 
parity gauge transform have identically the same slicing 
for magnetic-parity perturbations of Schwarzschild [al- 
though the same is not true for Kerr]; see Sec. II C\ . For 
electric-parity perturbations, the two gauges have differ- 
ent slicings and spatial coordinates. 

For all the perturbations that we visualize in this pa- 
per, the tendexes and vortexes show quite weak gauge 
dependence. See, e.g., Sec. |III[ where we present results 



from both gauges. The results in Sees. [TV] and [VJ are all 
computed in RWZ gauge. 

For Kerr black holes, there is no gauge analogous to 
RWZ; but the IR gauge and the CCK procedure that 
underlies it are readily extended from Schwarzschild to 
Kerr. In this extension, one constructs the metric per- 
turbation from solutions to the Teukolsky equation (see 
App. [B| for the perturbations to the Weyl scalars ^0 an d 
^4, in an identical way to that for a Schwarzschild black 
hole described above. Our results in this paper for Kerr 
black holes, therefore, come solely from the IR gauge. 



B. Sketch of computational methods 

This section describes a few important aspects of how 
we calculate the tidal and frame-drag fields, and their 
vortex and tendex structures which are visualized and 
discussed in Sees. |III[ \LV\ and|V| 

We find it convenient to solve the eigenvalue problem 
in an orthonormal basis (orthonormal tetrad) given by 
the four-velocities of the Kerr-Schild (KS) or Eddington- 
Finkelstein (EF) observers, and a spatial triad, e?a, car- 
ried by these observers. 

The background EF tetrad for the Schwarzschild 
spacetime, expressed in terms of Schwarzschild coordi- 
nates, is 



e1°) 



1 



y/l + 2M/r 



y/l + 2M/r 



ld_ 

r 06 ' 



e{ 0) = 



Id 2M d 

a 2 dt r dr 

d 2M d 
dr a 2 r dt 
1 d 



r sin( 



(2.4) 



[cf. Eqs. (4.4) of Paper II, which, however, are writ- 
ten in terms of the EF coordinate basis rather than 
Schwarzschild]. The background orthonormal tetrad for 



KS observers (in ingoing Kerr coordinates {£, r, 0, (/>}; see 
Paper II, Sec. VI C) is 



U^=Hd;- 2 -^d r . 



(0) _ VA . aH 



el 



el ' 



£ 1 
A sin 9 



9^ 



d } , (2.5a) 



where we have defined 
£ = r 2 + a 2 cos 2 , 



A= (r + a z Y - a\r z - 2Mr + a 2 ) sin 2 



(2.5b) 
(2.5c) 
(2.5d) 



[see Eq. (B2) of Paper II]. 

When the black hole is perturbed, the tetrad 
{u^ , , e^ , e } acquires perturbative corrections 
that keep it orthonormal with respect to the metric 

9nu = 9$ +^/izy- We choose the perturbative corrections 
to the observers' 4-velocity so as to keep it orthogonal 
to the space slices, i.e., so as to keep u = — aVt. (Here 
a = dr /dt, differentiating along the observer's world line, 
is the observer's lapse function.) A straightforward cal- 
culation using the perturbed metric gives the following 
contravariant components of this u: 



1 — -r ■ 



\1) 



2 '*oo '""(o) v - 



(2.6) 



where = h^u^u^y and is the four- velocity of 
the background observers. 

We choose the perturbative corrections to the spatial 
triad {ej} so the radial vector stays orthogonal to sur- 
faces of constant r in slices of constant t, the direction 
continues to run orthogonal to curves of constant in 
surfaces of constant r and t, and the <\> vector changes 
only in its normalization. 

When written in terms of the unperturbed tetrad and 
projections of the metric perturbation into the unper- 
turbed tetrad, the perturbation to the tetrad then takes 
the form 



= 




(2.7a) 


4 11 = 


^'hrdf "■f J 4 e (0) ' 


(2.7b) 


=1" = 


2 09 § "6>0 e (O) ' 


(2.7c) 




--h--ei 0) 
2 H ' 


(2.7d) 



where i is summed over f, 6, and 0, and A is summed 
over only and <j). 

In Appendices |A| (RWZ gauge) and [C] (IR gauge), we 
give the details of how we compute the components 



13 ij ij 



13 IJ IJ 



(2.8) 
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of the tidal and frame-drag field in this perturbed tetrad. 
The background portions £~? and are the station- 
ary fields of the unperturbed black hole, which were 
computed and visualized in Paper II. The perturbative 
pieces, and B^V are the time-dependent, perturba- 
tive parts, which carry the information about the quasi- 
normal modes, their geometrodynamics, and their grav- 
itational radiation. 

As part of computing the perturbative £~V = S£^, 

and = 5B^ for a chosen quasinormal mode of a 

ZJ 3 

Kerr black hole, we have to solve for the mode's Weyl- 
scalar eigenfunctions and and eigenfrequency 
lj. To compute the frequencies, we have used, through- 
out this paper, Emanuele Berti's elegant computer code 
[T5], which is discussed in [13^ and is an implementation 
of Leaver's method [13] . To compute the eigenfunctions, 
we our own independent code (which also uses the same 
procedure as that of Berti). In App. [Cj we describe 
how we extract the definite-parity (electric or magnetic) 
eigenfunctions from the non-definite-parity functions. 

To best visualize each mode's geometrodynamics and 
generation of gravitational waves in Sees. |III[ |IV[ andjV] 
we usually plot the tendex and vortex structures of the 

perturbative fields £~V and B~V . However, when we com- 

ij ij 

pare our results with numerical-relativity simulations, it 
is necessary to compute the tendex and vortex structures 
of the full tidal and frame-drag fields (background plus 
perturbation), because of the difficulty of unambiguously 
removing a stationary background field from the numer- 
ical simulations. As one can see in Figs. [15] and [26] in 
this case much of the detail of the geometrodynamics and 
wave generation is hidden behind the large background 
field. 

In either case, the tendex and vortex structure of the 
perturbative fields or the full fields, we compute the field 
lines and their eigenvalues in the obvious way: At se- 
lected points on a slice, we numerically solve the eigen- 
value problem 

E~V S = \\] (2.9) 



C. Gauge changes: Their influence on tidal and 
frame-drag fields and field lines 

For perturbations of black holes, a perturbative gauge 
change is a change of the spacetime coordinates, x a = 
x a + that induces changes of the metric that are of 
order the metric perturbation; when dealing with defi- 
nite parity perturbation, we split the generator of the 
transform ^ a into definite electric- and magnetic-parity 
components. The gauge change has two parts: A change 
of slicing generated by £°, and a change of spatial coor- 
dinates 

F=t + f° x j ' =a^ + C j . (2.11) 

Here all quantities are to be evaluated at the same event, 
V, in spacetime. 

Because £° is a scalar under rotations in the 
Schwarzschild spacetime — and all scalar fields in 
Schwarzschild have electric parity — for a magnetic-parity 
£° vanishes, and the slicings for magnetic-parity 
quasinormal modes of Schwarzschild are unique. For 
these modes, all gauges share the same slicing (see Ap- 
pendix JdJ. 2 

1. Influence of a perturbative slicing change 

For (electric-parity) changes of slicing, the new ob- 
servers, whose world lines are orthogonal to the new 
slices, t' =const, move at velocity 

Av = -aVf (2.12) 

with respect the old observers, whose world lines are or- 
thogonal to the old slices t =const). Here V is the gra- 
dient in the slice of constant ?, and a = (dr/dt) is the 
lapse function, evaluated along the observer's worldline. 
In other words, Av is the velocity of the boost that leads 
from an old observer's local reference frame to a new 
observer's local reference frame. Just as in electromag- 
netic theory, this boost produces a change in the ob- 
served electric and magnetic fields for small Av given by 
AB = Av x E and AE = —Av x B, so also it produces a 



for the three eigenvalues A and unit-normed eigenvectors 
V- of and similarly for By, and we then compute 
the integral curve (tendex or vortex line) of each eigen- 
vector field by evaluating its coordinate components V 3 ' 
in the desired coordinate system (KS or EF) and then 
numerically integrating the equation 



dx 3 
ds 



V 3 



(2.10) 



where s is the proper distance along the integral curve. 



2 In the Kerr spacetime, however, there are magnetic-parity 
changes of slicing, because £° no longer behaves as a scalar un- 
der rotations. To understand this more clearly, consider, as a 
concrete example, a vector in Boyer-Lindquist coordinates with 
covariant components £^ = (0, 0, X l e m , X l ( ^ n )f(r)e~' lU}t , where 
X 1 ™ are the components of a magnetic-parity vector spherical 
harmonic [see Eq. ( |C26a| )]. This vector's contravariant compo- 
nents are ^ = (g^X^, 0, g 00 X l e m , X l ™)f(r)e~ %UJ \ where 
gt<j> ^ g86 ^ anc j g(p(f) contr aval i&rit components of the Kerr 

metric (which have positive parity). The vector has magnetic 
parity and a nonvanishing component £°; therefore, it is an ex- 
ample of a magnetic-parity gauge-change generator in the Kerr 
spacetime that changes the slicing. 
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change in the observed tidal and frame-drag fields given 

by 



AB = (Av x £) s , A£ = -(Avx&f 



(2.13) 



(e.g., Eqs. (A12) and (A13) of [15], expanded to linear or- 
der in the boost velocity). Here the superscript S means 
symmetrize. 



2. Example: Perturbative slicing change for Schwarzs child 
black hole 



For a Schwarzschild black hole, because the unper- 
turbed frame-drag field vanishes, AS is second order in 
the perturbation and thus negligible, so the tidal field is 
invariant under a slicing change. By contrast, the (fully 
perturbative) frame-drag field can be altered by a slicing 
change; AB = (Av x S) s is nonzero at first order. 

Since the unperturbed tidal field is isotropic in the 
transverse (0, (j>) plane, the radial part of Av produces 
a vanishing AB. The transverse part of Av, by contrast, 
produces a radial-transverse AB (at first-order in the per- 
turbation). In other words, a perturbative slicing change 
in Schwarzschild gives rise to a vanishing AS and an 
electric-parity AB whose only nonzero components are 



AB^ 



ABfr and AB, 



AS, 



(2.14) 



For a Schwarzschild black hole that is physically unper- 
turbed, the first-order frame-drag field is just this radial- 
transverse AB, and its gauge-generated vortex lines make 
45 degree angles to the radial direction. 



(2.15)] on the line's eigenvector, and then (ii) compute 
the integral curve of the altered eigenvector field. 

Far simpler than this brute-force approach is to note 
that the tendex line, written as location V(s) in the 
slice of constant f as a function of spatial distance s 
along the curve, is unaffected by the coordinate change. 
Therefore, if the old coordinate description of the ten- 
dex line is (s) — xi[P(s)], then the new coordinate 
description is x j (s) = x^[P(s)] +&\P(s)\\ i.e., x J (s) = 
x^(s) + &[xi(s)]. In other words, as seen in the new 
(primed) coordinate system, the tendex line appears to 
have been moved from its old coordinate location, along 
the vector field £ J , from its tail to its tip; and similarly 
for any vortex line. 



4- Example: Perturbative spatial coordinate change for a 
Schwarzschild black hole 



Because the frame- drag field of a perturbed 
Schwarzschild black hole is entirely perturbative, it 
is unaffected by a spatial coordinate change. This, 
together with AB = for magnetic-parity modes implies 
that the frame- drag field of any magnetic-parity mode of 
Schwarzschild is fully gauge invariant! 

By contrast, a spatial coordinate change (of any parity) 
mixes some of the background tidal field into the pertur- 
bation, altering the coordinate locations of the tendex 
lines. 

As an example, consider an electric-parity (2,2) mode 
of a Schwarzschild black hole. In RWZ gauge and in the 
wave zone, the tidal field is given by 



3. Influence of perturbative change of spatial coordinates 

Because S and B are tensors that live in a slice 
of constant t, the perturbative change of spatial co- 
ordinates, which is confined to that slice, produces 
changes in components that are given by the stan- 
dard tensorial transformation law, Si>j'(x k [V]) = 
£ pq (x k [P])(dx q /dx % )(dx p /dxi ). To first order in the 
gauge-change generators x k , this gives rise to the follow- 
ing perturbative change in the tidal field 

~ —£ij\k£, k — £ik£, k \j — £jk£, k \i j (2.15) 

and similarly for the frame-drag field B. Here the sub- 
script I denotes covariant derivative with respect to the 
background metric, in the slice of constant t. The two ex- 
pressions in Eq. (2.15) are equal because the connection 



coefficients all cancel. 

The brute-force way to compute the influence of a spa- 



tial coordinate change x J 



£ J on the coordinate 



shape x^ (s) of a tendex line (or vortex line) is to (i) solve 
the eigenequation to compute the influence of ASij [Eq. 
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A 




f - 
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A 
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\a 






cos 


ujv 2 





2M 



O 



3 [2(<£ - (f> ) - u(t - r,)] 



A 



(2.16) 



where A is the wave amplitude 

Focus on radii large enough to be in the wave zone, 
but small enough that the wave's tidal field is a small 
perturbation of the Schwarzschild tidal field. Then the 
equation for the shape of the nearly circular tendex lines 
that lie in the equatorial plane, at first order in the wave's 
amplitude, is 



1 dr 



2rA 



ii' ic c \ onr cos[2((/> - O ) - u(t - r*)] 
r dqb (S^ - Srr) 3Muj 

(2.17) 

(an equation that can be derived using the standard per- 
turbation theory of eigenvector equations). Solving for 
r(4>) using perturbation theory, we obtain for the tendex 
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line's coordinate location 

r(<M) = r o + p(0,r o ,t) , 
Ar 



p{<t>,r ,t) 



3Muj 



(2.18) 

sin[2(</> - O ) - - r G *)] . 



Here r is the radius that the chosen field line has when 
(j) = (j) . Notice that the field line undergoes a quadrupo- 
lar oscillation, in and out, as it circles around the black 
hole, and it is closed — i.e., it is an ellipse centered on the 
hole. The ellipticity is caused by the gravitational wave. 
As time passes, the ellipse rotates with angular velocity 
d<j)/dt = uo j '2, and the phasing of successive ellipses at 
larger and larger radii r Q is delayed by an amount corre- 
sponding to speed-of-light radial propagation. 

Now, consider an unperturbed Schwarzschild black 
hole. We can produce this same pattern of elliptical 
oscillations of the equatorial-plane tendex lines, in the 
absence of any gravitational waves, by simply changing 
our radial coordinate: Introduce the new coordinate 



-£ r , where £ r = p(</>, r, t) , 



(2.19) 



with p the function defined in Eq. ( |2.18 ). In 
Schwarzschild coordinates, the equatorial tendex lines are 
the circles r = r Q =constant. In the new coordinate sys- 
tem, those tendex lines will have precisely the same shape 
as that induced by our gravitational wave [Eq. ( 2.18| )]: 
r' = r Q + p((/>, r , t). Of course, a careful measurement of 
the radius of curvature of one of these tendex lines will 
show it to be constant as one follows it around the black 
hole (rather than oscillating), whereas the radius of cur- 
vature of the wave- influenced tendex line will oscillate. In 
fact, if we follow along with the tendex line and measure 
the tendicity along the line, we find that the tendicity of 
the line is unchanged by the change in coordinates. To 
be explicit, consider the tendicity, which we denote A^, 
along one of the lines r = r Q . Enacting the coordinate 
transform on the tendicity but continuing to evaluate it 
along the perturbed line, we have the identity 



A^(r ) . 



A^(r ■ 



P-p) 
(2.20) 



Nevertheless, if one just casually looks at the 
Schwarzschild tendex lines in the new, primed, coordi- 
nate system, one will see a gravitational-wave pattern. 

The situation is a bit more subtle for the perturbed 
blac k hole . In this case, the tendex lines are given by 
Eq. (2.18), and we can change their ellipticity by again 
changing radial coordinates, say to 



,r,t) . 



(2.21) 



The radial oscillations Ar' of the elliptical tendex lines in 
the new (r',0, t) coordinate system will have amplitudes 
1+a times larger than in the original (r, (/>, t) coordinates, 
and in the presence of the gravitational waves it may not 
be easy to figure out how much of this amplitude is due 



to the physical gravitational waves and how much due to 
rippling of the coordinates. 

On the other hand, the tendicities of these tendex lines 
are unaffected by rippling of the coordinates. They re- 
main equal to = = M/r 3 + (A/r) cos[2(0 — (j) ) — 
u(t-r*)] =M/(r / ) 3 + (A/r / )cos[2(0-0 o )-cj(t-<)] at 
leading order, which oscillates along each closed line by 
the amount A£^ = (A/r) cos[2(0 — <p ) — uo(t — r*)] that 
is precisely equal to the gravitational- wave contribution 
to the tendicity. Note that in this example, even with- 
out evaluating the tendicity along the perturbed lines to 
cancel the coordinate change, the change in the tendic- 
ity due to the coordinate change enters at a higher order 
than the contribution from the gravitational wave. 

Therefore, in this example, the tendicity and corre- 
spondingly the structures of tendexes capture the grav- 
itational waves cleanly, whereas the tendex-line shapes 
do not do so; the lines get modified by spatial coordi- 
nate changes. This is why we pay significant attention 
to tendexes and also vortexes in this paper, rather than 
focusing solely or primarily on tendex and vortex lines. 



III. (2,2) QUASINORMAL MODES OF 
SCHWARZSCHILD AND KERR BLACK HOLES 



In Sec. |I C 6[ we described the most important fea- 
tures of the fundamental, (2,2) quasinormal modes of 
Schwarzschild black holes. In this section, we shall ex- 
plore these modes in much greater detail and shall extend 
our results to the (2,2) modes of rapidly spinning Kerr 
black holes. For binary-black- hole mergers, these are the 
dominant modes in the late stages of the merged hole's 
final ringdown (see, e.g., [16]). 



A. Horizon vorticity and tendicity 

We can compute the horizon tendicity £nn and vortic- 
ity Sat at [or equivalent ly ^2 = ^(£nn~\-H3nn)] using two 
methods: first, we can directly evaluate them from the 
metric p ertur bations, and s econd , we can calculate them, 
via Eq. ( |E12[ ) in the form (E14), from the ingoing- wave 
curvature perturbation ^0, which obeys the Teukolsky 
equation (App. [B] ). For perturbations of Schwarzschild 
black holes, both methods produce simple analytical ex- 
pressions for the horizon quantities; they both show 
that the quantities are proportional to a time-dependent 
phase times a scalar spherical harmonic, e~ %UJt Yi m [see, 



e.g., Eq. (E19)]. For Kerr holes, the sim plest formal ex- 
pression for the horizon quantities is Eq. ( |E14 ), and there 
is no very simple analytical formula. Nevertheless, from 
these calculations one can show that there is an exact du- 
ality between £nn and Sat at in ingoing radiation gauge 
for quasinormal modes with the same order parameters 
(n, /,ra) but opposite parity, for both Schwarzschild and 
Kerr black holes; see App. |E 2| For Schwarzschild black 
holes in RWZ gauge, there is also a duality for the horizon 
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quantities, although it is complicated by a perturbation 
to the position of the horizon in this gauge; see Appen- 
dices. [A4] and [A5] for further discussion. 

In Fig. [9j we show SEnn and SBnn for the (2,2) 
modes with both parities, of a Schwarzschild black hole 
(upper row) and a rapidly rotating Kerr black hole (bot- 
tom row). 

The duality is explicit in the labels at the top: the 
patterns are identically the same for SEnn (tendexes) of 
electric-parity modes and SBnn (vortexes) of magnetic- 
parity modes [left column] ; and also identically the same 
when the parities are switched [right column] The color 
coding is similar to Fig. [3] above (left-hand scale). The 
red (light gray) regions are stretching tendexes or coun- 
terclockwise vortexes (negative eigenvalues); the blue 
(dark gray), squeezing tendexes or clockwise vortexes 
(positive eigenvalues). 

For the Schwarzschild hole, the electric-parity ten- 
dex pattern and magnetic-parity vortex pattern (upper 
left) is that of the spherical harmonic F 22 (#,0), and the 
perturbative electric-parity vorticity and magnetic-parity 
tendicity vanish (upper right). 

For the rapidly spinning Kerr hole, the electric-parity 
tendexes and magnetic-parity vortexes (lower left) are 
concentrated more tightly around the plane of reflection 
symmetry than they are for the Schwarzschild hole, and 
are twisted; but their patterns are still predominantly 
Y 22 . And also for Kerr, the (perturbative) electric- 
parity vorticity and magnetic-parity tendicity have be- 
come nonzero (lower right), they appear to be pre- 
dominantly Y" 32 (#,(/>) in shape, they are much less con- 
centrated near the equator and somewhat weaker than 
the electric-parity tendicity and magnetic-parity vortic- 
ity (lower left). 



Equatorial-plane vortex and tendex lines, and 
vortexes and tendexes 



As for the weak-field, radiative sources of Paper I, so 
also here, the equatorial plane is an informative and sim- 
ple region in which to study the generation of gravita- 
tional waves. 

For the (2,2) modes that we are studying, the SEjk of 
an electric-parity perturbation and the SBjk for magnetic 
parity are symmetric about the equatorial plane. This 
restricts two sets of field lines (tendex lines for electric- 
parity SEjk] vortex lines for magnetic-parity SBjk) to 
lie in the plane and forces the third to be normal to 
the plane. By contrast, the electric-parity SBjk and 
magnetic-parity SEjk are reflection antisymmetric. This 
requires that two sets of field lines cross the equatorial 
plane at 45° angles, with equal and opposite eigenvalues 
(tendicities or vorticities), and forces the third set to lie 
in the plane and have zero eigenvalue (which makes them 
have less physical interest). 

In this section, we shall focus on the in-plane field lines 
and their vorticities and tendicities. 
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FIG. 9: (color online). Perturbative horizon tendicities SEnn 
and vorticities SBnn for the (2,2) quasinormal modes with 
electric and magnetic parities (see column labels at the top). 
The top row is for a Schwarzschild black hole, a = 0; the 
bottom for a rapidly spinning Kerr black hole, a/M — 0.945. 
The color intensity is proportional to the magnitude of the 
tendicity or vorticity, with blue (dark gray) for posi tive an d 
red (light gray) for negative. For discussion, see Sec. Ill A of 
the text. 



1. Magnetic-parity perturbations of Schwarzschild black 
holes 



In Sec. |IC 6| and Figs. [2] and [3j we discussed some 
equatorial-plane properties of the magnetic-parity (2, 2) 
mode. Here we shall explore these and other properties 
more deeply. Recall that for the magnetic-parity mode, 
the frame-drag field, and hence also the vortex lines and 
their vorticities, are fully gauge invariant. 

In Fig.[l0j we show six different depictions of the vortex 
lines and their vorticities in the equatorial plane, each 
designed to highlight particular issues. See the caption 
for details of what is depicted. 

The radial variation of vorticity is not shown in this 
figure, only the angular variation. The vorticity actually 
passes through a large range of values as a function of 
radius: from the horizon to roughly r = AM ~ 1.5A 
(roughly the outer edge of the near zone), the vorticity 
rapidly decreases; between r c± 4M and 12M (roughly 
the extent of the transition zone), it falls off as 1/r; and 
at r > 12M (the wave zone), it grows exponentially due 
to the damping of the quasinormal mode as time passes. 
(The wave field at larger radii was emitted earlier when 
the mode was stronger.) In the figure, we have removed 
these radial variations in order to highlight the angular 
variations. 



By comparing the left panels of Fig. 10 with Fig 9 of 
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FIG. 10: (color online). Three representations of the vortex lines and vortexes in the equatorial plane of a Schwarzschild 
black hole perturbed by a magnetic-parity (2, 2) quasinormal mode. The bottom panels span a region 56M on each side, 
and the top panels are a zoom-in of the lower panels, 30M on each side. All panels show positive-vorticity lines as solid and 
negative-vorticity lines as dashed. In all panels, blue (dark gray) corresponds to positive and red (light gray) to negative; the 
intensity of the color indicates the strength of the vorticity at that point normalized by the maximum of the vorticity at that 
radius (darker shading indicates a larger strength and lighter, weaker). Similarly, in all panels, the central circle surrounded 
by a narrow white line is the horizon colored by its vorticity as described above. Left column: Vortex lines colored and shaded 
by their scaled vorticity. Middle column: Negative vorticity coloring the plane with black vortex lines. Right column: Vorticity 
with the larger absolute value coloring the plane and black vortex lines. For discussion of this figure, see Sec. |III B 1| 



Sec. VI D of Paper I, we see a strong resemblance be- 
tween the vortex lines of our (2, 2), magnetic-parity per- 
turbation of a Schwarzschild black hole, and those of a 
rotating current quadrupole in linearized theory. As in 
linearized theory, when the radial (or, synonymously, lon- 
gitudinal) vortex lines in the near zone rotate, the effects 
of time retardation cause the lines, in the transition and 
wave zones, to collect around four backspiraling regions 
of strong vorticity (the vortexes) and to acquire pertur- 
bative tendex lines as they become transverse-traceless 
gravitational waves. The most important difference is 
that, for the black- hole perturbations, the positive vor- 
tex lines emerge from the blue, clockwise horizon vortexes 
and spiral outward (and the negative vortex lines emerge 
from the counterclockwise horizon vortexes) rather than 
emerging from a near-zone current quadrupole. 

Although the left panels of Fig. [10] highlight most 
clearly the comparison with figures in Paper I, the mid- 
dle and right panels more clearly show the relationship 
between the vortex lines (in black) and the vorticities, 
throughout the equatorial plane. In the middle panels 
(which show only the negative vorticity), the negative 



vortex lines that emerge longitudinally from the horizon 
stay in the center of their vortex in the near zone, and 
then collect onto the outer edge of the vortex in the tran- 
sition and wave zones. Interestingly, near the horizon, 
there are also two weaker regions of negative vorticity 
between the two counterclockwise vortexes, regions as- 
sociated with the tangential negative vortex lines that 
pass through this region without attaching to the hori- 
zon (and that presumably represent radiation traveling 
into the horizon). 



In the right panels of Fig. 10 (which show the in-plane 
vorticity with the larger absolute value), a clockwise vor- 
tex that extends radially from the horizon takes the place 
of the weaker region of counterclockwise vorticity. From 
these panels, it is most evident that the vortexes and vor- 
tex lines of opposite signs are identical, though rotated by 
90°. These panels also highlight that there are four spi- 
rals of nearly zero vorticity that separate the vortexes in 
the wave zone, which the spiraling vortex lines approach. 
All three vorticities nearly vanish at these spirals; in the 
limit of infinite radius, they become vanishing points for 
the radiation, which must exist for topological reasons 
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2. Gauge dependence of electric-parity tendexes for a 
Schwarzschild black hole 

In this subsection, we explore the gauge dependence 
of the (2,2) modes for a Schwarzschild black hole. In 
Sec. II C[ we showed that for magnetic-parity modes, all 
gauges share the same slicing. Therefore, to maximize 
any gauge dependence that there might be, we focus on 
the electric-parity (2,2) mode. 

Because the frame-drag field of the unperturbed 
Schwarzschild black hole vanishes, this mode's perturba- 
tive frame-drag field will be unaffected by perturbative 
changes of the spatial coordinates. Therefore, we focus on 
the perturbative tidal field 58 of the electric-parity mode, 
which is sensitive to both perturbative slicing changes and 
perturbative spatial coordinate changes. 

In Fig. [TT] we plot this field's perturbative equatorial 
tendexes and tendex lines for the electric-parity (2,2) 
mode in RWZ gauge (left panel) and IR gauge (right 
panel) — which differ, for this mode, in both slicing and 
spatial coordinates. The tendex lines for the two gauges 
were seeded at the same coordinate points, so all the 
differences between the panels can be attributed to the 
gauge differences. 

The two panels are almost identical. Therefore, these 
maximally sensitive tendexes and tendex lines are re- 
markably unaffected by switching from one gauge to the 
other. The primary differences are that (i) the tendex 
lines of IR gauge tend to be pulled closer to the horizon 
very near the horizon as compared to RWZ gauge (ii) the 
lines falling onto the attracting spiral are bunched even 
more tightly in IR gauge than in RWZ gauge; however, 
more lines reach the spiral in RWZ gauge in this figure, 
and (hi) the four tendex spirals wind more tightly in IR 
gauge, which is most easily seen by comparing the lower 
right and upper left corners of the two panels. 

One subtlety that must be remarked upon is that the 
central circle colored by the normal-normal component of 
the tidal field (surface tendicity) in the RWZ gauge (left 
panel of Fig. 11) is simply the surface r = 2M, and not 
the true event horizon. The location of the event hori- 
zon is affected by the perturbations in a gauge-dependent 
manner, as discussed by Vega, Poisson, and Massey [18]. 
We rely on the results of this article in the brief dis- 
cussion that follows. In RWZ gauge, the horizon is at 
th = 2M J r 5r{i 1 </>), where the function 5r can be solved 
for by ensuring that the vector tangent to the perturbed 
generators (in our case, / M = dx^/dt) remains null [18] . 
We give an expression for Sr in Appendix |A 5[ There 
we also discuss the correction to the horizon tendicity in 
RWZ gauge. One key result is that the horizon tendicity 
has the same angular distribution in RWZ gauge as in 
IR gauge (given by the Y 22 spherical harmonic), so that 
the normal-normal tendicity on the horizon and on the 
surface at r = 2M differ only by an amplitude and phase 



in RWZ gauge. Meanwhile, in IR gauge the horizon re- 
mains at th = 2M and so the colored central circle is in 
fact the horizon, colored by its horizon tendicity. 

However, the bulk tendexes and tendex lines are deter- 
mined completely independently of these horizon consid- 
erations, and so Fig. [TT] provides an accurate comparison 
of them in the two gauges. 



3. Duality and influence of spin in the equatorial plane 



In this subsection, we use Fig. 12 to explore duality and 
the influence of spin, for the fundamental (2, 2) mode. 

By comparing the left and center panels in the top row 



of Fig. 12, we see visually the near duality between the 
electric- and magnetic-parity modes for a Schwarzschild 
hole. This near duality is explored mathematically in 
App. [C] Specifically, the vortexes and their lines for the 
magnetic-parity mode (left) are nearly identical to the 
tendexes and their lines for the electric-parity mode. The 
only small differences appear in the size of the nearly 
zero-vorticity (or tendicity) regions, and the curvatures 
of the lines. 

For the fast-spinning Kerr black hole (the bottom left 
and center panels of Fig. 12), the near duality is still 



obvious, especially in the colored vortexes and tendexes; 
but it is less strong than for Schwarzschild, especially in 
the field lines. The vortex lines (on the left) continue 
to look like those of a Schwarzschild black hole, but the 
tendex lines (in the middle) curve in the opposite direc- 
tion, which makes some lines reach out from the horizon 
and connect back to it instead of spiraling away from the 
horizon. 

By comparing the top and bottom panels in the left 



and center columns of Fig. 12 , we see the influence of 
the background black hole's spin on the dynamics of the 
perturbative vortexes and tendexes. For fast Kerr (bot- 
tom), the vortexes and tendexes near the horizon look 
more transverse (less radial) than for Schwarzschild, be- 
cause the size of the near zone is much smaller. (The 
frequency of the waves is nearly twice that for a per- 
turbed Schwarzschild hole.) The higher frequency also 
explains why the spirals of the vortexes and tendexes are 
tighter. 

In the isometric embedding diagrams in the right col- 
umn of Fig. [12] we see that proper radial distance in 
the near zone is somewhat larger than it appears in the 
flat, planar drawing. Taking this into account, we con- 
clude that, aside from a few small differences, the qual- 
itative ways in which waves are generated for fast Kerr 
and for Schwarzschild are the same: two pairs of vortexes 
or tendexes emerge longitudinally from horizon vortexes, 
and twist into backward spirals that eventually form the 
transverse-traceless gravitational waves. 
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FIG. 11: (color online). The equatorial-plane, electric-parity tendexes and tendex lines of a (2, 2) perturbation of a Schwarzschild 
black hole in RWZ gauge (left panel) and IR gauge (right panel). The conventions for the lines, the coloring and the shading 
are identical to those in the right panels of Fig. [TU] 



4- Vortexes of electric-parity mode, and perturbative 
tendexes of magnetic-parity mode for a Schwarzschild black 

hole 



In Fig. [l3j we visualize the vortexes of the electric- 
parity (2,2) mode of a Schwarzschild hole. (By near du- 
ality, the perturbative tendexes of the magnetic-parity 
mode must look nearly the same.) 

As noted at the beginning of Sec. |IIIB} reflection anti- 
symmetry of the frame-drag field for this electric-parity 
mode dictates that through each point in the plane there 
will pass one zero-vorticity vortex line lying in the plane, 
and two vortex lines with equal and opposite vorticities 
that pass through the plane at 45 degree angles and are 
orthogonal to each other and to the zero-vorticity line. 



The vorticity plotted in Fig. [13] is that of the counter- 
clockwise, 45 degree line. For the clockwise line, the vor- 
ticity pattern is identically the same, but blue instead of 
red. 

There are again four regions of strong vorticity (four 
vortexes), which spiral outward from the horizon, becom- 
ing gravitational waves. In this case, the four regions look 
identical, whereas for the tendexes of this same electric- 



parity mode (middle column of Fig. 12) there is an al- 



ternation between blue and red. There is actually an 
alternation here, too, though it does not show in the fig- 
ure: The relative tilt of the lines (in the sense of the <ft 
direction) rotates, such that in one tendex, the red ten- 
dex lines pass through the plane with a forward 45 degree 
tilt on average, and in the next tendex, with a backward 
45 degree tilt; and conversely for the blue tendex lines. 



C. Three-Dimensional vortexes and tendexes 

In this section, we shall explore the 3-dimensional vor- 
texes and tendexes of the (2,2) modes of Schwarzschild 
and Kerr black holes, which are depicted in Fig. [l4j In 
the first subsection, we shall focus on what this figure 
tells us about the generation of gravitational waves, and 
in the second, what it tells us about duality. 



1. Physical description of gravitational- wave generation 



In Sec. II C 61 of the Introduction, we summarized in 



great detail what we have learned about gravitational- 
wave generation from our vortex and tendex studies. 
There we focused on the (2,2) magnetic-parity mode, 
and among other things we scrutinized the upper left 
panel of Fig. 14 (which we reproduced as Fig. [4|. Here, 



instead, we shall focus on the (2, 2) electric-parity mode 
as depicted in the right half of Fig. [14] 

We begin with the perturbative tendexes of the 
electric-parity, (2,2) mode of a Schwarzschild black hole 
(third panel on top row of Fig. 14). The 3-D tendexes 



emerge from the horizon as four deformed-cylinder struc- 
tures, two red (light gray) and two blue (dark gray). 
These are the extensions into the third dimension of the 
four near-zone, equatorial-plane tendexes of the center 



panel in Fig. 12 above. As we enter the transition zone, 
the four 3D tendexes lengthen vertically (parallel to the 
poles), and then as we enter the wave zone, they spi- 
ral upward and downward around the poles; they have 
become gravitational waves. They are concentrated near 
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FIG. 12: (color online). Vortexes and tendexes and their field lines in the equatorial plane for (2, 2) modes of Schwarzschild and 
Kerr black holes. The lines, the coloring and the shading are identical to those in the right panels of Fig. |10[ The upper row is 
for a Schwarzschild black hole (a = 0); the lower row, for a rapidly spinning Kerr black hole (a = 0.945); see labels on the left. 
The left column shows the vortex lines and vorticities for magnetic parity (which are gauge invariant for the perturbations of 
a Schwarzschild hole); the middle and right columns show the tendex lines and tendicities for the electric-parity mode in IR 
gauge; see labels at the top. In the right column, the equatorial plane is isometrically embedded in three-dimensional Euclidean 
space. The top panels are 24M across; the bottom, 14M. This figure elucidates duality and the influence of black-hole spin; 
see the discussion in Sec. IIII B 31 




FIG. 13: (color online). For an electric-parity (2,2) mode of 
a Schwarzschild black hole: the vorticity of the counterclock- 
wise vortex lines that pass through the equatorial plane at a 
45 degree angle. The clockwise vortex lines that pass through 
the plane have equal and opposite vorticity. By near duality, 
this figure also depicts the perturbative tendex structure for 
the magnetic-parity (2,2) mode. The conventions for coloring 
and shading are the same as in Fig. |10| Because the horizon 
vorticity is exactly zero for this mode, the horizon is shown as 
a white disk. The left panel, a region 30M across, is a zoom 
in of the right panel, which is 56M across. 



the poles because the (2,2)-mode gravitational waves are 
significantly stronger in polar directions than in the equa- 
tor. 

In this panel, we also see black tendex lines that emerge 
from the horizon and spiral upward and downward along- 
side the polar-spiraling tendexes, becoming nearly trans- 
verse at large radii — part of the outgoing gravitational 
waves. Of course, there are similar tendex lines, not 
shown, inside the spiraling tendexes. In addition, we also 
see tendex lines in the inner part of the wave zone that 
are approximately polar circles; these are also part of the 
outgoing waves. 

The top rightmost panel depicts the vortexes associ- 
ated with this electric-parity mode. The horizon vortic- 
ity vanishes, so the horizon is white. The vortexes near 
the horizon are dominated by the longitudinal-transverse 
part of the frame-drag field B LT , which interacts with 
S£ L and 5£ LT to maintain their joint near-zone structure 
as they rotate (cf. th e desc ription of the dual magnetic- 
parity mode in Sec. 



is also a B TT associated with the ingoing gravitational 



IC6). However, of course, there 
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FIG. 14: (color online). Three-dimensional vortexes and tendexes of the same four modes as are shown in Fig. |12| As there, 
so here , th e top row is for a Schwarzschild black hole, and the bottom for fast-spin Kerr, a — 0.945; the left two columns (one 
in Fig. 12) are for a magnetic-parity (2,2) mode (with the vortexes in IR coordinates), and the right two columns (one in Fig. 



12) are for electric parity in IR gauge. For each parity, the first column shows structures of the field that generates the waves 
JSB for magnetic parity; 58 for electric parity) and the second column shows structures of the other field (not included in the 
equatorial-plane drawings of Fig. 12). In each panel, the colored surfaces show the outer faces of vortexes (for SB) or tendexes 
(for SE), defined as the locations, for a given radius, where the largest-in-magnitude eigenvalue of the field being plotted (SB or 
SE) has dropped to a certain percentage (90, 85, 80 or 75) of its maximum for that radius; that percentage is shown alongside 
the colored surfaces. As in previous figures, the surface is red (light gray) if that largest-in-magnitude eigenvalue is negative 
and blue (dark gray) if positive. The off-white regions are surfaces where the largest-in-magnitude eigenvalue has dropped to 
15%, 20% or 25% of the maximum at that radius. In each panel the black lines are a few of the vortex lines (for SB panels) or 
tendex lines (for SE panels) that become transverse when they reach large radii, and thereby produce the tidal or frame-drag 
force of an emitted gravitational wave. For discussion of this figure, see Sec. |III C| 



waves. At large radii, in the out going- wave zone, the vor- 
texes, like the tendexes of the third panel top row, spiral 
upward and downward around the polar axis; they have 
joined with the tendexes to form the full gravitational- 
wave structure. 

For insight into how (we think) the near-zone tendexes 
of this electric-parity mode, extending radially out of the 
horizon, generate these outgoing gravitational waves, and 
how the ingoing waves, that they also generate, act back 
on them and drive their gradual decay, see the description 
of this mode's dual in Sec. II C 61 

For the rapid- spin Kerr black hole, the tendex and vor- 
tex structures (last two panels of second row of Fig. 14 ) 



are quite similar to those for the Schwarzschild black hole. 
The detailed differences are similar to those in the equa- 
torial plane (see discussion in Sec. IIIB3 above): smaller 
near zone and tighter spiraling for the tendexes because 
of the higher eigenfrequency; nonvanishing horizon vor- 
ticity with a predominantly F 32 (#,(/>) angular structure. 
In the near zone, the 3D vortexes seem to have acquired 
a longitudinal (radial) part, emerging from the Y 32 hori- 
zon vortexes (though this is largely hidden behind the 



off-white structures). Thus, for a Kerr black hole, one 
might intuitively describe the generation of gravitational 
waves as being produced by a superposition of near-zone 
tendexes that induce vortexes by their motions, and near- 
zone vortexes that induce tendexes by their motions. 
However, because the near-zone vortexes are weaker than 
the tendexes, the tendexes still play the dominant role 
for gravitational- wave generation in this electric-parity 
mode. 



2. Approximate duality 

By comparing the magnetic-parity left half of Fig. [14] 
with the electric-parity right half, we can visually as- 
sess the degree to which there is a duality between the 
modes in three dimensions. For the perturbations of 
Schwarzschild black holes (top row), the most notable 
difference between the magnetic- and electric-parity per- 
turbations is that the transition between the longitudi- 
nal near-zone and spiraling wave-zone vortexes of the 
magnetic-parity perturbation is more abrupt, and hap- 
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pens closer to the horizon than it does in the electric- 
parity perturbations. The re ason for this is encoded in 
Eqs. (C29), (C30), and (C31), but we do not have a sim- 



ple physical explanation for why this occurs. This differ- 
ence is magnified for perturbations of the rapidly rotating 
Kerr black hole (bottom row). Thus, the small breaking 
of duality quantified in App. [C] for Schwarzschild black 
holes seems to be more pronounced in three-dimensions 
than in two, and stronger for rapidly rotating black holes 
than for non-rotating ones. 

Nevertheless, the qualitative picture of wave genera- 
tion by longitudinal near-zone tendexes and vortexes is 
essentially dual for perturbations of the two parities. 



D. Comparing vortex lines of a perturbed Kerr 
black hole and a binary-black-hole-merger remnant 



As a conclusion to this section and a prelude to future 



work, in Fig. 15 we compare the vortex lines found us- 
ing our analytic methods to those found in a numerical 
ringdown of a fast-spinning Kerr black hole. 

More specifically, we compare an electric-parity, (2,2) 
quasinormal-mode perturbation of a Kerr black hole with 
dimensionless spin a/M = 0.945, to a ringing-down Kerr 
black hole of the same spin formed in a numerical simula- 
tion [19] of the merger of two equal-mass black holes with 
equal spins of magnitude 0.97 aligned with the orbital 
angular momentum. (Note that because of the symme- 
try of this configuration, during the ringdown there is no 
magnetic-parity (2, 2) mode excited.) 

For both the analytical and numerical calculations, the 
vortex lines are those of the full frame-drag tensor B^ = 

B 



(0) 



B$ , including the background part B^. J and the 



2(0) 



SBtf (as in App. 



F) 



perturbation B~~ _ ~ ~ , , 

lines look quite different from those for just SB, 



This makes the 
depicted 



in Figs. 12 and 14 



Near the equator, the vortex lines in both panels look 
like those of an unperturbed Kerr black hole (see Paper 
II). Closer to the axis of rotation, the background vortex 
lines become degenerate, and the perturbations break the 
degeneracy by picking the principal axes of the pertur- 
bative field. Correspondingly, near the rotation axis and 
at large enough radii to be in or near the wave zone, the 
vortex lines resemble those of transverse-traceless grav- 
itational waves, which are emitted symmetrically above 
and below the hole. 

Although the vortex lines from these similar physical 
situations were computed using very different methods 
and gauge conditions, the results are remarkably similar 
(see Fig. 15). This gives us confidence that our analytical 
methods can guide our understanding of the vortexes and 
tendexes in the late stages of numerical simulations. 



IV. SUPERPOSED (2, 2) AND (2, -2) 
QUASINORMAL MODES OF SCHWARZSCHILD 

A. Magnetic-parity superposed modes 



In Sec. |I C 9[ we summarized the properties of the 
quasinormal mode of Schwarzschild that is obtained 
by superposing the magnetic-parity (2, 2) and (2, —2) 
modes. Here we give details. The vortex lines and vor- 
ticities for this superposed mode are depicted in Fig. 16 
using the three types of visualizations in Fig. fTo] 



The left column of Fig. 16 shows the two families of 
vortex lines that lie in the equatorial plane, color coded 
by their vorticities. The solid-line family has predomi- 
nantly positive (clockwise) vorticity, but in some regions 
its vorticity becomes weakly negative (counterclockwise) . 
The dashed-line family has predominantly negative vor- 
ticity, but in some regions it is weakly positive. A rota- 
tion around the hole by angle tt/2 maps each family into 
the other. 

In the center column of Fig. [l6j the vortex lines are 
drawn black and the equatorial plane is colored by the 
vorticity of the dashed lines. To deduce the coloring for 
the solid lines, just rotate the colored plane (but not the 
lines) by tt/2 and interchange red (light gray) and blue 
(dark gray). By contrast with most previous figures, the 
radial variation of the vorticity is not scaled out of this 
figure; so in the wave zone (roughly, the outer half of 
right panel) the coloring oscillates radially, in color and 
intensity, in the manner of a gravitational wave. At large 
radii, there is also a growth of intensity (and saturation 
of the color scale) due to the waves emitted earlier having 
larger amplitude. 

In the right column of Fig.[l6j the vortex lines are again 
drawn black, and the equatorial plane is now colored by 
the larger of the two vorticities in amplitude. 



Together, the columns of Fig. 16 provide the following 
picture: For each family of lines, the equatorial-plane vor- 
texes form interleaved rings (dashed lines and red [light 
gray] vortexes for center column). Most of the family's 
vortex lines form closed, distorted ellipses that, when tan- 
gential, lie in a single vortex (red for dashed lines), and 
when more nearly radial, travel from one vortex to an- 
other. In the wave zone, these line and vortex structures 
grow longer tangent ially as they propagate outward, and 
they maintain fixed radial thickness. When one looks at 
both families simultaneously, focusing on the strongest 
at each point (right column), one sees vortexes of alter- 
nating red and blue vorticity (light and dark gray). The 
angular oscillations are those of a quadrupolar structure; 
the radial oscillations are those of a propagating wave. 

In the near zone, the vortex lines have rather sharp, 
right angled features associated with the quadrupolar na- 
ture of the near-zone perturbation; this is to be compared 
to the oscillating current quadrupole in linearized gravity 
(Paper I, Sec. VI, Fig. 15). There are multiple singular 
points in the lines, degenerate points where three lines 
cross with sharp bends, and where both families of vortex 
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FIG. 15: (color online). Left: Vortex lines of a a/M = 0.945 Kerr black hole perturbed by an electric-parity (2, 2) quasinormal- 
mode in IR gauge. Right: Vortex lines from a a/M = 0.945 ringing down Kerr black hole obtained from a numerical simula- 
tion [19] of two identical merging black holes with spins of magnitude 0.97 aligned parallel to the orbital angular momentum. 
In the simulation, we chose a late enough time that the common apparent horizon is essentially that of a single, perturbed 
black hole, and we computed the vortex lines using methods summarized in Ref. pQ. 



lines take on the same eigenvalue. (The third eigenvalue, 
that of the lines perpendicular to the plane, must then 
be minus twice the vorticity of these lines, in order for 
the sum of the eigenvalues to vanish). 

We shall discuss the dynamics of these vortex lines 
and vortexes in Sec. |IV Cj below, after first gaining insight 
into the electric-parity superposed mode (whose vortexes 
will teach us about this magnetic-parity mode's tendexes 
through the near duality). 



B. Electric-parity superposed mode 



texes have transitioned into outward traveling transverse 
waves, with regions of vanishing vorticity between the 
crests and troughs of each wave. The waves are strongest 
along the diagonals, though in the near zone the (LT) 
frame-drag field is strongest in the up, down, left and 
right directions. 

By (near) duality, the tendexes of the magnetic-parity 
superposed mode will have the same form as these 
electric-parity-mode vortexes. Accordingly, in the next 
section, we will use this figure to elucidate the magnetic- 
parity mode's dynamics — and by duality, also the dy- 
namics of this electric-parity mode. 



For the mode constructed by superposing electric- 
parity (2,2) and (2,-2) modes of Schwarzschild, as for the 
electric-parity (2,2) mode itself (Sec. IIIB4), symmetry 
considerations dictate that: (i) one family of vortex lines 
lies in the equatorial plane and has vanishing vorticity, 
(ii) two families pass through the equatorial plane at 45° 
angles, with equal and opposite vorticities, and (iii) the 
horizon vorticity vanishes. 



In Fig. 17 [analog of Fig. 13 for the (2,2) electric-parity 
mode], we show the vorticity of the family of counter- 
clockwise vortex lines, as they pass through the equato- 
rial plane. 

In the near zone of this figure, we see again a distinct 
quadrupolar structure, with four lobes of strong vorticity 
present near the horizon (four near-zone vortexes). Be- 
yond these near-zone lobes, there is a ring of vanishing 
vorticity, followed by an annulus where the cast-off vor- 
texes of a previous cycle have begun to deform into an 
annulus of stronger vorticity. In the wave zone, the vor- 



Dynamics of the magnetic-parity superposed 
mode 



We now turn to the dynamics of the magnetic-parity 
superposed mode, which we studied in Sec. |IV A| 

In Fig. 18 for this mode we show, in the equatorial 
plane, the time evolution of (i) the vortex lines and their 
vorticities near the black hole (top row), and (ii) on a 
larger scale that extends into the wave zone, the mode's 
vortexes (middle row) and perturbative tendexes (bot- 
tom row). The four panels in the top row are stills from 
a movie at [20]. To be absolutely clear, the vortexes and 
tendexes are those of the same magnetic-parity mode. 
As in the center column of Fig. [l6j the top and middle 
rows show only one family of vortexes, that for the dashed 
vortex lines which have predominantly negative vorticity; 
and as in Fig. [TtJ the bottom row shows only the tendic- 
ity of the negative-tendicity perturbative tendexes that 
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FIG. 16: (color online). The vorticities and vortex lines in the equatorial plane of a Schwarzschild black hole, for the fundamental 
magnetic-parity (2, 2) mode superposed on the fundamental magnetic-parity (2, —2) mode, depicted using the same three 
visualization techniques as in Fig. |1Q| Here, however, we do not scale the vorticity by any function, but the numbers on the 
vorticity scale on the left of the panels are equal to v / vorticrt}0^7 : (where r is radius), in units of the maximum value of this 
quantity, which occurs on the horizon at <j) — 37r/4 and <f> — 77r/4. The top panels cover a region 30M across, and the bottom 
panels are a zoom-out of the upper panels, 56M across. The central circle in all panels is the horizon as viewed from the polar 
axis, colored by its vorticity. Left column: The two families of vortex lines (one shown dashed, the other solid) with each line 
colored, at each point, by the sign of its vorticity (blue [dark gray] for positive, i.e., clockwise; red (light gray) for negative, 
i.e., counterclockwise), and each line has an intensity proportional to the magnitude of its vorticity. Center column: The same 
vortex lines are colored black, and the equatorial plane is colored by the vorticity of the dashed family of lines. Right column: 
The same as the center column, but the equatorial plane is colored by the vorticity with the larger magnitude. 



pass through the equatorial plane at 45°. Time t = (left 
panels of Fig. 18 ) is chosen at a moment when the horizon 



vorticity is maximum, whereas Figs. [16] and [17] are snap- 
shots at the slightly earlier time at ~ — tt/3 (which gives 
nearly the same vortex structures as the fourth column, 
after a rotation by tt/2). 

In interpreting Fig. [THJ especially the top row, we em- 
phasize that there is no unique way of following a sin- 
gle vortex or tendex line in time. The same is true of 
electric and magnetic field lines in Maxwell's theory (cf. 
[21]). While we hope to elucidate this issue in future 
work, here, in constructing the panels in the top row of 



Fig. 18, we have simply started the integration of the vor- 
tex lines from the same points at each time step, making 
no attempt to identify and follow individual lines from 
moment to moment. Correspondingly, in order to inter- 



pret Fig. 18 and gain insight into the dynamics of the 
superposed mode, instead of trying to follow individual 
lines, we will focus on the lines' evolving shapes, and the 
structures of the vortexes and tendexes and the equations 



governing their evolution on the horizon. 

As a foundation for understanding the near-zone dy- 
namics depicted in this figure, we write down explicit ex- 
pressions for the longitudinal and longitudinal-transverse 
parts of the frame-drag and tidal fields on the horizon: 
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Here the normalization is that of App. [A} lo = 
(0.37367— 0.08896i) /M is the mode's eigenfrequency, e^ B 
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FIG. 17: (color online). For the electric-parity, superposed 
(2,2) and (2,-2) fundamental modes of Schwarzschild: the vor- 
ticity of the counterclockwise vortex lines that pass through 
the equatorial plane at 45° angles. By near duality, this 
figure also depicts the perturbative tendex structure for the 
magnetic-parity superposed mode. The intensity scale of the 
red color (left edge of figure) is the same as that in the center 
column of Fig. |16| The left panel, a region 30M across, is a 
zoom-in of the right panel, which is 56M across. This figure 
is the superposed-mode analog of Fig. |13| 



is the Levi-Civita tensor on the horizon, 

y = Y 22 + Y 2 ~ 2 = V15/87T sin 2 cos 20 (4.1e) 

is this mode's scalar spherical harmonic, and D A is the 
covariant derivative on the unit 2-sphere (related to the 



covariant derivative on the horizon by D_% = 2MV 



C26a), and definition (2.1) of the EF time coordi 



A)' 



[Equations (|4lj) follow from Eqs. flA12| ), flA13| ), flA26|) of 
App. 
and ( 
nate. 



A| the vector-spherical-harmonic definitions (C22a) 



Equations (4.1) are the fields measured by Eddington- 
Finkelstein observers. The conservation law ( 1.4a| ) for 
longitudinal field lines threading the horizon (which we 
shall need below) involves, by contrast, the LT frame- 
drag field measured by Schwarzschild observers on the 
"stretched horizon" (very close to the event horizon). 
Since the Schwarzschild observers move outward with ve- 
locity v = (2M/r)N ~ N with respect to the EF ob- 
servers, and with 7 = — v 2 ~ l/(y/2a), the field 
they measure is B^ h = 7 (B LT - N x 5£ LT ). This field 
diverges as 1/a as the stretched horizon is pushed to- 
ward the event horizon; to remove that divergence, in the 
Membrane Paradigm [9]) we renormalize by multiplying 
with a: 
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(4.2) 



The s econ d line is obtained by inserting the EF fields 
( 4.1a ) and ( 4.1b ), and 7 = l/y/2a, into the first line. The 
conservation law for longitudinal vortex lines threading 
the horizon (actually, one of the Maxwell-like Bianchi 
identities in disguise) says that 



dB NN /dt + V A (-B? A )=0; 



(4.3) 



cf. Eq. (1.4a) and subsequent discussion. (In this 
Schwarzschild-perturbation-theory case, there are no 
small spin-coefficient terms to spoil the perfection of the 
conservation law.) The vortex-line density and flux ex- 



pressions (4.1a) and (4.2) do, indeed, satisfy this conser- 



vation law, by virtue of the fact that the 2-dimensional 
Laplacian acting on the quadrupolar spherical harmonic 

y gives D A D A y = -ey. 

Equations (4.1) and (4.3) tell us the follow ing: ( i) On 



and near t he hor izon, the LT fields B L T [E q. (4.1b)] and 



SS LT [Eq. kick], and also Bf A [Eq. (|4.2|)], all oscillate 
approxi mately out of phase with the longitudinal field 
B L [Eq. (4.1a)]. 3 Therefore, near-zone energy is fed back 
and forth between the L and LT fields as the black hole 
pulsates. 



(ii) The conservation law (4.3) says that, if we regard 



Sat at as the density of vortex lines of B threading the 
horizon, and as the flux of vortex lines (number 

crossing a unit length in the horizon per unit time), then 
these horizon-threading vortex lines of B L are conserved 
during the pulsation. More specifically: 

(iii) As the mode evolves in Fig. 18 from at = to 
at = 7r, the conserved vortex lines are pushed away from 
the center of each horizon vortex toward its white edges, 
and there the conserved lines from the red region (coun- 
terclockwise) annihilate with the conserved lines from the 
blue region (clockwise). The pushing is embodied in the 
vortex-line flux — which grows stronger during this 
evolution. 

(iv) As the mode evolves further from at = tt/2 to 
at = 7r, conserved vortex lines of B L are created in pairs 
(one clockwise, the other counterclockwise) at the white 
edges of the horizon vortexes, and move inward toward 
the center of each vortex. 

Turn, now, from the conserved vortex lines of B L pierc- 
ing the horizon to the 3D vortex lines outside the horizon, 
depicted in the top row of Fig. [18] Because these are lines 
of the full 3D frame-drag field B and not its longitudinal 
part B L , they do not obey a conservation law and there 
is no unique way of following individual lines from one 
panel to the next. However, their evolving shapes teach 
us much about the geometrodynamics of this superposed 
mode: 

At time at = (upper left panel), the horizon-piercing 
vortex lines of the full frame-drag field B are almost per- 
fectly radial, with clockwise (solid) tendex lines emerg- 



3 The longitudinal-transverse frame-drag field H LT lags approx- 
imately 1.04 « 7r/3 radians behind on the horizon, while 
the LT component of the tidal field, £ LT , lags approximately 
1.21 w 2tt/5 radians behind i3/viv- Most importantly for inter- 
preting Fig. |18| the two nonzero tendicities of the tidal field are 
(in the equatorial plane) + £|7, and on the horizon they 

lag nearly tt/2 radians behind Bnn; the damping of the per- 
turbation adds slightly to the phase lag, so that it is actually 
7r/2 + arctan[S(w)/sft(w)] ~ tt/2 + 0.234 out of phase with B NN . 
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FIG. 18: (color online). Time evolution of the equatorial vortexes (top and middle rows) and equatorial perturbative tendexes 
(bottom row) for the superposed (2, 2) and (2, —2) magnetic-parity mode of Schwarzschild in RWZ gauge. The color scale is the 
same as the center column of Fig. |16| and the gravitational- wave-induced exponential decay of the vorticity and tendicity has 
been removed. Top row: Equatorial vortex lines and their vorticity plotted in a region near the horizon (16M across) followed 
over time t. The real part of the eigenfrequency is denoted a, so the successive panels, left to right cover half a cycle of the 
mode's oscillation. Middle row: The vorticity of the equatorial vortex lines in the near, intermediate and beginning of wave 
zone (30M across) at the same time steps as the top row. Bottom row: Tendicity of the counterclockwise tendex lines passing 
through the equatorial plane (which is dual to the left panel of Fig. 17), plotted at the same time steps as the top row. 



ing from the two blue horizon tendexes, and counter- 
clockwise (dashed) tendex lines emerging from the two 
red horizon tendexes. As time passes, the horizon pierc- 
ing lines become less radial and the horizon vorticity de- 
creases (at = 7r/4) until the lines' angles to the horizon 
are almost all near 45° and the horizon vorticity vanishes 
(at = 7r/2). Note that the lines that lie precisely on the 
diagonals, and which contact the horizon radially in the 
middle panel (at = 7r/2), have zero vorticity where they 
strike the horizon. This latter fact allows them to have 
a more radial angle of intersection than almost all other 
lines. The near-horizon frame-drag field has evolved at 
this time from being predominantly longitudinal, B L , to 
being predominantly longitudinal-transverse, B LT , but 
with some small admixture of transverse-traceless ingo- 
ing waves, B TT . As time moves onward from at = ir/2 
to at = 7r, the vortex lines in the outer part of each panel 
reach around on the horizon and attach to a quadrant on 
the side rather than directly below themselves — a quad- 
rant that has newly acquired the color corresponding to 
the lines' own vorticity (blue for solid lines, red for dashed 
lines). 

At at = (upper left panel), the near-horizon, nearly 



circular vortex lines in each quadrant represent, predom- 
inantly, the transverse-isotropic part of B L and keep it 
trace free. As time passes and B L decreases, we can re- 
gard these lines as traveling outward, forming the dis- 
torted ellipses which become gravitational waves in the 
far zone. The manner in which these circular vortex 
lines are restored each cycle appears to be as follows: 
As discussed above, as the horizon vorticity oscillates 
through zero, formerly longitudinal lines are pushed away 
from their respective vortexes and become first more 
longitudinal- trans verse, and then attach to a different 
quadrant of the horizon; meanwhile the B L is being re- 
generated with opposite sign in each quadrant, which re- 
quires new transverse-isotropic lines of opposite vorticity. 
These lines run tangent to the horizon, hugging it while 
they cross through a vortex of strong opposite vorticity. 
At each edge of the vortex these isotropic-transverse lines 
link up with predominantly radial lines which have the 
same sign of vorticity (it appears that it is the degenerate 
points at these edges that allow for such a reconnection) . 
This deforms the highly distorted, nearly circular arcs, 
which then lift off the horizon and propagate away as the 
cycle progresses. 
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This entire evolution is being driven by the oscillatory 
turn-off and turn-on of the longitudinal part of the frame- 
drag field B L . 

Turn attention to the evolution of the mode's equato- 
rial vortexes (middle row) and tendexes (bottom row). In 
accord with our dis cussi on above of the evolution of the 
horizon fields, Eqs. ( |4.1| ), these panels reveal (see below) 
that the vortexes and tendexes oscillate out of phas e with 
each other. Near-zone energy (see footnote 1 in Sec. I C 8 ) 
gets fed back and forth between vortexes and tendexes in 
an oscillatory manner (though during this feeding, some 
of it leaks out into the transition zone and thence into 
gravitational waves). This oscillatory feeding enables the 
near-zone perturbative tendexes to store half of the oscil- 
lation energy (while the LT vortexes store the other half) 
when the longitudinal frame-drag field B L is temporarily 
zero; and then use that energy to regenerate B L . 



The evidence for this near-zone feeding, in Fig. 18 



the following: (i) The near-horizon vortexes are strongest 
along the diagonals, while the regions of strong near-zone 
tendicity always occur along the vertical and horizontal 
directions. Thus, the vortexes and tendexes tend to oc- 
cupy different regions, with a 7r/4 rotation between the 
patterns [as on e shou ld ex pect f rom the angular depen- 
dences in Eqs. (4.1a) and ( 4.1c[ )]. (ii) There is a c± tt/2 
phase difference in the time evolution of the vortexes 
and tendexes. At those times when the horizon vorticity 
and near- horizon vortexes are strongest, the near- horizon 
tendexes are weak. As the horizon oscillates through 
zero vorticity, the tendexes are reaching their maximum 
strength. 

A careful study of the phases of these time behaviors 
reveals that the dynamics are not precisely tt/2 out of 
phase, as can be seen clearly in the first panel in the 



bottom row of Fig. 18 though the horizon vorticity is at 
it's maximum, the tendicities have just oscillated through 
zero in this region and are beginning to regenerate. As 
mentioned in footnote 3 above, this additional phase lag 
is due to the mode's damping, and in radians its magni- 
tude is arctan[3(cj)/9?(cj)] ~ 0.234. 

The 7r/4 differences in spatial phase and c± tt/2 differ- 
ences in temporal phase are lost as the frame-drag and 
tidal fields travel outward through the intermediate zone 
and into the wave zone; an inspection of the outer edges 
of the time series plots shows bands of strong tendicity 
and vorticity in phase in time and space, propagating 
outward in synch. This must be the case, since for plane 
waves in linearized gravity, the vortex and tendex lines 
are in phase temporally and spatially (though the lines 
are rotated by 7r/4 with respect to each other at each 
event; see Paper I, Sec VI A). 

By scrutinizing the middle and bottom rows of Fig. [18} 
(which extend from the near zone through the interme- 
diate zone and into the inner parts of the wave zone), 
one can see visually how the oscillatory feeding of en- 
ergy between near-zone vortexes and tendexes gives rise 
to outgoing vortexes and tendexes that represent gravi- 
tational waves. 



V. (2,1) AND (2,0) QUASINORMAL MODES OF 
SCHWARZSCHILD 



In this section we will complete our study of the 
quadrupolar perturbations of Schwarzschild black holes. 
Specifically, we will explore the vortex and tendex struc- 
tures and the dynamics of the (2,1) and (2, 0) magnetic- 
and electric-parity perturbations of a Schwarzschild hole, 
in RWZ gauge. 



A. Vortexes of (2, 1) magnetic-parity mode and 
perturbative tendexes of (2, 1) electric-parity mode 



In Sec. |I C 7\ we summarized the most important 
properties of the (2,1) magnetic-parity mode of a 
Schwarzschild black hole. In this subsection and the next, 
we shall give additional details about this mode and its 
electric-parity dual. We begin with the vortex structure 
for magnetic parity. 

The horizon vorticity of the magnetic-parity (2,1) 
mode has an angular dependence given by the spheri- 
cal harmonic Y 21 (6, <j>) (of course). We display this hori- 
zon vorticity in the left panel of Fig. 19 There are four 
horizon vortexes, two of each sign, and vanishing horizon 
vorticity all along the equator. 

As we noted in Sec. I C 7[ this mode's symmetry dic- 
tates that the frame-drag field be reflection antisymmet- 
ric through the equatorial plane. As for the electric- 
parity (2,2) frame drag field, which also has this prop- 
erty (second paragraph of Sec. IIIB also Sec. IIIB4 and 
Fig. [l3| ), this implies that one family of vortex lines lies 
in the equatorial plane with vanishing vorticity, and two 
cross through that plane at 45° with equal and opposite 
vorticities. The negative vorticities of the crossing lines 



are plotted in the middle panel of Fig. [19j along with 
the projected horizon vorticity, as if the horizon were 
viewed from above. The positive- vorticity pattern of the 
other family of crossing lines is identical to this negative- 
vorticity pattern, since at each point the two lines have 
the same vorticity magnitude. 

The fact that there are just two spiraling vortexes in 
this figure, by contrast with four for the (2,2) modes, is 
guaranteed by the modes' azimuthal orders, m = 1 here 
and m = 2 for (2,2). 

The vortex structure outside the equatorial plane, de- 
picted in the right panel of Fig. [19} was discussed in Sec. 
I C 7 The two red (light gray), 3D vortexes are the same 
ones depicted in the middle panel. They actually extend 
across the equatorial plane (via the 45° vortex lines) into 
the region occupied by the blue (dark gray) vortexes; but 
we do not see them there in the 3D drawing because the 
blue vortexes have larger vorticity and we have chosen to 
show at each point only the largest-vorticity vortex. 
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FIG. 19: (color online). The (2,1) magnetic-parity horizon vorticity and vortexes. Left panel: The horizon vorticity for the 
(2, 1) magnetic-parity perturbation of Schwarzschild, colored as in Fig. [9] There are four horizon vortexes, two clockwise (blue 
[dark gray]) and two counterclockwise (red [light gray]). The horizon vorticity vanishes at the equator and the poles. Middle 
panel: Vorticity of the counterclockwise vortex lines passing through the equatorial plane, colored and normalized as in the 
middle and right hand columns of Fig. 10 and plotted in a region 24M across. The (blue [dark gray]) vorticity of the clockwise 
vortex lines has precisely this same pattern, because the two families of lines pass through the equatorial plane with the same 
magnitude of vorticity at each point. Right panel: Three dimensional vortexes colored and labeled as in Fig. [14] By near 
duality, this figure also represents (to good accuracy) the tendicity and tendex structure of the (2,1) electric-parity mode. 



B. Vortexes of (2, 1) electric-parity mode and 
perturbative tendexes of (2, 1) magnetic-parity mode 

Turn, next, to the vortex lines and vortex structure of 
the (2,1) electric-parity mode. [By near duality, the per- 
turbative tendex lines and tendex structure of the (2,1) 
magnetic-parity mode will be the same.] 

For this mode, with the parity reversed from the pre- 
vious section, the frame-drag field is symmetric under 
reflection through the equatorial plane rather than anti- 
symmetric. Therefore, there are two sets of vortex lines 
that remain in the equatorial plane, with the third set 
normal to it. In this sense the vortexes structures are 
analogous to those of the magnetic-parity, (2, 2) mode; 
and in fact they are strikingly similar, aside from having 
two arms rather than four. 

" The 



We show the vortexes and vortex lines in Fig. 20 



left and middle panels of Fig. 20 show the lines that re- 
main in the equatorial plane, along with color-intensity 
plots depicting the lines' vorticities. 



The left panel of Fig. 20 shows the counterclockwise 
lines and their vorticities. As in the case of the (2,2) 
mode, we see a spiraling region of strong vorticity which 
contacts the horizon, and an accompanying spiral of low 
vorticity. At the horizon, the frame-drag field is primar- 
ily longitudinal-transverse, and correspondingly its vor- 
tex lines enter the horizon at a (nearly) 45° angle. As 
for the (2,2) mode, there is a limiting spiral that all the 
outspiraling vortex lines approach, near the edge of the 
vortex. 

There is also a small region of strong vorticity near 
the horizon which forms a second spiral, opposite the 
primary spiral, although it quickly becomes weak; this 



second vortex coincides with the region of strong positive 
vorticity, as we see in the middle panel of Fig. 20 , and 



we think its existence is due to the frame-drag field at 
the horizon being primarily longitudinal-transverse. It 
also should be compared to the similar regions of strong 
negative vorticity near the positive horizon vortexes of 



the magnetic-parity (2, 2) mode in Fig. 10 



In the middle panel, we plot both the counterclockwise 
(dashed) and clockwise (solid) vortex lines, and we color 
each point by the vorticity that is strongest. We see 
two strong vortexes spiraling out to form gravitational 
waves, and we see that under a rotation through 180° 
the clockwise and counterclockwise vortex lines map into 
each other. 



Finally, in the right panel of Fig. [20j we show the vor- 
texes in three dimensions using the same conventions as 
the red (light gray) and blue (dark gray) sur- 



in Fig. 14 



faces are the locations where the vorticity of largest mag- 
nitude has fallen to 90% of its maximum at each radius. 
By contrast with the (2,1) magnetic-parity mode, where 
the 3D vortexes are antisymmetric through the equator 
(and so they flip colors; see Fig. 19), here they are sym- 
metric and so have the same color above and below the 
equatorial plane. 

By duality, for the (2,1) magnetic-parity mode, with its 
antisymmetric 3D vortexes (Fig. 19), the 3D perturbative 
tendexes are symmetric through the equatorial plane and 
have the form shown in this right panel of Fig. 20 
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FIG. 20: (color online). The (2, 1) electric-parity vortex lines, vorticities and vortexes in the equatorial plane. For this mode 
the horizon's vorticity vanishes, so the horizon is plotted as a white disk or sphere. Left panel: Counterclockwise vortex lines 
(dashed) and their vorticity (red [light gray] color) normalized as in Fig. 10 and plotted in a region 24M across. Middle 
panel: Both clockwise (solid) and counterclockwise (dashed) vortex lines, and the vorticity (color) of the line with the larger 
magnitude of vorticity, in a region 24M across. Right Panel: Three-dimensional vortexes colored and labeled as in Fig. |14| 
By near duality, this figure also represents to good accuracy the perturbative tendex lines, tendicities and tendexes for the 
magnetic-parity (2,1) mode. 



C. Vortexes of (2,0) magnetic- parity mode and 
perturbative tendexes of (2,0) electric-parity mode 



In Sec. |I C 8[ we described in detail the dynam- 
ics of the axisymmetric (2,0) magnetic-parity mode of 
Schwarzschild and the gravitational waves it emits — 
waves in which the vortex and tendex lines wrap around 
deformed tori. In this section and the next, we shall dis- 
cuss some other details of this mode and its dual, the 
(2,0) electric-parity mode. 

In Fig. [2l] we show the horizon vorticity for this 
magnetic-parity mode. Of course, it is proportional to 
the scalar spherical harmonic Y 20 (9, (/)). At this moment 
of time, there are clockwise vortexes (blue [dark gray]) 
in the northern and southern hemispheres, and a band- 
shaped counterclockwise vortex (red [light gray]) in the 
equatorial region. As time passes, the horizon vorticity 
oscillates, with red vortexes becoming blue and blue be- 
coming red in each half cycle, while also decaying expo- 
nent iallv^The cause of these oscillations, as we discussed 
in Sec. I C 8 is exchange of energy between c)B L (whose 



normal-normal component is the horizon vorticity) and 
SS LT (which we will visualize in the next section). 

As we discussed in Sec. |I C 8[ symmetries dictate that 
this mode have two families of vortex lines lying in planes 
of constant <\> and a third family consisting of az- 
imuthal circles of constant r and 6. In Fig.[6j we explored 
in detail the wave-zone wrap-around-torus shapes of the 
vortex lines, and their vorticity patterns. In the near 
zone, the line shapes and vorticities are somewhat more 
complex. We elucidate them in Fig. [22l where, to make 




FIG. 21: (color online). The horizon vorticities (5Bnn) of 
the quadrupolar, (2,0), magnetic-parity mode. As in Fig. 
[9] the color intensity is proportional to the magnitude of the 
vorticity with blue (dark gray) for positive and red (light gray) 
for negative. The arrow points along the polar axis. The 
vorticity oscillates sinusoidally in time, causing SBnn first to 
vanish and then to change sign. 



the figure more understandable and preserve some fea- 
tures lost in Fig. [6| we show the two families of vortex 
lines in separate panels, left and center. 

As for the superposed (2,2) and (2, —2) perturbations 
of Sec. [rV] each family takes on both positive and neg- 
ative vorticities, but is predominantly one or the other. 
And unlike the (2,2) mode and the superposed mode, 
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FIG. 22: (color online). Left and middle panels: Near- and transition-zone vortex lines and their vorticities in an <S</> plane of 
constant 0, for the axisymmetric (2,0) magnetic-parity mode of Schwarzschild in the near and transition zones (30M across). 
Left panel: The predominantly counterclockwise family of vortex lines. Middle panel: The predominantly clockwise family of 
vortex lines. Right panel: Vorticity of the axial lines normal to an plane, in near and wave zones (56M across). The color 
intensity in each panel gives the vorticity of the lines, scaled as in Fig. [16] 



the (2,0) line families do not map into each other after 
a 90° rotation; rather, they have distinct patterns (as 
one might expect, since their plane is rather than 
the equatorial plane). On the other hand, because of 
this mode's oscillating nature, the predominantly pos- 
itive lines are the same as the predominantly negative 
lines a half-cycle previous (with signs reversed). For this 
reason, we illustrate the two families at the same moment 
in time, the moment when the horizon vorticity reaches 
a maximum with blue (dark gray) near the poles and red 
(light gray) near the equator. 

One striking feature of Fig. [22] is a set of isolated points 
where six lines meet, three from each family (three in 
each panel). These are nodes (zeros) of the frame-drag 
field, as one can see from the fact that the coloring there 
is white. These are also points where, dynamically, the 
field lines can reconnect, changing their topologies. 

Let us focus on the near-horizon, predominantly nega- 



tive vortex lines (dashed lines) in the left panel of Fig. 22 
The lines that emerge from the counterclockwise horizon 
vortex in the equatorial region loop over the north or 
south pole of the black hole, and reconnect to the oppo- 
site side of that horizon vortex. We think that, as the 
mode oscillates, these lines will merge at the equator then 
slide off the horizon and form closed loops surrounding 
the hole, of the sort that we see in the outer parts of 
the lenticular blue (dark gray) region of the center panel, 
and these will then expand and deform and reconnect to 
form the set of wrap-around-deformed-torus lines of the 
left panel, which lie in the outer part of the transition 
zone and are becoming outgoing gravitational waves. 

Next focus on the near-horizon, predominantly positive 



22 



vortex lines (solid lines) in the middle panel of Fig 
The lines, that emerge from the clockwise horizon vortex 
in the north polar region, swing around the equator and 
descend into the south polar horizon vortex. We think 
that, as the oscillation proceeds, these lines will slide off 
the horizon and immediately form closed loops that wrap 



around deformed tori, which expand to become like those 
near the left and right edges of the left panel (outer part 
of transition zone), and then continue their expansion, 
becoming the gravitational- wave wrap-around-torus lines 
whose inner parts are at the left and right edges of the 
right panel. 

Notice, in the middle panel near the equator, two re- 
gions of weakly negative (pink [light gray]) vorticity, and 
their near-zone lines that appear to have just discon- 
nected from the horizon but are mostly radially directed. 
And notice similarly the pink regions near the left and 
right edges of this panel, again with vortex lines that are 
traveling roughly radially. These pink regions are actu- 
ally toroidal, because of the rotation symmetry around 
the vertical axis. In the outer transition zone and the 
wave zone, they are the regions in which this family's 
wrap- around- tor us, gravitational- wave vortex lines are 
crossing over from one clockwise vortex (wave crest) to 
another. This feature of crossover lines with weakly re- 
versed vorticity appears to be a robust feature of os- 
cillatory modes. For other examples, see the weakly 
blue regions in the left panel, and see the superposed 
(2,2) and (2,-2) mode in Fig. 16 , where the dashed vor- 



tex lines, with predominantly counterclockwise (red [light 
gray]) vorticity, become weakly blue (darker gray) in the 
crossover regions. 



The right panel of Fig [22] shows the vorticity of the 
axial lines (constant r, 9 circles) in both near zone and 
wave zone. Near the horizon, these lines are largely part 
of the transverse, isotropic piece of the longitudinal field 
SB L ; they have opposite color to the horizon vortexes at 
the horizon, as they must, in order to keep 5B L trace- 
free. Near the horizon, these lines also contain a smaller 
component of the ingoing-wave transverse-traceless field 
8B TT . In the wave zone, they are fully out going- wave 



SB 



TT 
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FIG. 23: (color online). For the (2,0) electric-parity mode 
of Schwarzschild in RWZ gauge: the vorticity of the coun- 
terclockwise vortex lines that pass through the plane S</> of 
constant 0. The color intensity (scale on left) is scaled as in 
Fig. [16] Left panel: 30M across, showing the near and inter- 
mediate zones and beginning of the wave zone. Right panel: 
56M across. 



D. Vortex lines of (2,0) electric-parity mode and 
perturbative tendex lines of (2,0) magnetic-parity 
mode 

For the (2,0) electric-parity mode of Schwarzschild in 
RWZ gauge, the only nonzero components of the frame- 
drag field are 5B ^ and SB^. Near the horizon, where 
decomposition into longitudinal, longitudinal-transverse, 
and transverse-traceless parts is meaningful, SB L van- 
ishes (and hence the horizon vorticity vanishes), SB-i is 

the sole component of £B LT , and SB^ is the sole com- 
ponent of 5B TT . 

By (near) duality, the same is true for the (2,0) 
magnetic-parity mode of the last subsection, with SB re- 
placed by S£. 

Because the only nonzero components are SB-^ and 
SBq^ and because of the axisymmetry, there is a family 
of zero-vorticity vortex lines which lie in a plane of 
constant </>, and the other two sets of vortex lines have 
equal and opposite vorticity and pass through at 45 



degree angles. In Fig. 23 , we show in the vorticity of 



the counterclockwise lines that pass through it. A plot 
for the clockwise lines would be identical, but with blue 
changed into red. 

Notice the remarkable absence of structure in the near 
zone. All we see is toroidal vortexes separated by circular 
null surfaces and a polar null line. (Recall the axisym- 
metry around the vertical polar axis). The absence of 
structure is presumably due to the fact that this mode 
is sourced by the longitudinal perturbative tendex field, 
and not by this frame-drag field (though its longitudinal- 
transverse part plays a key role of periodically storing 
near-zone energy during the oscillations; cf. the discus- 
sion of the dual mode below). The vorticity vanishes 
along the polar axis because of axisymmetry and the 
fact that the radial-radial part of the frame-drag field 
vanishes. 

For greater insight into this frame-drag field, we show 



FIG. 24: (color online). Vortex lines for the (2,0) electric- 
parity mode of Schwarzschild in RWZ gauge. Left panel: Two 
vortex lines of each sign in the near zone and innermost part 
of transition zone. The positive (blue lines) and negative (red) 
lines are identical but wind their tori in opposite directions. 
Right panel: Two counterclockwise vortex lines in the transi- 
tion zone, with the vorticity shown in a semitransparent slice 
of constant as a density plot. The vortex lines are plotted 
in black rather than red in this panel to aid the eye. 



vortex lines in the near zone and innermost part of the 
transition zone. These vortex lines wind densely around 
axisymmetric deformed tori. 

Note that the large torus is attached to the horizon. Its 
vortex lines intersect the horizon at the approximately 
45° angles characteristic of the longitudinal-transverse 
part of the field, which is these lines' dominant compo- 
nent. 

For the dual, magnetic-parity (2,0) mode, this torus 
depicts the perturbative tendex lines of the near zone, 
and those lines predominantly belong to the longitudinal- 
transverse part of the tidal field, SS LT . This is the part 
that stores the mode's near-zone oscillation energy when 
SB L is passing through zero and its perturbative vortex 
lines are detached from the horizon (see the discussion of 
this mode's dynamics in Sec. IC8). Immediately after 



this snapshot, these tendex lines' 58 begin to regen- 
erate the near-horizon longitudinal frame-drag field SB L 
and its horizon vorticity. As it does so, these tendex lines 
and their torus (presumably) detach from the horizon 
and expand outward into the transition then wave zone, 
becoming the tendex-line component of a gravitational- 
wave torus like those displayed in Fig. [7] above. 

The small torus in the left panel of Fig. [24] encircles the 
equatorial point on the innermost node of the field (in- 
nermost white circle in Fig. 23). It is also the innermost 
torus shown in Fig. [7] above. 



in the left panel of Fig. 24 several of its three-dimensional 



In the right-hand panel of Fig. 24, for the (2,0) electric- 
parity mode we show two counterclockwise vortex lines 
in the transition zone. By chance, the larger of the se- 
lected lines nearly forms closed orbits, and so even after 
wrapping its torus four times it appears as a thin ribbon. 
While it is difficult to tell with this nearly closed line, ex- 
tended integration reveals that it does wrap a (deformed) 
torus. 

Notice that both lines (both tori) in the right panel 
straddle the second null of the frame-drag field (second- 
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from-center white circle in Fig. 23). 



VI. CONCLUSIONS 

Although the theory of black-hole quasinormal modes 
is roughly half a century old, most past studies of them 
have focused on their mathematical properties, their 
eigenfrequencies and emitted gravitational waves, and 
their excitation by various physical processes. Aside 
from a geometric-optics interpretation of high-frequency 
modes (see, e.g., [22-24 ), little was known, before this 
paper, about their geometrodynamic properties — e.g., 
the structure and dynamics of their spacetime curvature 
in the near zone and transition zone, and how the near- 
zone curvature generates gravitational waves. 

In this paper we have used our new vortex and ten- 
dex tools to explore, in depth, the geometrodynamics of 
the quadrupolar modes for Schwarzschild and Kerr black 
holes. Most importantly, we have discovered that: 

(i) There is a near duality between electric-parity 
modes and magnetic-parity modes, in which tendexes get 
mapped into vortexes and vortexes into tendexes. 

(ii) The electric-parity (/,ra) = (2,2) and (2,1) modes 
are generated by near-zone, longitudinal vortexes that 
extend out of the horizon and rotate (four tendexes for 
m = 2; two for m = 1). The vortexes' rotation gener- 
ates outgoing and ingoing gravitational waves; the ingo- 
ing waves act back on the longitudinal vortexes, gradually 
pushing them off the horizon, which results in the mode's 
exponential attenuation. 

(iv) By (near) duality, the electric-parity (2,2) and 

(2.1) modes are generated and attenuated in the same 
way, but with near-zone, longitudinal tendexes rather 
than vortexes playing the central role. 

(v) The magnetic-parity (2,0) mode and superposed 

(2.2) & (2,-2) mode are generated by near-zone, lon- 
gitudinal vortexes that extend out of the horizon and 
oscillate between clockwise and counterclockwise vor- 
ticity. In these oscillations, energy is fed back and 
forth between the longitudinal vortexes and longitudinal- 
transverse, near-zone tendexes that do not penetrate the 
horizon. In each oscillation, as the horizon vorticity 
passes through zero, the longitudinal vortex lines slide 
off the horizon and reconnect to form toroidal vortexes 
that travel outward, becoming gravitational waves; and 
the near-zone tendexes then regenerate the longitudinal 
vortexes (with reversed vorticity), thereby triggering the 
next half cycle of oscillation. 

(vi) The electric-parity (2,0) mode and superposed 
(2,2) & (2,-2) mode exhibit these same geometrodynam- 
ics, but with the roles of the vortexes and tendexes re- 
versed. 

In future papers, these quasinormal-mode insights will 
be a foundation as we explore the geometrodynamics of 
merging binary black holes using numerical simulations. 

While all analytic approximations fail near the time 
of merger, black-hole perturbation theory does approxi- 



mate a binary-black- hole spacetime well in some epochs: 
the merged hole during its ringdown, each tidally de- 
formed hole during inspiral, and each perturbed hole 
during the initial relaxation that causes spurious "junk" 
gravitational radiation. Before exploring the fully non- 
linear vortex and tendex structures in simulations, we 
are likely first to compare numerical vortex and tendex 
structures during these epochs with the corresponding 
perturbative results (as in Fig. 15). Such comparisons 



will allow us to determine to what degree the insights 
we have gained from our perturbative studies can also be 
applied to numerical simulations — particularly the rela- 
tive insensitivity of vortexes and tendexes to changes in 
gauge and slicing. 

Building on these comparisons, our future work will 
then include initial explorations of the fully nonlinear 
geometrodynamics of the warped spacetimes present in 
binary-black- hole simulations. For example, Kamaret- 
sos, Hannam, and Sathyaprakash [25] have recently ob- 
served relationships between the properties (masses and 
spins) of the initial holes in a binary-black-hole merger 
and the particular quasinormal modes that are excited 
in the remnant (modes which generate the ringdown por- 
tion of the gravitational waves). By examining the vor- 
tex and tendex structures of a variety of binary-black- 
hole mergers, we hope to gain insight into the origin of 
such relationships. Also, following Dennison and Baum- 
garte's recent exploration [7 of the vortex and tendex 
structures in approximate, perturbative initial data, we 
intend to explore the vortex and tendex structures of 
constraint-satisfying binary-black-hole initial data, which 
could give insight into the initial perturbations (and 
the corresponding spurious "junk" gravitational radia- 
tion) that appear in all currently used, binary-black-hole 
initial-data schemes. 

Ultimately, we plan to use vortexes and tendexes to ex- 
plore the geometrodynamics of binary-black-hole space- 
times throughout the entire simulated inspiral, merger, 
and ringdown. We expect that these tools will provide in- 
sights into the behavior of these spacetimes and perhaps 
also motivate new ways of constructing phenomenologi- 
cal waveform templates for use in gravitational- wave data 
analysis. 
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Appendix A: Quasinormal Modes of a Schwarzschild 
Black Hole in Regge- Wheeler Gauge 

In this appendix, we review the Regge- Wheeler-Zerilli 
(RWZ) formalism for black-hole perturbations, and we 
discuss the calculations that underlie the results reported 
in Sees. |III[ [TV] and [V] for quadrupolar perturbations of 
non-spinning black holes in the RWZ gauge. 



1. Regge- Wheeler-Zerilli formalism 

Here we review the equations governing quasinor- 
mal modes for a non-rotating black hole in the Regge- 
Wheeler-Zerilli gauge [27[[28]. 4 We write the metric in 
Schwarzschild coordinates with a small perturbation h^, 



ds 2 



-a 2 dt 2 



dr 2 



a 2 = 1 - 2M/r . 



- r 2 (d0 2 + sin 2 Odcf) + h^dx^dx" , 

(Al) 



The components of h^ v obey separable differential equa- 
tions, and importantly h^ v can be split into definite- 
parity perturbations (electric and magnetic) which do not 
couple to each other. 

For magnetic-parity perturbations, the only nonzero 
components of h^ v in Regge- Wheeler gauge are 



h tA = h {r)e- mt XT{e, 4>) , h rA = h^e-^X^i 

(A2) 

Here uo is the mode's complex QNM eigenfrequency, and 
X 1 ™ is the magnetic-parity vector spherical harmonic on 
the unit two- sphere, 



- esc 0Y l 



X 1 ™ = sm0Y lm f , (A3) 



with Y lrn (0,(j)) denoting the scalar spherical harmonics. 
Regge and Wheeler [27j [29] showed that the radial parts 
of the metric perturbation, ho(r) and fti(r), can be ex- 
pressed in terms of a single scalar radial eigenfunction 
Q(r) as 



/>o = - — (rQ), r , fti = ^, (A4) 



which satisfies the eigenequation 
Q,r*r* WQ = V Q (r)Q , V Q (r) = a 2 < 1(1 + ]) 



Here r* is the tortoise coordinate 
dr 



2Mln(a 2 r/2M) , 



6Af\ 
73" J ■ 

(A5) 



(A6) 



which goes to +00 far from the hole and —00 at the hole's 
horizon. This eigenequation must be solved subject to 
the boundary conditions of outgoing waves at infinity, 
Q ~ e lujr * as r* — >• +00, and ingoing waves at the horizon, 
Q ~ e~ lujr * as — )► —00. 

For electric-parity modes, the nonzero components of 
hftv in RWZ gauge are [28] 

h tt = a 2 H (r)e- i " t Y lrn 
h tr = H^e-^Y 1 ™ 



h r 



Ho(r) 



a* 

„2 ( 



h AB = r 2 n AB K(r)e- iujt Y lrn . 

(A7) 



Here &ab denotes the metric on the unit 2-sphere. We 
can write the metric perturbation functions in terms of 
the Zerilli function Z(r) as 5 



K 
Hi 
Ho 



A(A + l)r 2 + 3XMr + 6M : 



21 



-iuj 



r 2 (Ar + 3M) 
Ar 2 - 3AMr - 3M 



Z + aZ q 



2 1 



(r-2M)(Ar + 3M) 



Z — iuorZ , 



Ar(r - 2M) - u 2 r 4 + M(r - 3M ) 



(r-2M)(Ar + 3M) 
(A + l)M-o; 2 r 31 



K 



Hi • 



iujr(Xr + 3M) 
Here we have used Zerilli 's notation 



A=-(Z-l)(Z + 2). 



(A8) 



(A9) 



(A10) 



The Zerilli function satisfies the eigenequation 

Z, r ^r* + u 2 Z = V z (r)Z , 
where 

2 T2A 2 (A + l)r 3 + 6A 2 Mr 2 + 18AM 2 r + 18M 3 ' 
Vz ^~ a { r 3 (Ar + 3M) 2 

(AH)" 

The slices of constant Schwarzschild time t do not in- 
tersect the black hole's horizon, so in performing our 
3+1 split, we use slices of constant Eddington-Finklestein 



There are many errors in the original paper of Regge and 
Wheeler [27], most of which were corrected by Edelstein and 
Vishveshwara [29]. We use the corrected equations without fur- 
ther comment. 



5 In Hi we have corrected a term in the numerator of the fraction: 
the last term, — 3M 2 , was incorrectly written as — 3M by Zerilli, 
an error that should be obvious on dimensional grounds. 
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time t = t + 2Mln(r/2M - 1). Written in Schwarzschild 
coordinates, t he pe rturbe d te trad for the EF observers is 
given by Eqs. ( |2.4| ) and ( |2.7| ). For any chosen mode, we 
compute the frame-drag and tidal fields by (i) computing, 
from the metric-perturbation components h^ vi the per- 
turbation 5R a p 7 s to the Riemann tensor (same as Weyl) 
in Schwarzschild coordinates; (ii) projecting the total 

Riemann tensor R a (3 7 5 = Ra^s+^R^s (where R^ s is 
the unperturbed Riemann tensor) onto the perturbed EF 
tetrad; (iii) reading off S &i = R &m and % = \e m R m 
and splitting them into their unperturbed and perturbed 
parts. 



2. Magnetic- parity (2,m) mode: Frame-drag field 

We first focus on the (2, m) quadrupolar modes for 
magnetic-parity perturbations. Carrying out the above 
computation, expressing the answer for the frame-drag 



field in terms of the Regge- Wheeler function Q(r) and 
the electric-parity scalar, vector, and tensor harmonics 
(see discussion in App. C3), and simplifying using Eq. 
(A5), we obtain: 



rr 


= SR 


-Bl(m) 


rA 


= SR 


B 2 (m) 


AB 


= 5R 


B 



-iujtyr2m 



; 5 i(m)^ y2m 



2(m) 



(r) = 



12Q 

Q 1 

luor 6 

UMujQ + 2a 2 rQ' 
iujr 3 a 2 y/l + 2M/r 



(A12a) 
(A12b) 



(A12c) 
(A12d) 

(A12e) 



B 



3(m)( r ) = - ■ 5 Aj—^n i [ 3 « 2 ( r - M )( r " + 4M ") + 4iMwr 2 (r - 3M) - rV(r 2 + 4M 2 )] Q 
iujr°or\r + 2M ) V 

+ ra 2 [(r - 3M)(r 2 + 4M 2 ) + 4iMcjr 3 ] Q') , (A12f) 



where a prime denotes a derivative with respect to r, 



Fj m and Y^g are given by Eqs. flC23| ), and is the 
Kronecker delta. 



We have solved the Regge- Wheeler equation (A5) nu- 



merically for the most slowly damped, quadrupolar nor- 
mal mode. When the numerical solution is inserted into 
the above expressions for , numerical errors cause 

ab 

problems with delicate cancellations in the transverse- 
traceless and radial-transverse components near the hori- 
zon. To deal with this, we have derived the following 
asymptotic formula for Q(r) near the horizon, r*/M <C 
-1: 



Q 



Y 



-2iMoj 



3Y 



(1 - AiMuj)e 
9iMujY 2 
"(1 -4iMu)(l -2iMu)e 2 

3(1 + 12iMu + 40M 2 uj 2 )Y 3 
' 2(1 - MMuj)(1 - 2Miuj)(3 - 4iMu)e 3 



0(Y 4 ) 
(A13) 



ities for the fundamental (2,2) quasinormal mode. We 
illustrate them in Figs. [2j [3| [2] and [lO} 

For our superposition of the (2, 2) and (2, —2) modes 
we can simply sum the (2, 2) and (2,-2) harmonics in 
the above expressions. We plot the vortex lines for the 
resulting frame-drag field in Figs. [8j [l6j and the top row 
of Fig. 18 We use the (2,1) harmonics for generating 
the vortexes of the magnetic-parity, (2,1) perturbations 
that are illustrated in Fig. [l9] Finally, we use the (2,0) 
harmonics to produce the vortexes and vortex lines of 
the (2, 0) magnetic-parity perturbation. We note that 

y?° = YJii = for this mode - This means that B[V is 

4> Ucp ab 

block-diagonal, and the vortex lines split into a pair of 
lines which remains in a slice of constant (j) and a single, 
axial line that runs in circles of constant (r, 6). In a 
slice of constant 0, we illustrate the two sets of vortex 
lines in the slice and their voticity together in Fig. [6] and 



separately in the left and middle panels of Fig. 22 We 



also plot the vorticity of the axial lines in a slice in the 
right panel of Fig. [22) 



where Y 



r*/2M 



Inserting this into Eqs. (A12), we 



find, of course, that all components of B^ are finite at 
the horizon. 

Using Eqs. ( A12| ) for the frame-drag field, our analytic 
formula (A13) for Q(r) near the horizon, and our nu- 
merical solution for Q(r) at larger radii, and the (2,2) 
harmonics, we compute the vortex lines and their vortic- 



3. Electric-parity (2,m) modes: Frame- drag field 

Carrying out the calculation described at the end of 
Sec. |A 1| using the electric-parity metric perturbation 
( A7 ) , expressing the result in terms of the Zerilli function 



Z with the aid of Eqs. ( A8 ), and simplifying using Zerilli's 
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differential equation (A10), we obtain for the frame-drag where 



field of a (2,m) electric-parity perturbation 



rr 



(i) 



o, 



B (1 ) = 5R [B^e'^Xf 1 ] 



fA 

B { }\ 

AD 



-iujt y2ml 

x ab\ 



(A14a) 
(A14b) 
(A14c) 



(e) 



[6M 2 a 2 



5. 



2(e) 



: (2r + 3M)] Z - 2Mra 2 (2r + 3M)Z' 



2r 5 a 2 v /l + 2M A 



(A14d) 



1 



2r 4 a 4 (2r + 3M)(r + 2M) 



( [-12Ma 2 (M 2 + 4r/3i) - iu(r 2 + 4M 2 )/3 2 + 4Mcj 2 r 3 (2r + 3M)] Z 



ra 2 [4M/3 2 + iur(2r + 3M)(r 2 + 4M 2 )] Z' 
I 



(A14e) 



Here a prime denotes a derivative with respect to r, and 
X 1 ™ and X l £g are the magnetic-parity vector and tensor 
spherical harmonics given by Eqs. ( |C26| | ( see discussion 
in App. C3). We have defined here for convenience the 



functions 



Pi 



r 2 + Mr + M 2 
2r + 3M 



h = (2r 2 -6Mr-3M 2 ). 

(A14f) 

We note again that the horizon vorticity, B$ N , van- 
ishes. With this B^ we can again compute the eigen- 
vector fields and eigenvalues for the perturbed spacetime, 
and from them compute the vortex lines. We use these 
expressions to calculate the vortex lines and their vortic- 
ities generated by electric-parity perturbations. In order 
to compute the vortex lines for these modes, once again 
we expand Z around the horizon in terms of Y = e r */ 2M 
up to 0(Y 3 ). We use this series to match to a numerical 
solution of the Zerilli equation subject to ingoing-wave 
boundary conditions. Because the Zerilli potential V z is 
more complicated than the Regge- Wheeler potential Vq , 
the coefficients of the expansion of Z in powers of Y are 
lengthy, but easily computed using algebraic computing 
software such as Mathematica. For this reason, we do 
not give the coefficients here. 

For an electric-parity (2,2) perturbation, the only set 
of vortex lines that are confined to the equatorial plane 
have vanishing vorticity (and are of less physical interest). 
Instead, we used the above frame-drag field to compute, 



and then plot in Fig. 13, the vorticity of one of the sets 
of vortex lines that pass through the equatorial plane at 
a 45 degree angle: the set with negative vorticity. Just 
as with the magnetic-parity modes, we superpose a (2, 2) 
perturbation with a (2, —2) perturbation by a simple sum 
of the harmonics. The vorticity of these lines passing 
through the equatorial plane (the analog of Fig. 13) is 



plotted in Fig. [17 

For the vortex lines of the (2, 1) mode, there is a re- 
flection symmetry about the equatorial plane, which im- 
plies that there are two sets of vortex lines confined to 
the plane, with a third normal to it. We illustrate the 
vortex lines and 3D vortexes of this mode in Fig. |20j Fi- 
nally, when we use the (2,0) harmonics, we note that 
X T = x Tz = x fi = 0- While this means that the 

v 00 <p<p 

frame drag field is simple, it is not block-diagonal and its 
nonzero vortex lines pass through all three dimensions. 
There is a single set of axial vortex lines with zero vortic- 
ity, and two sets with equal and opposite vorticity that 
wind around deformed tori. We illustrate the vorticity 
in a slice of constant <fi for the negative lines in Fig. [23[ 
In addition, we illustrate some 3D vortex lines in Figs. [7] 
and [2! 



4. Electric-parity (2, 2) mode: Tidal field 

To help understand the slicing dependence of our 
results, we compare fields generated by electric-parity 
perturbations, because the slicings are identical for all 
magnetic-parity perturbations. In particular, we focused 
on the perturbed tidal field for the electric-parity, (2, 2) 
mode. Carrying out the calculation of this mode as abov e 
when using the electric-parity metric perturbation (A7), 
we obtain 
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where Z is a function which obeys the same Regge- 
Wheeler equation (A5) as Q and can be built from the 



Zerilli function Z as, (see, e.g., Ch. 4, Eq. (156) of [30] ) 



A 2 (A + 1) 
3M 



3Ma 2 



r(Ar + 3M)_ 



A 



(A16) 



for integers I > 2. This implies there is an exact du- 
ality between Eff for electric-parity perturbations and 
Brr magnetic-parity perturbations [in fact for any (Z,ra) 
mode] in RWZ gauge. This follows from the facts that 
these radial-radial components have the same time, ra- 
dial, and angular dependence (but not necessarily the 
same amplitude and phase). However, we can fix the 
relative normalization of the Regge- Wheeler function Q 
and Zerilli function Z such t hat Q = —o jZ/8, in which 
case we have for Eqs. (A12d) and ( |A15a[ ) 



(A17) 



Substituting Z into Eqs. (A15b) — ( |A15c ) does not illus- 
trate the near-duality of the other components of £^ 
and in an obvious manner, so we leave these equa- 
tions in terms of Z. 

As we discuss in the next section, however, the exact 
duality of sQ) and B$ does not immediately correspond 
to an exact duality of the horizon tendicity and vortic- 
ity. This happens because in RWZ gauge, the electric- 
parity perturbations deform the horizon, which changes 
the horizon tendicity. 

For the electric-parity, (2, 2) perturbation, the tidal 
field is symmetric about the equatorial plane, and there 
are two sets of tendex lines that remain in the equato- 
rial plane (just as the vortex lines of the (2, 2) magnetic- 
parity mode did). The tendex lines are illustrated in the 
left-hand panel of Fig. 11 



5. 



Perturbed horizon and horizon tendicity for 
electric-parity modes 



We discuss here the correction to the position of the 
horizon and its influence on the perturbed horizon ten- 
dicity for the electric-parity (2, 2) modes. First, we cal- 
culate the correction to the horizon position Sr using the 



same procedure as that of Vega, Poisson, and Massey 
[18] . The horizon generators, Z, for the perturbed space- 
time are given by 



(A18) 



with an overdot represents a derivative with respect to t. 
The functions SO and S<j) change the location of individ- 
ual generators, but do not alter the shape of the surface 
defined by the instantaneous horizon. We will not treat 
them here, but they are described in [18 . By requir- 
ing that the generators remain null to first order in the 
perturbation, we find 



Sr - 4MSf = 2Mhu . 



(A19) 



For I R ga uge, hu 
Eq. (IA19I) is Sr = 



= and the only physical solution of 
0. Magnetic-parity RWZ perturba- 
tions also have hu = 0, and, therefore, the coordinate 
location of the horizon does not change in this gauge ei- 
ther. For electric-parity perturbations in RWZ gauge, 
we use the fact that hu = hu on the horizon to solve for 
the perturbation to the horizon's shape. For a general 
electric-parity perturbation of indices (Z,m), Eqs. (A7) 
allow us to write 



Sr = 3? 



-2iMu 



1 lm 



lim 

^2M 



r-2M 



(A20) 



where k = (4M) _1 is the horizon's unperturbed surface 
gravity. We evaluate these quantities on the horizon us- 
ing the near- horizon expansion of the Zerilli function Z, 
and they are finite. 

The perturbation to the position of the horizon cor- 
rects the perturbative horizon tendicity in two ways: 
First, the background horizon tendicity £Q\ when eval- 



uated at r = 2M 
Sr, 

AO) 



Sr, becomes, through first order in 



S^\r = r H ) 



1 



8M 3 



Sr . 



(A21) 



4M 2 

Next, we recall that = 4 
of constant r through perturbative order. Now that the 



ej 1 ^ is normal to surfaces 
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horizon's surface is deformed, however, the normal to 
the horizon N is no longer precisely the same as . It 
receives a correction such that 

N» = ^ v V v (r + Sr) 

=ef^ + + 5N" - (snJ^) e<°>" , (A22) 

where N = + N^ is a normalization factor and 
SN^ = (^V u 6r)/N^ deforms N away from . Note 

that the leading-order normal remains 

deformation of the horizon normal produces additional 

modifications to the horizon tendicity, 



£nn =£^{2M + 8r)N»N» 
AO) , cW , 3 



-2£ 



8M 3 



Sr 



2£^SN» 



(A23) 



where, as usual, sQ) includes the effects of both the per- 
turbation to the tidal field and to (and where all 
quantities are evaluated at the unperturbed horizon po- 
sition r = 2M). The new contributi ons [t he last three 
terms on the right-hand side of Eq. (A23)] come from 
the displacement of the position of the horizon Sr and 
the deformation to the normal SN. 

In RWZ gauge, the (£, r) components of SN vanish, 
although SN does have angular components; this means 
that the deformation to the normal to the horizon SN 
does not affect the horizon tendicity in RWZ gauge. [To 
show this, note first that when the deformation to the 



normal has no (t, r) components SN v e, 



(O)i/ 



0. Then 



observe that the (projected) spatial tidal field £^I is di- 
agonal and that SN has only angular components; there- 

and all terms involving SN in 
Eq. (A23) vanish as well.] Only the shifted coordinate lo- 
cation of the horizon, changes the horizon tendicity, and 
we find 



fore, the term S^SN^ 



(1) 
rr 



8M 3 



Sr . 



(A24) 



From Eq. (A20), we see that the angular distribution of 
in RWZ gauge is the same as in IR gauge [it is 



£ 



r 22 (#, </>)]. 

With the angular dependence of the horizon tendicity 
well understood, let us focus on the amplitude and time 
dependence of the horizon tendi city. Using a notation 
analogous to that in Eq. (A15a), we write the horizon 
tendicity in the form 
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Plot illustrating the contributions to the ampli- 
[Eq. ( |A25 )] of the perturbed horizon tendicity 



in RWZ gauge for the electric-parity, (2,2) perturbation. 
Plotted agianst t are the amplitude contributions from £^) 
(dashed line), and from the perturbative shift of the horizon 
generators [dot-dashed line; see Eqs. (A21) and ( |A24 )]. The 
time-dependent amplitude of the total perturbed horizon ten- 
dicity is the solid line. 



(and phase) of £^\ and the other from the correction to 
the radial perturbation of t he gen erators [second term on 
the right hand side of Eq. (A24)]. We plot these contri- 
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as a func- 



butions to E 1{e) e~ lUJt of Eq. (|A25| in Fig. 
tion of t and normalized by the maximum of the (per- 
turbed) horizon tendicity. We also plot the amplitude 
of total perturbation to the horizon tendicity, Ei^e~ lu}t 
(the sum of the two contributions). The two contribu- 
tions are of roughly the same magnitude, but are out 
of phase. The influence of the change in horizon posi- 
tion (dot-dashed line) is slightly larger than sQ) (dashed 
line). 

That £jp N differs from £iV only by an amplitude and 
phase means that, in some sense, the duality between 
the horizon tendicity and vorticity (which is exact in IR 
gauge) is still intact; however, they are no longer related 
by the simple phase shift of i. In fact, we could choose a 
different normalization between the Regge- Wheeler func- 
tion Q and the Zerilli function Z than we did in App. |A4] 
to restore this duality relation, but this would only hold 
for t he ho rizon tendicity and vorticity [and the duality in 
Eq. (A17) would be more complicated, with a complex 



amplitude replacing the factor of i] 



Magnetic-parity, superposed (2, 2) and (2, —2) 
modes: Tidal field 



4 = B 1(e) (r = 2M)e- <wi V 22 



(A25) 



for some am plitud e ^i( e ) [we can do this because both 
terms in Eq. ( A24 ) have the same time dependence] . This 
amplitude has two contributions: one from the amplitude 



In order to understand better the dynamics and in- 
teraction of the tendexes and vortexes, we compute the 
perturbations to the tidal field that must accompany the 
frame-drag field of a magnetic-parity perturbation for the 
superposed (2, 2) and (2, —2) modes. These tidal pertur- 
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bations are much like the frame-drag perturbations of an 
electric-parity metric perturbation, as expected by the 
near-duality. Their odd parity ensures that they must be 
reflection antisymmetric about the equatorial plane. If 
we consider the tendex lines of £^ alone, there must be 
a zero tendicity set of tendex lines in the plane, and two 
sets which pass through the plane at 45° with equal and 
opposite tendicity. We are also assured that £ 



(i) 

NN 



0. 

From the calcultion described at the beginning of Sec. 
A 1 1 above, we obtain: 
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We illustrate the tendicity of the predominantly negative 
tendex lines in the equatorial plane in the time series of 
Fig. [l8] (bottom row) , which shows the evolution over a 
half period of oscillation of the metric perturbation. 



Appendix B: Teukolsky's Equation and Black-Hole 
Perturbations in the Newman-Penrose Formalism 



When a is taken to zero, we recover the Kinnersley tetrad 
for Schwarzschild spacetime in the Schwarzschild coordi- 
nate basis [Eq. (4.1) of Paper II]. The Teukolsky equation 
also requires the NP spin coefficients, certain contrac- 
tions of covariant derivatives of the tetrad above given by 
Eq. (4.1a) of [32] (though with the opposite signs because 
of differing metric-signature conventions). The nonzero 
spin coefficients in this tetrad are 



The results in this appendix appear in many places in 
the literature (see, for example, Teukolsky's paper [31]). 
We summarize them here because we will need them in 
Apps.[C]and[E] 

Teukolsky's equation relies on the Newman-Penrose 
(NP) formalism using Kinnersley's tetrad, which is the 
principal complex null tetrad in the Schwarzschild and 
Kerr spacetimes. For Kerr, in the Boyer-Lindquist coor- 
dinate basis {<9t, d ri do, d^} [Eq. (6.1) of Paper II], this 
tetrad's contravariant components are given by 

^(r 2 + a 2 ,A,0,a), 
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21! 



(r 2 +a 2 ,-A,0,a) : 
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rw 



V2(r + ia cos 9) 



{ia sin 0,0,1,2 esc $) , (Bl) 



with the final leg given by m*, the complex conjugate of 
fa. Here 

A = r 2 + 2Mr + a 2 , E = r 2 + a 2 cos 2 . (B2) 
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a = 7r — (3* , 

r - M 

7 = " + ^^' 
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(B3) 



The Weyl scalars ^0 an d ^4 are defined in terms of 
the Weyl tensor by ^ = C fIl/p J^m u l p m a and ^4 = 
C fIupcr n fI m* u n p m* cr . These both vanish in the back- 
ground when using the Kinnersley tetrad, and are gauge 
invariant at first order in the perturbation theory [31], 
consequently. At that perturbative order, they satisfy 
decoupled, linear, second-order partial-differential equa- 
tions. Teukolsky's big breakthrough [31 was to show 
that, when those equations are re-expressed in terms of 



ip2 = ^0 and ip-2 = p ^4 



(B4) 
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they take on a unified form (the Teukolsky equation) that 
depends on the spin- weight s — +2 for ^2 and s = —2 for 
V?_2, and that is separable; i.e., it has a solution of the 
form ^ = sRimU^sSimU^e^-^. The Teukol- 
sky equation implies for the radial function s Rimuj the 
following ordinary differential equation (in vacuum) 



dr 



s-\-l d s Rlrnuj 



dr 



K 2 - 2is(r - M)K 



+4isujr + 2amuj 



2 2 
a uj 



?Ri 



(B5) 



where s Ai m is a separation constant that is a function 
of au [i.e., s Ai m = 3 Ai m (au)], and K = (r 2 + a 2 )u - 
am. The radial function has the symmetries s Ri m uj = 
(-1)™ sRi-m-u* • The angular function, _ 2 Sj ma; (0), 
(called the spin-weighted spheroidal harmonic) satisfies 
the ordinary differential equation (in vacuum) 



= csce-Une— 



+ (a uj cos 



— 2aa;5 cos — 2ms cot esc — s 2 cot 2 



- m 2 esc 2 
(B6) 



This angular function has the symmetries 
sSi m M ~ 0) = (-lY m+l ^- s S lmuj (0) and s Sf TOW (0) 



(-1) 



s^l—m- 



«(0), where we are using a phase con- 



vention such that the angular functions agree with the 
usual convention for spin-weighted spherical harmonics 
in the limit that spin parameter, a, goes to zero. 

It is often useful, in working with the perturbation 
equations, to change variables from the separation con- 
stants s Ai m to Chandrasekhar's [30] 



s Xlm = s Ai m + S + \s\ 



2amuj + a 2 uo 2 



(B7) 



which are the same for positive and negative spin weights, 
±s. 



specific perturbations of the Weyl curvature scalars ^0 
and ^4 in the Schwarzschild spacetime. Ori then derived 
a similar result for Kerr black holes, using a frequency- 
domain calculation [35]. Yunes and Gonzalez were the 
first to explicitly compute the metric of a perturbed Kerr 
black hole from the Hertz potential [36 , and Keidl, Fried- 
man, and Wiseman were the first to use the procedure 
to calculate the metric perturbation from a static point 
particle in the Schwarzschild spacetime [37] . More re- 
cently, Keidl, Shah, and their collaborators articulated 
a formalism for computing the gravitational self-force 
of a point particle in the Schwarzschild or Kerr space- 
times using the metric constructed from a Hertz poten- 
tial [38]. They were then able to compute the conserva- 
tive piece of the self-force from this metric perturbation 
in the Schwarzschild spacetime [39]. In the first article 
[38] , they gave a concise summary of constructing met- 
ric perturbations from a Hertz potential, and they called 
this process the Chrzanowski-Cohen-Kegeles (CCK) pro- 
cedure or formalism (names we will also adopt). 



In the first part of this appendix, we will review the 
CCK formalism in a similar manner to how Keidl sum- 
marized it in [38 . While the metric we ultimately com- 
pute in this paper is nearly identical to that described by 
Chrzanowski [33 , we find it helpful to put Chrzanowski's 
original calculation into the context of the more recent 
work on the CCK procedure. Furthermore, we review 
the CCK procedure here, rather than simply referring 
the interested reader to [38], because there are several 
differences between our calculation and that set forth in 
[38] : we use a metric of the opposite signature, we cal- 
culate the metric corresponding to quasinormal modes 
with complex frequencies, we construct the metric in a 
different radiation gauge, and (like Chrzanowski's origi- 
nal calculation) we are interested in metric perturbations 
of definite parities. 



Appendix C: The Chrzanowski-Cohen-Kegeles 
Procedure and the Ingoing-Radiation-Gauge Metric 

In this appendix, we will review the formalism used 
for computing the ingoing-radiation-gauge (IR gauge) 
metric, using what is known as the Chrzanowski-Cohen- 
Kegeles (CCK) procedure. We will also connect the 
CCK procedure to Chrzanowski's original calculation of 
definite-parity harmonics, which we find useful for our 
calculations. 

Although Chrzanowski conjectured that "the concep- 
tual benefits of having found the perturbed Kerr met- 
ric potentials surpass the usefulness of these potentials 
for doing future computations" [33], the procedure he 
helped to formulate has found several applications in 
the past few years. Lousto and Whiting [34] revisited 
Chrzanowski's construction and found explicit expres- 
sions for computing the Hertz potential corresponding to 



Because the CCK formalism relies heavily on the 
Newman- Penrose formalism and Teukolsky's equation for 
perturbations of Weyl curvature scalars, we review these 
in Appendix [B] In the second part of this appendix, we 
will describe how to use the CCK procedure to com- 
pute definte-parity metric perturbations corresponding 
to quasinormal modes. In the third part of this appendix, 
we compute the metric perturbations in a notation in 
which they can be compared more easily with those of 
the RWZ formalism (a calculation originally performed 
by [33]), and we also give explicit analytial expressions 
for the tidal and frame-drag fields for (2,2) perturbations, 
which highlight a near duality between the perturbative 
pieces of these fields for perturbations of opposite pari- 
ties. In the final part, we summarize how we numerically 
calculate the IR gauge metric perturbations that we use 
in the visualizations in Figs. [4] [12] [14] [15] and 26 



41 



1. The CCK procedure 



The purpose of the CCK procedure is to construct 
a metric perturbation, h^ Vl from a given solution to 
Teukolsky's equation, either ip2 — or ^-2 = P _4 ^4 
(see Appendix [B] for a summary of the Teukolsky for- 
malism). As part of the calculation, it is necessary to 
relate the solutions of the Teukolsky equation to a Hertz 
potential from which the m etric perturbation is directly 
cons truct ed [see Eq. (C5) for the general relationship, 
Eq. (Cll) for the relationship for the radial functions for 



their harmonics, and Eq. (C18) for the relationship of the 



radial functions of definite-parity perturbations]. 

The CCK procedure can construct a metric in either 
ingoing-radiation (IR) gauge, 



or outgoing-radiation gauge 



h^n v 



0. 



\1V 



0. 



(CI) 



(C2) 



for Schwarzschild and Kerr black holes. Here V and n v 
are two vectors of a Newman-Penrose null tetrad [f or ou r 
calculations, we will use the Kinnersley tetrad, Eq. ( |B1[ )], 
and is the background Schwarzschild or Kerr metric. 
Because our goal is to compute vacuum perturbations 
of Kerr that are regular on the future event horizon, we 
will construct the metric perturbation in IR gauge, and 
we will be able to compute it by algebraically inverting 
a differential relationship between the harmonics of the 
Hertz potential and those of [the result is in Eq. 

den} ]- 

The Hertz potential is tensor with the same symme- 
tries as the Riemann tensor, whose double coordinate 
divergence is a harmonic coordinate metric. Stewart [40] 
showed that in Type D spacetimes, there is sufficient 
gauge freedom that one can represent the independent 
degrees of freedom of the perturbative part of the Hertz 
potential as a single complex scalar; furthermore, if one 
applies a coordinate transformation from harmonic gauge 
into IR gauge, the Hertz potential, which we will denote 
by ^h-, is a solution of the vacuum Teukolsky equation 
for scalars of spin weight s = — 2 (the same as ^-2)- 
One can then construct a metric perturbation from the 
Hertz potential by applying several differential operators 
to 

V = {-Z M k(« + a* + 30-T)(5 + 4/3 + 3T) 
-m^m v (T) - p + 3e - e*)(D + 3p + 4e) 
+Z (/i ra I/) [(D + p* - p + e* + 3e)(<5 + 4(3 + 3r) 
+(<S + 3(3 - a* - - r)(D + 3p + 4e)]} 9 H 
+c.c. (C3) 



the complex conjugate of the entire expression, so that 
the metric perturbation is real. 

When computing perturbations of black holes, it is 
helpful to be able to relate a given Hertz potential ^ h 
to a specific perturbation of the Weyl scalar ^4. It is 
possible to do this by computing the components of the 
perturbative Riemann tensor from the metric perturba- 
tions (C3 that correspond to the Weyl scalar 



^4 



The result can be expressed compactly as 



1 



^-2 = 5 (£ f4 ^ - 12Md t * H ) 



(C4) 



(C5) 



see, 



e.g., Table I of the paper by Keidl [38], where our 
& is their C). We have used the shorthand that jC^ 4 = 
jC^_ 1 CqC\C2 where 

C\ = — (do + s cot — i esc Od^) + ia sin 0d t . (C6) 

In general, solving for the Hertz potential that corre- 
sponds to a perturbed Weyl scalar ^4 involve s inv erting 
the fourth-order partial differential equation (C5); how- 



ever, when and ijj-2 are expanded in harmonics in the 
frequency domain, it is possible to perform the inversion 
algebraically. 

The algebraic inversion can be completed by expanding 
in harmonics, 



Imuj 



Imoj 



where - 2 ^ mu; (r) and _ 2 ^ mw ^) satisfy Eqs. ([B5] 



(C7) 
and 

(|B6|. In the IR gauge, the Hertz potential is a solution to 
the Teukolsky equation with spin s = — 2; consequently, 
it can also be expanded in the same harmonics 



H 



/ H 



,i(m0— tot) 



Imuj 



(C8) 

The radial function of the Hertz potential's harmonics 
a l so satisfies the vacuum Teukolsky radial 
equation, but because it is not the same radial function 
as in the harmonics of ^-2, we denote it with a different 
function. The radial functions of the harmonics of ^ h 
and ip-2 can be related by substituting Eqs. (C7) and 



identity 



(C8) into Eq. (C5) and using the Teukolsky- St arobinsky 



2^lmuj 



D-2Si ri 



(C9) 



(Eq. (59) of Ch. 9 of [30], after noting that our C\ is 
equivalent to — C s there), and the identity 



-iSimu - (-l) m 2^_ 



(CIO) 



(see, e.g., Eqs. (93) and (94) of [37]). The differential 
operators are defi ned by D = and 6 = m^V^. The 

last term in Eq. (C3), denoted by "c.c," means to take 



Then, it is necessary to equate the full radial function 
for a given angular and time harmonic of the Hertz po- 
tential to the radial functions _2^ mw of ^-2- After this 
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relationship is inverted, the individual radial harmonics 
of the Hertz potential can be written as 



-2^h 



D* 2 + 1UM 2 uj 2 



(Cll) 



The constant D* D z * mcj Di_ m _ u * is defined by 

L> 2 = A 2 (A + 2) 2 -8A(5A + 6) (a 2 cj 2 - aracj) + 96Aa 2 cj 2 
+144(a 2 cj 2 - amw) 2 , (C12) 



where A is the separation constant (B7) used by Chan- 
drasekhar [30 (a choice of the separation constant that is 
the same for both the angular and the radial equations). 
Although the Teukolsky-Starobinsky identities are usu- 
ally derived assuming real frequencies, they have been 
shown to hold for complex frequencies as well (for a re- 
cent derivation, see [41]). 

The general description of the CCK formalism is now 
complete: (i) For a vacuum perturbation of ^4, we can 
find the Hertz potential that corresponds t o th is per- 
turbation by expanding \I> h in harmonics [Eq. ( |C8[ )], and 
expressing the radial functions of th is exp ansion in terms 



of those of -0_ 2 = p \I>4 via Eq. (Cll). (ii) From the 



resulting we c an then compute the metric perturba- 
tions via Eq. (C3). 

Next, we will show that if we choose the radial function 
of the Hertz potential to correspond to a perturbation of 
with definite parity, then the result of this calcula- 
tion is equivalent to Chrzanowski's original calculation of 
definite-parity metric perturbations [33] . 



2. Definite-parity harmonics and Chrzanowski's 
calculation 



Let us now connect this CCK procedure with 
Chrzanowski's original calculation of metric perturba- 
tions of definite parities. We shall begin by showing 
that, if the perturbations have definite parity (electric 
or magnetic), then the Hertz potential must itself trans- 
form as —> ±(— 1)^# under parity (the plus and 
minus correspond to electric- and magnetic-parity per- 
turbations, respectively), and its radial harmonics must 



satisfy _ 2 X* 



l—m- 



±(-iy 



iXimuj. In turn, this im- 



plies that the radial harmonics of must satisfy the 
same relationship _ 2 ii*_ m _ w * = ±(-l) m _ 2 ^ m w 

To deduce these relationships, we d iscuss the parity 
of the terms that appear in Eq. ( |C3| ). The Newman- 
Penrose tetrad and spin coefficients of the background 
spacetime transform in several different ways under par- 
ity: I and ft have positive parity, and fa does not have 
a definite parity, m — » — m*. Similarly, the differential 
operator D = has positive parity, and S = m M V M 

again does not have a definite parity, S — > — 5*. Three of 
the nonzero spin coefficients map to their complex con- 
jugates under parity (p — » p*, \i — >• /i*, and 7 — >• 7*), and 
the remaining four spin coefficients become minus their 



complex conjugates under parity (a — )► —a*, (3 — > — /?*, 
7r — » — 7r*, and r — > — r*). These relationships hold true 
for both Schwarzschild and Kerr, although in the former 
case, the spin coefficients are real and, therefore, have 
definite parity. 

When applying a parity transformation to the pertur- 
bative metric tensor, h^ v dx^dx v where h^ v is given by 



Eq. (C3), we can show that the tensor differential opera- 



tor in Eq. (C3) becomes its complex conjugate by using 
the parity transformations for the spin coefficients, NP 
tetrad, and differential operators above. As a result, the 
metric perturbation will have either electric or magnetic 
parity when the Hertz potential transforms as 



(C13) 



under parity. The plus sign corresponds to an electric- 
parity perturbation, and the minus sign describes a 
magnetic-parity perturbation. The condition this implies 
on the harmonics is also quite simple, which we can de- 
termine by applying a parity transfor matio n to the Hertz 
potential expanded in harmonics [Eq. ( |C8[ )] and equating 
it to its complex conjugate. Then using the properties of 
the Teukolsky angular functions 



Or-0) = (-1) 

(-1) 



ra+Z 



s^l — m—uj* ($) 



(C14) 
(C15) 



(see Appendix [B]) and equating the radial function of 
each time and angular harmonic, we obtain the following 
condition on its radial functions, 



- 1- 



±(-1) 



(ci6) 



Similarly by substituting Eq. (Cll) into the expression 



above, we find an analogous relationship for the radial 
function of the Weyl scalar ^_2, 



-2^1 



±(-1) 



2^mw • 



(C17) 



For these definite-parity perturbations, the relationship 
between the radial functions of the Hertz potential and 
^_2, Eq. ( Cll| ) also simplifies, 



±8(L>* ± UiMuj)- 1 _ 2 Ri 



(C18) 



namely, for definite-parity perturbations, the radial func- 
tions of and ^ h differ by only a complex constant. 
Because Eq. (C18) shows that the two radial functions 
_2^ mw and -2Rimuj differ only by a constant multiple, 
we will express both and ip-2 in terms of the radial 
function of -0_2, -2^ mw , for simplicity. 

In the next part (and also for all other IR gauge calcu- 
lations in this paper), we will compute a metric pertur- 
bation that corresponds to a perturbation of of the 
form 



= ±i(D* ±12iMcu)- 2 R lmiV e i ^-^ . 2 S h 

+ -(-l) m (D t 12iMLJ*)- 2 Rtmuje~ i{m ' t> ~ OJ ' t) 

(C19) 
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The corresponding Hertz potential is 

\& H = „ R, p i(m<t>-wt) 



-2Ri m ue^- ut > - 2 Simu (C20) 
±(-l) m _ 2 i?r row2 e- i ( m ^* f )_ 2 5 i _ m _ w ,. 



We choose the prefactors on the modes of ip-2 so as to 
make the Hertz potential (and, therefore, the metric) as 
simple as possible. Furthermore, this choice gives the 
same definite-parity metric as that of Chrzanowski (when 
we take the real part of his expressions). 



3. Definite-parity CCK metric perturbations and 
tidal and frame-drag fields for Schwarzschild black 
holes 

In the first two parts of this section, we will cal- 
culate electric- and magnetic-parity perturbations of 
Schwarzschild black holes in IR gauge. Because 
Chrzanowski performed this calculation in Table III of 
reference [33 , and our results agree with his, we do not 
go into great detail describing the calculations; instead, 
we aim show the results here so as to be able to compare 
with the RWZ formalism in Appendix [A] In the third 
part, we will compute the tidal and frame-drag fields 
corresponding to these metric perturbations and show a 
near duality of the tidal and frame-drag fields of opposite 
parity perturbations for the (2, 2) mode. 



a. Electric-parity metric perturbations 

We begin this part by comparing the metric produced 
by the CCK procedure to that of the RWZ formalism. We 
will write the RWZ metric using the covariant notation 
described by Martel and Poisson [42] . Martel and Poisson 
write the electric-parity perturbations as 



k( e ) _ h lrn v lrn 
n ab — / . n ab Y •> 



Ira 



Im 



\ -lm\rlm 
Ja 1 B •> 



Im \ 
AB •> 



(C21a) 
(C21b) 
(C21c) 



where the lowercase indices run over the radial and time 
coordinates (e.g., a, b = t,r), and uppercase indices run 
over the angular coordinates as before, A, B = 6, <j>. The 
angular functions Y lm are scalar spherical harmonics, 
Y B m are the electric-parity Regge- Wheeler harmonics, 
and Y^jj are transverse-traceless, electric-parity tensor 
harmonics; the term Qab is the metric on a 2-sphere. 
The vector and tensor harmonics are defined by 



Im 



Y 



1 AB 



rim 



D A Y l 
D A D B 



-i(i + i)n 



AB 



Y 



Im 



(C22a) 
(C22b) 



where D A is the covariant derivative on a 2-sphere. 

Because the Schwarzschild spacetime is spherically 
symmetric, we can see, intuitively, that the CCK 
metric, Eq. ( |C3[ ), corresponding to an electric-parity 
quasinormal-mode perturbation [the plus sign in Eq. 
(C20)] will have a relatively simple form. The angular 
operators acting on the Hertz potential in Eq. (C3) be- 



come the spin- weight raising and lowering operators, and 
the angular functions become the spin-weighted spheri- 
cal harmonics; furthermore, and when the spin- weighted 
harmonics are combined with the appropriate factors of 
m and m* the angular functions become proportional to 
the scalar, vector, and tensor harmonics described above. 
When performing the calculation, we will need to use 
the following identities, which can be found, for exam- 
ple, by adapting Eqs. (2.22a) and (2.38e) in the review 
by Thorne [43] to the notation used here, 



Y 1 / 



1(1 + 1) 



(-iYi m m A - iYi m m* A ) , (C23a) 



D 



Ya% = ^(-2Yi m m A m B + 2 Y lm m* A m* B ) . (C23b) 

The Teukolsky-Starobinsky constant for spin-weighted 
spherical harmonics is D = (I + 2)!/(Z — 2)!. We can 
then find that the metric coefficients are given by 



l (e) 
HA 



(C24a) 



a Ka ~ 21(1 



+ l)a 2 1 



dr* ' 



2 Rim 



7 (e) 
AB 



2a 2 , „ 

H -2Rlm 

r 



d 



y-lm —iujt 



(C24b) 



(iujr 2 - M)-^_ 2 R lm - [ln 2 a 2 



dr 



iw(-3r + IM) - r z w z }- 2 Ri 



y-lm —iujt 
AB e 



(C24c) 



In the last equation we have used the radial Teukolsky 
equation to eliminate the second-derivative term, and we 
have defined ji 2 = (I — l)(l + 2) (which is also equal to 
1(1 + 1) - s(s + l) for s = -2). 

There are a few noteworthy differences between the 
IR gauge electric-parity perturbations, and the electric- 
parity RWZ-gauge metric. The CCK metric has a strictly 
angular part of the perturbation which is proportional 
to the transverse-traceless harmonics, and the trace por- 
tion of the angular block vanishes; conversely, the angular 
block of the RWZ metric perturbation has a trace part, 
but no transverse-traceless perturbation. The hff part 
of the metric perturbation also has a simpler relation- 
ship with the h[f and h r e J components in IR gauge than 
in RWZ gauge; one reason for this is that the IR gauge 
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metric has electric-parity vector perturbations, whereas 
the RWZ metric sets these to zero. Finally, the IR gauge 
metric is finite on the future event horizon for ingoing ra- 
diation. One can see this by noting that both -2 Rim and 
d-2Rim/dr* scale as a 4 e~ zujr * near the horizon, which 
will cancel any negative powers of a 2 in the expressions 
for the metric coefficients. The same is not as manifest 
for the RWZ perturbations (see Appendix [A] for more 
details on the RWZ formalism). 



Because they have the same radial dependence as the 
electric-parity metric, the magnetic-parity peturbations 
will also be well-behaved on the future event horizon. 

The major difference between the RWZ formalism's 
magnetic-parity metric and the IR gauge metric is that 
in IR gauge, the transverse-traceless metric perturbation 
is no longer required to be zero. 

c. Tidal and frame- drag fields of the (2,2) mode 



b. Magnetic-parity metric perturbations 



The magnetic-parity perturbations are given by 



ri ab — U ' 



lm 



7 (m) 
AB 



lm 
AB •> 



(C25a) 
(C25b) 

(C25c) 



lm 



where the magnetic-parity harmonics are defined by 



In this part, we calculate the tidal and frame-drag 
fields for a (2,2) mode in IR gauge of both electric and 
magnetic parities. We find an interesting near duality be- 
tween the tidal and frame-drag fields of opposite-parity 
peturbations that we noted in Sees. |I C 3| |IIIB 3| and 
IIIIC 21 

We compute the tidal and frame-drag fields from the 
metric by evaluating the com ponents of the Weyl tensor 
and its dual in the tetrad ( |2.4|) i ncluding the perturbative 
corrections to the tetrad (|2.7a|)-(2.7d). We find that for 



an electric-parity mode, the tidal and frame-drag fields 
can be written as 



lm 



x 



AB 



-e A B D B Y 
-\{e A c D B + e B c D A )D c Y l 



(C26a) 
(C26b) 



££' e) =23?[£ I(e) (r)y^e— «] 



and cab is the Levi-Civita tensor on a unit 2-sphere. 
As in th e previous part, we can compute the CCK met- 
ric ( |C3| ), which is relativitely simple for a Schwarzschild 
black hole. The reason for the simplification is the same, 
but we will need the following two identities that relate 
the spin- weighted spherical harmonics to magnetic-parity 
vector and tensor harmonics 



Xlm 
A 



Xlm __ 
AB - 



(-iYi m m A + iYi m m* A ) , (C27a) 



D 



-i—(-2Yimm A m B - 2Yi m m* A m* B ) . (C27b) 



These relationships can found in Eqs. (2.22b) and (2.38f) 
of [43]. The magnetic-parity metric perturbations have 
the same radial and time dependence as the electric- 
parity perturbations for the vector and tensor parts, 



_ h (m) _ 2 h (m) _ VP 
~ a ^A ~ 2l{l + 1} 



2oS 
r 



d 

dr* 



iuj ■ 



lR\m 



2 Rim 



Xlm ^ — iujt 
A 6 



(C28a) 



7 (m) 
AB 



• iu(-3r + 7M) 



(iujr 2 - M) — _ 2 Rim ~ \\^ 2 o? 



r 2 uj 2 }- 2 Rim 



AB e 



(C28b) 



Eff =2K[E u{e) {r)Yfe-^} 



rA 

c(l,e) 
C AB 



o(l,e) 

rA 
AB 



25R 



-^E I(e) (r)5 AB Y 



22 



^II(e)WF^ ]e 



AB 



25R[B I(e) (r)Xfe— *] , 
2®[B II(e) (r)Xf 6 e- i " t ] 



(C29a) 
(C29b) 



(C29c) 

(C29d) 
(C29e) 
(C29f) 



The symbol S^ B is the Kronecker delta function, and the 
traceless property of £ requires that the radial function 
in front of the Kronecker delta must be minus one-half 
that of £$ [i.e., -(l/2)£ 1(e) (r)]. 

For the magnetic-parity perturbation, the frame-drag 
and tidal fields are 



o(l,m) 

ff 
B (lm) 

rA 



B 



(l,m) 
AB 



(l,m) 



rA 
AB 



: 23?[S II(m) (r)F| 2 e-- t ] 
-m[{-\B l{m) {r)5 AB Y 
B m{m) {r)Yf B )e-^] , 

■-2U[E l{m) {r)X 2 ?e- iM }, 
■■m[E ll{m) {r)Xf B e- iM ) 



22 



(C30a) 
(C30b) 

(C30c) 
(C30d) 
(C30e) 
(C30f) 



Interestingly, the radial functions of the tidal and 
frame-drag fields of the opposite-parity perturbations are 
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nearly identical 

Bl(m)(r) = ^l(e)O) 
£ll(m)0) = ^II(e)(0 - 1 



(C31a) 



M^/3(r + 2M) 



r 5 a 4 v // 2r 



6/ 



-(2a + iuor)- 2 R 22 + r-—- 2 R 22 



(C31b) 



£i(m)0) = -i5 I(e) (r) +2 



5 a 4 v / 2 



M v / 3(r + 2M) 



(C31c) 



r 5 a 4 y2r 



-(2a 2 + iujr)_ 2 R 22 + r-—_ 2 R 22 

(C31d) 

^ii (m) (r) = -iSii ( e)(r) - ^T7^ Mw - 2i?22 ' 

(C31e) 



In fact, there is an exact duality of the radial-radial com- 
ponents, which implies that the horizon vorticity of a 
magnetic-parity perturbation is the same as the horizon 
tendicity of an electric-parity perturbation. For com- 
pleteness, we list the expressions for the radial functions 
for the electric-parity perturbations, which are lengthy, 
but will be needed in the next appendix. 



E l(e) (r) = (r 2 (r - 3M + iujr 2 )-^_ 2 R 22 + [-5r 2 + WMr - 12M 2 - iujr 2 (4r - 9M) + rW}_ 2 R 22 } , 

T Oi ^ CIV* J 

(C32a) 

E u(e) (r) = - 1 r -(r 2 [3r 2 + 6M 2 + iwr 2 (r - 3M) + r W}-^-. 2 R 22 

+ [(-9r 3 + 18Mr 2 - 12M 2 r + 24M 3 ) - iwr 2 (8r 2 - 16Mr + 18M 2 ) + 2wV(4r - 9M) + ir 6 w 3 ]_ 2 i? 22 j , 

(C32b) 

£ III(e) (r) = r5(r + / M)Q4 ^ M-2r 2 + 3Mr + 3M 2 + rW)±_ 2 R 22 + [6(r 2 + 4M 2 ) 

+ iw(4r 3 - HMr 2 + 12M 2 r + 12M 3 ) - w 2 r 2 (4r 2 - 4Mr - 9M 2 ) - 3iu 3 r 5 a 2 + r 6 w 4 ]_ 2 i? 22 j , (C32c) 

B Ue) (r) = ; = -jr 2 [9M -iu>r(r-3M) + r 3 u> 2 }^-- 2 R 22 

I(e)W rV 3r ( r + 2M )" 4 I dr * 

+ [-24Mra 2 + iojr(12M 2 - 25Mr + 5r 2 ) - w 2 r 3 (4r - 9M) - ir 5 w 3 ]_ 2 i? 22 | , (C32d) 
B u{e) (r) = r4(r + 2 ^ )Q4 ^ {^(-2r 2 + 3Mr + 3M 2 + rV) ±. 2 R 22 

+ [-24M + 2icjr(2r - 7M) + cj 2 r(-4r 2 + 4Mr + 9M 2 ) - 3icj 3 r 4 a 2 + r 5 uj 4 ]_ 2 R 22 1 . (C32e) 



I 

From these expressions, it is clear that the tidal and responding powers of a in the denominators of these func- 
frame drag-fields are regular on the horizon, because, as tions. 
noted above -2 Rim an d d- 2 Ri m /dr* scale as a 4 e~ lur * 
near the horizon; consequently, they will cancel the cor- 
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4. Analytical and numerical methods for 
computing metric perturbations and tidal and 
frame-drag fields in IR gauge 

The procedures for calculating the metric perturba- 
tions and their tidal and frame-drag fields are identi- 
cal for Schwazschild and Kerr black holes; however, be- 
cause the analytical expressions for the Newman- Penrose 
quantities, the angular Teukolsky function, and the met- 
ric derived from these mathematical objects are signifi- 
cantly simpler for Schwarzschild black holes, the amount 
of work we can perform analytically differs for rotating 
and non-rotating black holes. Even for Schwarzschild 
black holes, however, we will not be able to compute all 
aspects of the metric perturbation analytically. We cal- 
culate the the least-damped I = 2, m = 2 quasinormal- 
mode frequencies for both Schwarzschild and Kerr black 
holes using the Mathematica notebook associated with 
[13], an implementation of Leaver's method [14]. Simi- 
larly, we compute the radial Teukolsky functions -2Rimu 
corresponding to a quasinormal-mode solution for both 
Schwarzschild and Kerr black holes numerically. We 
compute it in two ways, which give comparable results: 
we solve the boundary-value problem for a quasinormal 
mode solution to the radial Teukolsky equation, Eq. (B5 ), 



using a shooting method, and we compare the result with 
a series solution given by Leaver [14] (as is also done in 
the notebook of [13]). For Kerr black holes, the numerical 
solution requires the angular eigenvalue, s ^im associated 
with the quasinormal mode frequency, which we again 
compute from the implementation of Leaver's method in 

The most significant difference between the calcula- 
tions of quasinormal modes of Schwarzschild and Kerr 
black holes arises from differences in the Teukolsky angu- 
lar function, an d th e angular operators used in comput- 
ing the metric (C3). First, the spin- weighted spheroidal 
harm onics in the expression for the Hertz potential, Eq. 
(C20), reduce to spin- weighted spherical harmonics for 



Schwar zsch ild black holes. Second, the angular operators 
in Eq. (C3) reduce to spin- weight lowering operators, in 



the non-spinning limit. As a result, the metric perturba- 
tion can be expressed, analytically, in terms of electric- or 
magnetic-parity scalar, vector, and tensor spherical har- 
monics of a single for Schwarzschild black holes. For 
perturbations of Kerr black holes, there are not these 
additional simplifications. First, we must calculate the 
spin-weighted spheroidal harmonics numerically, which 
we do using a series solution put forward by Leaver [14] 
(the same method as that implemented in [13 ]). Second, 
the angular operators are no longer the spin-weight low- 
ering operators. The metric perturbation computed from 
these functions, therefore, is not nearly as simple as that 
of the Schwarzschild limit. In fact, for our calculations 
with spinning black holes, we find it easier to work with a 
numerical fit to the analytical expression for the metric. 

Once we calculate the metric perturbation, we con- 
struct the perturbation to the Weyl tensor in the same 



way for both rotating and non-rotating black holes. We 
can then calculate the tetrad components of the tidal 
field, and frame-drag field, using the background 



tetrad in Eq. (2.4 for Schwarzschild holes or Eq. (2.5a 



for Kerr holes and its perturbative corrections in Eqs. 
(2.7a)-(2.7d). From these fields, we can solve the eigen- 



value problem and compute tendex and vortex lines, and 
their corresponding tendicities and vorticities. 



Appendix D: Relationship Between 
Regge-Wheeler-Zerilli and Ingoing-Radiation 
Gauges 

In this appendix, we construct generators of infinites- 
imal coordinate transformations between RWZ and IR 
gauges, for both magnetic- and electric-parity perturba- 
tions of Schwarzschild black holes. 



1. Magnetic-parity gauge transformation 

In this part, we compute the gauge-change generator 
that transforms the magnetic-parity metric in IR gauge 
to the same metric in Regge- Wheeler gauge. We show, 



as noted in Sec. II C that this infinitesimal magnetic- 
parity coordinate transformations does not change the 
time function that specifies the slicing (into surfaces of 
constant t). In addition, perturbative changes of the spa- 
tial coordinates will not alter the coordinate (or tetrad) 
components of the frame-drag field; therefore, the fields 
in both gauges will be equal. 

The calculation that shows these facts is relatively 
straightforward. Regge and Wheeler showed in Eq. (17) 
of [27] that, beginning in any gauge, it is possible to re- 
move the transverse-tracel ess pa rt of the magnetic-parity 
metric pe rturb ation [Eq. ( |C25c ) in the notation used in 
Appendix C 3 by an infinitesimal coordinate transforma- 
tion of the form 



r(m) 



1 



2 ^ 



ftr(o,o,x 



lm\ 



(Dl) 



where X^ m is a magnetic-parity, vector spherical har- 
monic. This follows from the fact that the perturbation 
to the metric transforms under this change of coordinates 
by 



(D2) 



(where | denotes a covariant derivative with respect to the 
background metric, and parenthesis around the indices 
means the expression is symmetrized) and from the defi- 
nition of th e magn etic-parity, transverse-traceless tensor 
harmonics ( |C26b| ). The result can also be found from 
Eqs. (5.5) and^6) of [42]. 

For a multipolar perturbation with indices (/, m) in IR 
gauge, the function - \ £ Zm h l 2 m /2 is given by the radial 
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function in Eq. (C28b) multiplied by e~ iut , and the full 



coordinate transformation vector is therefore 



(iivr 2 - M) — - 2 Rir, 



St 


_ t (m) _ 
Sr 


o, 


SA 


1 


H 



(D3a) 



iuj(-3r + 7M) - rV]_ 2 i^ m 



XTe 



(D3b) 



A short calculation can verify that ft^ and fe r ^ are the 
only nonzero components of the metric after this trans- 
formation (the same as in RWZ gauge) and they are given 

by 



+ l)Dra 4 



\u 2 ra 2 + {iMlo + r 2 ui 2 )] 



x —- 2 Rim + [a 2 (a 2 + iojr) - w 2 r(3r - 7M) 



icu 3 r 3 }- 2 Rim 



'ir A = 3 j 



a e 



(D4a) 



+ l)Da 6 



(r-3M + iwr 2 ) — _ 2j R, 



+ [-a 2 {\n 2 + 3a 2 ) - iw(4r - 9M) + tu 2 r 2 }_ 2 R lm 



a e 



}• 



(D4b) 



It is not immediately apparent, however, that this gauge 
is RWZ gauge, because it is expressed in terms of the 
radial function of ^-2, (-2-Rzm)? rather than the Regge- 
Wheeler function Q. 

To show that this transformation did bring the metric 
into Regge- Wheeler gauge, it is necessary to use the re- 
lationship between Q and _2-Rz m given in, e.g., Eq. (319) 
of Ch. 4 of [3D] 6 



Q 



Va*(D-l*Ma,) { (r " 3M + ^^-^ 
+ \-c?{\[i 2 + 3a 2 ) - iw(4r - 9M) 

+ r 2 w 2 ]_ 2 i? ;m l. 



(D5) 



6 Aside from several differences in notation (the radial function 
used by Chandrasekhar, Z^~\ is related to the Regge- Wheeler 
function by Q = iujZ^~\ and his radial function for the Teukol- 
sky equation, Y_2 is related to that of this paper by _2-Rj m = 
r 3 Y_2), there is one additional subtle point about using this 
equation. This equation is expressed as a relationship between 
Y+2 (proportional to the radial function of \Po) and Z^~\ Be- 
cause the time dependence of is given by e+ tcrt in |30| . then 
the Y+2 there is equivalent to y_£ 2 °f with a time dependence 
given by e~ tujt . In addition, because y_2 satisfies the same equa- 
tion as y+2 5 then this equation is valid for y_2 when a is replaced 
by u. 



Afte r substituting this relationship into Eqs. (A4) and 
(A2) and taking its imaginary part — so that the RWZ 
metric is real and is expressed in terms of -2 Rim — it 
becomes apparent that the transformation brings the IR 
gauge metric into RWZ gauge. 

Because the gauge change from IR to RWZ is generated 
by a strictly spatial £ ( m ) and because B is a strictly first- 
order quantity for perturbatio ns of Schwarzschild holes, 
the frame-drag analog of Eq. (2.15) guarantees that the 



frame-drag field must be identically the same in the two 
gauges: 



IRG 



RW 



(D6) 



This can be confirmed explicitly for the (2,2) mode by 
substituting Eq. (D5) into the tetrad compo nents of the 
RWZ frame-drag field in Eqs. flA12a| )-( |AT2f| ) and finding 
that the y are i dentic al to the IR gauge frame-drag field 
of Eqs. dCgOb, PHY and dCggb- 



When thought of as abstract tensors without reference 
to any coordinate system, it is also the case that the tidal 
fields are equal, 



7IRG 



7 RW 



(D7) 



Because there is a background tidal field, perturbative 
differences in the coordinates enter into the components 
of the ti dal fie ld and the components are no longer equal; 
see Eq. (2.15). Therefore, visualizations of the tendex 



lines or tendicity in the two coordinate systems (when 
the coordinates are drawn as though they were flat) look 
different. 



2. Electric-parity gauge transformation 

In this part, we construct an infinitesimal generator of 
an electric-parity coordinate transformation that brings 
the electric-parity IR gauge to RWZ gauge. The trans- 
formation changes the time function (and hence how we 
define the slicing) in addition to the spatial coordinates. 
This implies that neither the frame-drag fields nor the 
coordinate components of the tidal field will equal in the 
two gauges (but the ti dal fie ld wr itten without coordi- 
nates will be); see Eqs. (2.14) and (2.15), respectively. 

The gauge-change generator that connects the two 
gauges is somewhat more complex for the electric-parity 
perturbations than it was for the magnetic-parity ones. 
The transformation can be found by using Eq. (19) of 
[27] or Eqs. (4.6)-(4.9) of [42]. The general approach 



to find the transformation is to use £ A 



(e) 



to remove the 



transverse-traceless part of the IRG metric, and then use 
£/ e ^ and to annul the transverse metric coefficients. 
After a short calculation, it is possible to express the 
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generator as 



radial Teukolsky function - 2 R 22 , we find 



t ( e ) = 

?< 2/i 2 r 

d 



r 2 a 4 \ 



r 2 [/Ta 2 - 4iuj(M - r 2 



x —- 2 Rim + [/i 3 ra 4 + iujr 2 fi 2 a 2 (fi - 2) 



4rMu;V(3r - 7M) - 4iu 3 r% 2 R, 



im 



s^rlm „ — tout 



t(e) = 



r 3 a 6 { [ 1 



-} 

(D8a) 

V(^ + 2) +4iw(r-3M) 



- 4r 2 cj 2 ] — _ 2 ^ m + {-r/i 2 a 2 (/i + 2)(2ra 2 + ir 2 u) 

- 2r 2 (/i + 2) - 2iur[(fi 2 + /i - 2)r 2 + 2Mr(/i - 8) 
+ 24M 2 ] + 4cj 2 r 3 (4r - 9M) + 4^ 3 r 5 }_ 2 i^ 



t(e) 



'a 4 I 



(D8b) 



(iwr 2 - M)-f-_ 2 i?; m - [i/i 2 a 2 



iw(-3r + 7M) - r 2 w 2 ]_ 2 ^ 



y-lra —iojt 
A & 



}■ 



(D8c) 



where we used Eq. (B5) to reduce second-order radial 



derivatives to first-order ones. 

To confirm that this gauge-change generator does bring 
the IR-gauge metric to the RWZ metric, we again use 
the relation between the Zerilli function and the radial 
Teukolsky function encoded in Eq. (319) of Ch. 4 of [30]: 



2r[iujr 2 (n 2 r + 6M) 



r 2 a 4 (D + l2iMuu)(n 2 r + 6M) 

d 



(/i 2 r 2 - 3/i 2 Mr - 6M 2 )]-— _ 2 i^ m + {a 2 (/i 2 r + 6M) 2 
2(3a 2 + iwr)[/ii 2 r 2 - 3/ii 2 Mr - 6M 2 



iojr 2 (p 2 r + 6M)}}- 2 Rh 



(D9) 



This allows us to confirm that the IR metric was brought 
to RWZ gauge through the transformation vector £ ( e ) . 

With this expression, we can also compare the frame- 
drag fields in the two gauges for the (2, 2) mode. By 
expressing the radial functions for the frame-drag field 
of a (2,2), electric-p arity m ode i n RWZ gauge [Bi( e \(r) 
and B 2 ( e )(r) of Eqs. (A14d) and (A14e)] in terms of the 



£i(e)0) = £i( e )0) 



2r 5 a 4 v"2r(r + 2M) 



-4ra 2 



- 2iu(2r 2 - bMr + 6M 2 ) + 3r 2 u; 2 (r - 3M) 
+ ir 4 uj 3 }_ 2 R 22 - r 2 a 2 (r 2 u 2 + 3zMcj - 2) 



2^22 



dr : 

B 2 (e)(r) = B u{e) (r), 



(DlOa) 
(DlOb) 



where Bi( e \(r) and i?n( e )(0 are the equivalent radia l 
functions of the IR-gauge frame-drag field in Eqs. (C32). 



Because the functions £>i( e )(0 and Bj^{r) determine the 
radial dependence of the transverse part of the frame- 
drag field (and i?2(e)( r ) and £?n( e )(0 do the same for 
the transverse-traceless par t), we see a particu lar illus- 
tration of the result of Eq. (2.14) of Sec. II C 2 namely, 
a change in slicing from an electric-parity gauge change 
will induce a change in the longitudinal-transverse com- 
ponents of the frame-drag field (but not the longitudinal 
or transverse-traceless parts). 



Appendix E: Horizon Tendicity and Vorticity 
Calculated from the Weyl Scalar ^ 

In this appendix, we modify a calculation by Hartle [6 
to compute the horizon tendicity and vorticity by apply- 
ing differential operators of the background spacetime to 
a perturbation of ^0 = ^2 (where ip 2 satisfies Teukolsky's 
equation; see Appendix |B|) . Using this result, we derive 
the duality between the horizon vorticity and tendicity 
of opposite-parity perturbation mentioned in Sec. |III A[ 
for both Schwarzschild and Kerr black holes. We also re- 
late the horizon quantities to the complex curvature and 
show that they are proportional for Schwarzschild holes 
and differ only by the product of spin coefficients X^cr^ 
for Kerr holes. This proves these claims made in Sec.|IB| 



Constructing a hypersurface-orthogonal tetrad 
on the horizon 



As in Hartle's calculation, we must work in a NP tetrad 
in which the null vector I is tangent to the horizon, n is 
normal to the horizon, and m and its complex conjugate 
lie in the instantaneous horizon (constant v = t + r for 
Hartle, though we will use constant t). This NP tetrad 
must also satisfy additional constraints 



1 ,t ^ 
u = -^( l + n ) 

iV=-L(f-n) 

m = — (e 2 + tes) 
V2 



(Ela) 
(Elk) 
(Elc) 
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[wit h the associated non-null tetrad given by Eqs. (2.5a) 
and ( |2.7| )], which ensure that the slicing vector u associ- 
ated with this NP tetrad is hypersurface-orthogonal on 
the horizon and the spatial basis vectors are tied to our 
coordinate system in the desired way. 

To describe the unperturbed NP tetrad, it is useful 
to first construct Hartle's tetrad, which can be obtained 



from Kinnersley's tetrad (Bl), by a boost followed by a 



null rotation about / (also called class III and I transfor- 
mations, respectively): 



2(r 2 +a 2 ) 



la sin c 



y/2(r + ia cos 0) 



n H 



2(r 2 +a 2 )^ 



ia sin 



A 

ia sin 



n K 



y/2(r + ia cos 6) 



\/2{r + ia cos 6) 
a 2 sm 2 r 



(E2a) 
(E2b) 

(E2c) 



The quantites A and E are defined in Eq. (B2 ). Then, we 



can construct an unperturbed tetrad from Hartle's tetrad 
using the following spin-boost transformation (also called 
class III): 



(0) 



m 



(o) 



(o) 



iV z _1 n H , (E3a) 



where 
N t 



S + 2Mr 



j@ r + ia cos ( 



(E3b) 



2S ' ■ VZ 
One can verify that the resulting orthonormal tetrad 



tr<°> 

-(0) 



1 
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(i(0) + ™(0)) 

1 



1 
1 
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(ra (0 ) ^ 

(™(0) 



0(0) - ™(0)) , 



(0)J> 



(E4a) 
(E4b) 
(E4c) 
(E4d) 



is exactly the ingoing-Kerr tetrad ( 2.5a| ), when evaluated 
on the horizon, though away from the horizon it is not. 

For the NP null tetrad to correspond, on the hori- 
zon, to the hypersurface-orthogonal {{?, e? , e^, e ?} of Eqs. 



(2.5a) and (2.7) via Eqs. (El), we must choose the per- 



turbative corrections to the tetrad to satisfy 



1 



n (i) 



m (1) 



2v^ 
1 

2^2 
1 



[( ft 66 S (0) - ^ffe f (0) ) - 2%e ( o) 
[(^66^(0) + ^rfe 



2 % e 7o) 



2^(0)] , 

(E5a) 

2^(0)]. 
(E5b) 



2^ 



r(0) 



) + (E5c) 



Because eV is normal to surfaces of constant r in slices 
of constant t through perturbative order, we will need to 
choose our gauge so that the coordinate position of the 
horizon does not move from the constant value r = r+. 
Although Poisson [44] has shown that there are a wide 
class of gauges that satisfy this property (horizon-locking 
gauges), for our calculation, we find it convenient to work 
in an ingoing radiation gauge based on the unperturbed 
tetrad vector Z( )> 

V^(o)= > ^V = °- (E6) 
On the horizon, these gauge conditions imply that 



6 oo 



°0r ' 



(E7a) 



and that the null vector I undergoes a perturbative boost, 



*(i) o W 



66H0) • 



(E7b) 



[To derive this, one should split —2h^e^ into a sum 



of two terms —2(h^e^ + ^qa^(o)) anc ^ use ^ ne relation 
in Eq. (El).] In addition to keeping the horizon at a 
constant coordinate position r = r + (see [44]), using this 
gauge condition allows us to calculate the perturbation 
to ^2 in a much simpler way, as we describe in the next 
subsection. 



2. Computing the horizon tendicity and vorticity 
from ^0 



Although the explicit expressions for spin coefficients 
in this tetrad are somewhat lengthy (and, as a result, we 
do not give them here), through a direct calculation we 
can verify that on the horizon 



7 (0) 



Ho) = > e (o) £ 



(E8a) 



Moreover, because in this ingoing radiation gauge the 
perturbation to the vector I can be obtained by apply- 
ing a boost (Class III) transformation to the tetrad, the 
perturbed value of n will also vanish, 



«(i) 



0. 



(E8b) 



From Eq. (310a) of Ch. 1 of [30 (which describes the com- 
ponents of the Riemann tensor in the Newman-Penrose 
formalism), we see that the perturbation to the spin- 
coefficient p satisfies an equation 

D (o)P(i) = 2€( )P(i) , 

where e( ) > 0. If is not zero, then the separated 

solution to this equation, p^ = f(r^6)e~ l ^ UJt ~ rn ^\ im- 
plies the constraint that 2e( ) + iNi(uj — mco> + ) = [here 
uj+ = a/(2Mr + ) is the horizon angular velocity]. This 
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condition is not satisfied for constant frequencies lj, so 
the perturbation to the spin coefficient must vanish: 



P(i) 



0. 



(E8c) 



From these conditions on the spin coefficients, and the 
fact that ^q ^ = = 0, we can write the Bianchi 

identities (see, e.g., Eqs. (321a) and (321b) of Ch. 1 of 
[30]) as 



(D - 2e)tf i = (<5* + 7T - 4a)^ , 



(i) 



(S* + 2tt - 2a)*i , 



(E9a) 
(E9b) 



where we have dropped the superscripts indicating per- 
turbative orders on all differential operators and spin co- 
effients (because they are all background quantities) and 
the Weyl scalars \I/o and (because they are strictly 
perturbative quantities). Note that we do not need the 
term of the form D^)^^ , because on the horizon the dif- 
ferential operator D( X ) contains only time and azimuthal- 

angle derivatives, but the background Weyl scalar 
is only a function of r and 0. 

By applying the differential operator (D — 2e) to the 
second Bianchi identity and using the identity (valid on 
the horizon) that £*D - D£* = (a + /3* - tt)D, we find 
that 



(D - 2e)D^ 



(i) 



: [5* +3(7r-a) -/T)]((T +7T-4a)^ 
(D-2e)(A* ), (E10) 



Using Geroch-Held-Penrose [45] notation, and the equa- 
tion for a component of the Riemann tensor (Eq. (310g) 
of Ch. 1 of [30]) restricted to the horizon 

DA - S*tt = -2eA + + a - /T ) , (Ell) 

we find that 

ft)*^ = (S'S' + 4ttS' + 2tt 2 - Ab)^ • (E12) 

Note that the \I/o here is related to that which satisfies 
Teukolsky's equation in the Kinnersley tetrad by 



at2 2i& a2 
^ = till fL^K = A ^K 

(r 2 +a 2 ) 2 ~ 



(E13) 



Starting from a modal solution for the Kinnersley \I/o 
(denoted by ^/ m£J ), then we see that the corresponding 
perturbation to ^2 is given by 



(i) _ 9'9' + 4ttS' + 2tt 2 + X(iNiQ + 4e) 



2lr, 



NiQ(2ie - N t n) 



l moo J 1 



(E14) 



where ^ = — mco> + , and where we have used the fact 
that D = Niipi + ^ + 9t) on the horizon. 

A Weyl scalar, ^ formed from the superposition of 
modes V$ lmuJ ± (-l) m ^z_ m -u;*> with radial functions 
that obey 2^1-m-uj* — ±( — l) m 2^* mu; 7 transforms un- 
der parity as \Er -> ±( — 1)^0- Tne perturbation of ^ 2 



formed from superimposing Eq. (E14) for the individ- 
ual modes of ^ above also transforms under parity as 
^2 ±(-l) l ^2' Using the relation 2^ 2 = S NN ^iB NN 
and taking the real and imaginary parts of ^2, it becomes 
clear that Enn and Bnn have definite parity. Moreover, 
it is not difficult to see that £nn of an electric-parity 
mode is equal to iBnn of a magnetic-parity mode of 



K 

Imuj ' 



and £nn of a magnetic-parity mode is —iBnn 



of an electric-parity mode of i^^i muj - 

This demonstrates a perfect duality between electric- 
parity modes and magnetic-parity modes, on the horizon 
of a Kerr black hole. 



Relationship between ^2 and the complex 
curvature 



As a final part of this appendix, we discuss how the 
relationship between the complex curvature and ^2, 



1 



(K + iX) = -tf 2 + [ip - Act , 



(E15) 



simplifies for perturbations of Schwarzschild and Kerr 
black holes in the tetrad and gauge discussed in the sec- 
tions of this appendix above. First, we note that the spin 
coefficient A has as its unperturbed value on the horizon 



Mr + uj\ sin 2 ( 



-2i& 



(0) Ni(r-ia cos 9f 

x [4Mr + + (r+ - M) (r+ - ia cos 0)] , (E16) 

where we have made use of the fact that on the horizon 
r++a 2 = 2Mr + . For a Schwarzschild black hole, uj+ 
vanishes, and, therefore, the background values of all four 
spin coefficients /i, p, A, and a all vanish. Through first- 
order in perturbation theory, therefore, 



n = -2£ 



AT AT • 



X = -2B 



NN • 



(E17) 



[We briefly digress here to note that for Schwarzschild 
black hole, the spin coefficient ir also vanishes, and Eq. 



(E12) reduces to 



(i) 



For a modal solution, Eq. (E14) also simplifies to 



(i) 

2 In, 



4Ma 



cj(i- 4MujV2) 



2RlmXlr, 



(E18) 



(E19) 



where a 2 = 1 - 2M/r, D = (I + 2)1/ (I - 2)!, u = tt 
(because uj+ = 0), and the radial function of ^2, 2^Zm 5 
is evaluated at the horizon r = 2M. We have also used 
the fact that e = 1/(8M) in this tetrad. Because the 



spin coefficients vanish in Eq. (E15) for this perturbed 



Schwarzschild hole, the above expression is equivalent to 
minus one quarter of the complex curvature.] 
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For a Kerr black hole A( ) ^ 0, and we must compute 
the perturbation to a. It satisfies the differential equation 



(D - 2e)a (1) = ^ 



(E20) 



[Eq. (310b) of Ch. 1 of [30] specialized to our tetrad and 
gauge]. For a modal solution of a^), we can solve this to 
find 



iNiSl + 2e 



(E21) 



which implies that the perturbation to a does not vanish. 
Thus, for a Kerr black hole, 



1 



(7i(i) + zAf ( i 



-^2^ - A (0)°"(l) , 



(E22) 



so the horizon tendicity and vorticity are no longer ex- 
actly equal the horizon's intrinsic and extrinsic scalar cur- 
vatures. 



Appendix F: Vortex and Tendex Lines of (2,2) 
Perturbations of Schwarzschild and Kerr Black 
Holes with the Background Frame-Drag and Tidal 
Fields 

In this appendix, we show the tendex and vortex lines 
of Schwarzschild and Kerr black holes when we plot a 
small (2, 2) perturbation of either electric or magnetic 
parity on top of the background tidal or frame-drag fields 
in Fig. [26] Specifically, we plot the vortex and tendex 
lines of 



£(0) +6 £(D 



S (0) + e£ (1) , (Fl) 



where £^ and B^ are the stationary, unperturbed 
background fields (visualized in Paper II), £^ and B^ 
are the perturbations (visualized by themselves in Fig. 
14), and e is a constant that sets the scale of the pertur- 
bation. To describe the strength of the perturbation, we 
will compare the perturbative horizon tendicity or vortic- 
ity to the background tendicity (for Schwarzschild holes) 
and the tendicity or vorticity (for Kerr black holes). 
For the Schwarzschild black holes in Fig. [26] we chose 
^nn/^nn " ^ x 10 -4 for the electric-parity perturba- 
tions and SjvTv/^iviv ~ 2 x 10 -4 for the magnetic-parity 
perturbations. For electric-parity perturbations of Kerr 



holes, we chose Z^n/Z^n 



(i) /P (o) 



3.5 xlO- 3 and B%/B K NN 



2(0) 



3 x 10 , and for the magnetic-parity perturbations the 
ratios we selected were S^n/^nn ~ ^.5 x 10 -3 and 
^nn/^nn ~ 5 x 10 -3 . We anticipate that these im- 
ages may be useful for comparing with the results of 
numerical-relativity simulations, for which it is more diffi- 
cult to separate a spacetime into a stationary background 
and dynamical perturbations. 

In the top panels of Fig. [26] are electric-parity pertur- 
bations, and the bottom panels are magnetic-parity per- 
turbations. The left column of images are tendex lines 



of Schwarzschild black holes, the center column are ten- 
dex lines of Kerr black holes of spin a/M = 0.945, and 
the right column are the corresponding vortex lines of 
the perturbed Kerr black holes. The lines are colored by 
the sign of their tendicity or vorticity (blue [dark gray] 
for positive and red [light gray] for negative) and the 
horizons are colored by their tendicity or vorticity. The 
transparent spheres are placed in the figures to help guide 
the eye, and do not indicate any feature of the vortexes 
or tendexes. 

In these figures, we must choose an amplitude for the 
perturbation (described in the first paragraph above). 
For the all the black holes, we make the perturbation 
sufficiently small that one cannot see the effect of the 
perturbation in either the horizon tendicity, or the red 
(light gray) radial tendex lines. For the Kerr holes, we 
also require that the amplitude of the perturbation is 
less than the difference of the tendicities of the two non- 
radial tendex lines at the equatorial plane and around 
the radius at which the angular lines reach closest to the 
horizon. With this choice, the angular tendex lines will 
retain some features of the unperturbed lines before they 
become more distorted by the perturbation in the regions 
near the poles. 

First, we will describe the tendex lines of 
the Schwarzschild black holes. An unperturbed 
Schwarzschild black hole is spherically symmetric, the 
tendicity on a sphere of constant radius is constant, 
and, therefore, any direction tangent to the sphere is a 
valid tendex line. For a weakly perturbed Schwarzschild 
black hole, although the perturbation may be small, 
the perturbation restricted to a sphere of constant r 
completely determines the variation in the tendicity, 
and, furthermore, it will determine the directions of 
the tendex lines. This is analogous to degenerate 
perturbation theory in quantum mechanics, in which the 
eigenstates of the perturbing Hamiltonian restricted to 
the subspace spanned by the degenerate eigenstates are 
treated as the unperturbed states within the degenerate 
subspace. In directions that are not degenerate, however, 
its effects are negligible. 

We can now use these facts about degeneracy to under- 
stand the tendex lines in the angular direction. The tidal 
field in the strictly angular directions, Eq. (C29c) will de- 
termine the structure of the tendex lines on the sphere. 
The angular dependence is determined by the transverse- 
traceless, electric-parity tensor harmonic (for the top-left 
panel), because the trace term in Eq. ( |C29c[ ) is propor- 
tional to the identity and will not lift the degeneracy 
of the tendex lines. We would expect, therefore, that 
the tendex lines in the angular direction would resemble 
those of transverse-traceless, I = 2, m = 2, gravitational 
waves generated by a time- dependent mass quadrupole. 
These were shown in [2j [17] , and the pattern of the lines is 
nearly identical. The tendicity along the lines is quite dif- 
ferent from those of a gravitational wave, because for the 
perturbed Schwarzschild black hole, the tendicity is pri- 
marily determined by the constant unperturbed value on 
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FIG. 26: (color online). Tendex and vortex lines of Schwarzschild and Kerr black holes (of spin a/M = 0.945) perturbed by a 
(2, 2) mode of either electric or magnetic parity, without removing the background tidal or frame-drag fields. The tendex lines 
and vortex lines are colored by the signs of their respective tendicities and vorticities (blue [dark gray] for positive and red [light 
gray] for negative). The horizons are colored and shaded by their vorticities or tendicities, and the transparent spheres have 
no physical significance, but they help to add perspective to the figures. The top panels are electric-parity perturbations and 
the bottom panels are magnetic-parity ones. Left column: Tendex lines of Schwarzschild black holes. Middle column: Tendex 
lines of Kerr black holes. Right column: Vortex lines of Kerr black holes. 



the sphere. Nevertheless, the tendex lines on the sphere 
show a striking similarity to those of gravitational waves 
at infinity. 

For the magnetic-parity perturbation (the bottom-left 
panel), the tendex lines are determined by an I = 2, 
m = 2, magnetic-parity tensor harmonic; consequently, 
we would expect that the lines would resemble those of 
transverse-traceless gravitational waves at infinity, pro- 
duced by a time-dependent, current-quadrupole source. 
Those lines were shown in [17] , and they appear identical. 
Once more, though, the value of the tendicity along the 
lines is set by the background Schwarzschild black hole 
for the lines in bottom-left panel of Fig. 26 (unlike the 
tendicity of the lines studied in [T7]). 

The degeneracy between the angular tendex and vortex 
lines can also be used to explain the tendex and vortex 
lines in the middle and right columns of Fig. |26j respec- 
tively. For both the tendex and vortex lines, when the 
lines are near the equatorial plane (0 = tt/2) they resem- 



ble the unperturbed lines, but as they head toward the 
poles, they begin to become perturbed. This happens 
because the perturbation is small compared to the differ- 
ence in the eigenvalues near the equatorial plane, and the 
perturbations have little effect on the tendex or vortex 
lines. Near the poles, however, the background vorticities 
and tendicities in the angular directions become degen- 
erate (see the discussion at the end of Apps. A and B of 
Paper II), and the perturbation restricted to the degener- 
ate subspace controls the lines' directions. In the vicinity 
of the poles, the degenerate subspace is a plane parallel 
to the equatorial plane, and the perturbative tendex lines 
must form a regular grid around these points. When we 
combine this observation with the parity of the pertur- 
bation, we see that the lines at the opposite poles must 
be either parallel or orthogonal. Thus, these few simple 
constraints combine to explain the relatively simple pat- 
tern of the vortex and tendex lines of the perturbation 
plus the background frame-drag and tidal fields. 
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